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At very high energies or small values of Bjorken x, the density of par- 
tons, per unit transverse area, in hadronic wavefunctions becomes very 
large leading to a saturation of partonic distributions. When the scale 
corresponding to the density per unit transverse area, the saturation 
scale Qs, becomes large {Qs ^ ^qcd), the coupling constant becomes 
weak {ots(Qs) ^ 1) which suggests that the high energy limit of QCD 
may be studied using weak coupling techniques. This simple idea can be 
formalized in an effective theory, the Color Glass Condensate (CGC), 
which describes the behavior of the small x components of the hadronic 
wavefunction in QCD. The Green functions of the theory satisfy Wilso- 
nian renormalization group equations which reduce to the standard hn- 
ear QCD evolution equations in the limit of low parton densities. The 
effective theory has a rich structure that has been explored using ana- 
lytical and numerical techniques. The CGC can be applied to study a 
wide range of high energy scattering experiments from Deep Inelastic 
Scattering at HERA and the proposed Electron Ion Collider (EIC) to 
proton/deuterium-nucleus and nucleus-nucleus experiments at the RHIC 
and LHC colliders. 
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The Color Glass Condensate and High Energy Scattering in QCD 3 

1. Outstanding Phenomenological Questions in High Energy 
QCD 

1.1. Introduction 

The advent of a new generation of high energy colUder experiments, 
beginning with HERA and the Tevatron in the early 1990's to RHIC 
and LHC in the new millenium (and others in the planning stages) have 
launched a new era in the study of the strong interactions. Questions 
which have been around since the early days of the strong interactions, 
such as the behavior of cross-sections at high energies, the universality 
of hadronic interactions at high energies, the nature of multi-particle 
production and the possibility of creating thermalized states of strongly 
interacting matter have acquired fresh vigor. For instance, it is often 
believed that little could be learned about the high energy limit of QCD 
since the physics is assumed to be entirely non-perturbative. On the 
other hand, we have learned from HERA that paxton densities are large 
at high energies or equivalently, at small values of Bjorken x. The large 
densities of small x, or "wee" partons suggest, as we will discuss at length 
later, that semi-hard scales may be present which allow one to describe 
the physics of this regime using weak coupling techniques. QCD at high 
energies can thus be described as a many-body theory of partons which 
are weakly coupled albeit non-perturbativc due to the large number of 
partons. We will call this system a Color Glass Condensate (CGC), for 
the following reasons: 

• "Color" , since the gluons are colored. 

• "Glass" because of the strong analogy of the system to actual glasses. A 
glass is a disordered system which evolves verj' slowly relative to natural 
time scales: it is like a solid on short time scales and like a liquid on 
much longer time scales. Similarly, the partons of interest are disordered 
and evolve in longitudinal momentum in a manner analogous to a glass. 

• "Condensate" because it contains a very high density of massless gluons 
whose momenta are peaked about some characteristic momentum. In- 
creasing the energy forces the gluons to occupy higher momentum states 
(due to repulsive interactions) causing the coupling to become weaker. 
The gluon density saturates at a value of order l/os ^ 1, corresponding 
to a multiparticle Bose condensate state. 

We will argue in the following that the Color Glass Condensate is 
the effective theory describing high energy scattering in QCD. We will 
outline the rich structure of the theory and discuss how it provides in- 
sight into outstanding conceptual issues in QCD at asymptotically high 
energies. The theory can be applied to study phenomena at a number 
of existing and upcoming high energy collider facilities. We will discuss 
applications of the CGC to study the initial conditions and equilibra- 
tion in heavy ion collisions and to describe heavy ion phenomenology. 
We will discuss applications to Deep Inelastic Scattering, to the current 
experiments at HERA but also for future experiments being discussed 



February 1, 2008 20:44 



WSPC/Trim Size: 9in x 6in for Review Volume 



QGPS'final'web 



The Color Glass Condensate and High Energy Scattering in QCD 

at HERA and at Brookhaven. In addition, we will discuss applications 
of the CGC to proton/deuteron- nucleus scattering experiments planned 
at RHIC and at the LHC and to peripheral nucleus-nucleus scattering 
experiments at RHIC. 

1.2. Light cone kinematics and dynamics 

The appropriate kinematics to discuss high energy scattering are hght- 
cone (LC) coordinates. Let z be the longitudinal axis of the collision. 
For an arbitrary 4-vector = (w", u^, v^, ii^) {v^ = Vz, etc.), we define 
its LC coordinates as 

+ _ 1 / , 3n - _ 1 / 3n _ / 1 2^ /TIN 

=~^y'" )' " ^71^^ = )■ (i-i) 

In particular, we shall refer to — {t + z)/\/2 as the LC "time", and 
to x~ = {t — z)/\/2 as the LC "longitudinal coordinate". The invariant 
scalar product of two four-vectors reads: 

p ■ X — p x^ + p^x — pj_ ■ xj_ , (L2) 

which suggests that p~ should be interpreted as the LC energy, and p'^ as 
the (LC) longitudinal momentum. In particular, since p^ = {1/V2)(E± 
Pz) with E = (m^ + P^)^^^, the light cone dispersion relation takes the 
form 



1 p']^+ rr? _ 1 

2 ^ 2 ~ 



(1.3) 



where the transverse mass m_|_ is defined as = p\_ -f rr? . The mo- 
mentum space rapidity is further simply as; 

.-im^^ln^. (1.4) 



2 p- 2 



m 



These definitions are useful, among other reasons, because of their sim- 
ple properties under longitudinal Lorentz boosts: p'^ — > Kp"*", p~ — > 
(1/k)p~, where k is a constant. Under boosts, the rapidity is just shifted 
by a constant: y ^ y + k. 

The utility of light cone kinematics is not merely that of a conve- 
nient coordinate transformation. The Hamiltonian dynamics of quantum 
field theories quantized on the light cone have several remarkable fea- 
tures as was first pointed out by Dirac Firstly, the LC Hamiltonian 
p~ (which is the generator of translations in the light cone time x'^) can 
be written in the form p~ — p^ + V, where p^ is the free Hamiltonian, 
corresponding to a complete set of non-interacting Fock eigenstates, and 
V is the light cone potential. Second, the LC vacuum is trivial, namely, 
the vacuum state is an eigenstate of both the free and the full Hamil- 
tonian. As a consequence of these properties, multi-parton Fock states 
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The Color Glass Condensate and High Energy Scattering in QCD 5 

can be constructed as eigenstates of the QCD Hamiltonian. Thus, in LC 
quantization (and in the Ught-cone gauge — 0) the quark-parton 
picture of QCD becomes manifest. Finally, we note that the apparent 
non-relativistic structure of the light cone Hamiltonian suggested by the 
dispersion relation in Eq. 11.31 is not accidental but is a consequence 
of an exact isomorphism between the Galilean subgroup of the Poincare 
group and the symmetry group of two dimensional quantum mechan- 
ics 1=1 Thus in LC quantization, the Rayleigh-Schrodinger perturbation 
theory with off-shell energy denominators can be used instead of the 
more usual Feynman rules. For a more detailed discussion of the light 
cone formalism and its application to high energy scattering, we direct 
the reader to Ref. 1^. 

1.3. High energy behavior of total cross-sections 

We now return to the outstanding phenomenological questions we men- 
tioned in the introduction. Clearly, computing total cross- sect ions as 
E ^ (X is one of the great unsolved problems of QCD. Unlike processes 
which are computed in perturbation theory, it is not required that any 
energy transfer become large as the total collision energy i5 — > oo. Com- 
puting a total cross-section for hadronic scattering therefore appears to 
be an intrinsically non-perturbative procedure. In the 60's and early 
70's, Regge theory was extensively developed in an attempt to under- 
stand the total cross-section. The results of these analyses were, to our 
mind, inconclusive, and at any rate, certainly cannot be claimed to be 
understood from first principles in QCD. 

On the basis of very general arguments invoking unitarity, analyticity 
and crossing, Froissart has shown that the total cross-sec tion for the 
strong interactions grows at most as fast as In^ i? as i? — > oo^^^^. Several 
questions arise in this regard. Is the coefficient of In^ E universal for all 
hadronic processes? Can this coefficient be computed from first principles 
in QCD? How do we understand the saturation of the unitarity limit 
dynamically in QCD? Or is the Froissart bound an intrinsically non- 
perturbative phenomenon? 

1.4. Multi-particle production in QCD 

Can we compute N{E), the total multiplicity of produced particles as 
a function of energy in QCD ? By this we mean not only the multiplic- 
ity of particles in jets (which is fairly well understood in perturbative 
QCD) but also the total number of particles, at least, for semi-hard 
momenta. Consider the collision of two identical hadrons in the center 
of mass frame, as shown in Fig. The colliding hadrons are ultrarela- 
tivistic and therefore Lorentz contracted in the direction of their motion. 
Furthermore, we assume that the typical transverse momenta of the pro- 
duced particles is large compared to ^qcd ■ We know from experiments 
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The Color Glass Condensate and High Energy Scattering in QCD 

that the leading particles (valence partons) typically lose only some fi- 
nite fraction of their momenta in the collision. The produced particles, 
which are mostly mesons, are produced in the "wake" of the nuclei as 
they pass through each other. In light cone coordinates, the right mov- 



large p small p laige p 

Fig. 1. A hadron-hadron collision. The produced particles are shown as circles. 



ing particle ("the projectile") has a 4-momentum p^* — {p^ , 0^) with 
p^ ~ \/2pz and ~ M^/2\/2pz (since pz ^ M and m± = M, with 
M = the projectile mass). Similarly, for the left moving hadron ("the 
target"), we have = pj^ and p^ = p^. The invariant energy squared 
is s = (pi -|-p2)^ = 2pi ■ p2 — 2pj^P2^ ~ 4p2, and coincides, at it should, 
with the total energy squared [Ei + £2)^ in the center of mass frame. 

Consider a pion produced in this collision and which is moving in 
the positive z direction. For such a pion, we define the longitudinal mo- 
mentum fraction, or Feynman x, as ; 



Pt 



(right mover). 



(1.5) 



which implies m±/\/2p'^ < x < 1. The rapidity of the pion is then 



1, P 

y = oln 



= 4ln^ 



Ptv 



— Yproj 



-In- 



In 



M 



(1.6) 



(yproj = ln(V2p^/M) ~ ln(y's/M)), and lies in the range < y < 
Yproj + ln(A//m_L). For a left moving pion (pj < 0), we use similar 
definitions where p~^ and p~ are exchanged. This gives a symmetric 
range for y, as in Fig. |21 All the pious are produced in a distribution of 
rapidities within this range. 
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The Color Glass Condensate and High Energy Scattering in QCD 7 

In Fig. 121 dN/dy is the number of produced particles (say, pions) per 
unit rapidity. The leading particles are shown in the solid line and are 
clustered around the projectile and target rapidities. For example, in a 
heavy ion collision, this is where the nucleons would be. In the dashed 
line, the distribution of produced mesons is shown. 



dN 
dy 




Fig. 2. The rapidity distribution of particles produced in a hadronic collision. 



Several theoretical issues arise in multiparticle production. Can we 
compute dN/dy? Or even dN/dy at y = ("central rapidity") ? How 
does the average transverse momentum of produced particles {p±) be- 
have with energy? What is the ratio of produced strange/nonstrange 
mesons, and corresponding ratios of charm, top, bottom etc at y = 
as the center of mass energy approaches infinity? Does multiparticle 
production as s — > cxd at y = become simple, understandable and 
computable? 

Note that y = corresponds to particles with pz = or — 
mj^/\/2, for which x = m±/ (\/2p'^) = m±/^ is small, x ^ 1, in 
the high-energy limit of interest. Thus the multiparticle production at 
central rapidity reflects properties of the small-x degrees of freedom in 
the colliding hadron wavefunctions. 

There is a remarkable feature of rapidity distributions of produced 
hadrons, generally referred to as either Feynman scaling or Limiting 
Fragmentation. If we plot rapidity distributions of produced hadrons at 
different energies, then as function of y — jproj, the rapidity distributions 
are to a good approximation independent of energy. This is illustrated in 
Fig. 13 where the rapidity distribution measured at one energy is shown 
with a solid line and the rapidity distribution at a different, higher, 
energy is shown with a dotted line. In this plot, the rapidity distribution 
at the lower energy has been shifted by an amount so that particles of 
positive rapidity begin their distribution at the same yproj as the high 
energy particles, and correspondingly for the negative rapidity particles. 
This of course leads to a gap in the center for the low energy particles 
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8 The Color Glass Condensate and High Energy Scattering in QCD 

due to this mapping. 

This means that as we go to higher and higher energies, the new 
physics is associated with the additional degrees of freedom at small 
rapidities in the center of mass frame (small-x degrees of freedom) . The 
large x degrees of freedom do not change much. This suggests that there 
may be some sort of renormalization group description in rapidity where 
the degrees of freedom at larger x are held fixed as we go to smaller 
values of x. We shall see that in fact these large x degrees of freedom act 
as sources for the small x degrees of freedom, and the renormalization 
group is generated by integrating out degrees of freedom at relatively 
large x to generate these sources. With this understanding, one may be 



dN 
dy 




yproj ypi'oj 



Fig. 3. Feynman scaling of rapidity distributions. The two different lines correspond to 
rapidity distributions at different energies. 



able to compute the number and distribution of particles produced at 
central rapidities. 

1.5. Deep Inelastic Scattering 

In the previous section, we discussed hadron-hadron collisions in which 
large numbers of particles are produced. Here we will discuss Deep In- 
elastic Scattering (DIS) of a lepton scattering off a hadronic target 
Fewer particles are produced in DIS, so this provides a relatively clean 
environment to study QCD at high energies. In Fig.|l]is shown the car- 
toon of a DIS experiment. 

To describe quark distributions, it is convenient to work in a reference 
frame where the hadron has a large light-cone longitudinal momentum 
P"*" S> M ("infinite momentum frame"). In this frame, one can describe 
the hadron as a collection of constituents ("partons"), which are nearly 
on-shell excitations carrying some fraction x of the total longitudinal 
momentum . (The correct mathematical formulation of this picture 
involves the light cone quantization mentioned previously: the hadron 
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The Color Glass Condensate and High Energy Scattering in QCD 9 



electron 




Fig. 4. Deep inelastic scattering of an electron on a hadron. 

can then be expressed in a "quark-parton" basis.) Thus, the longitudinal 
momentum of a parton is — xP^, with < x < 1. 

For the struck quark in Fig.|l] this x variable ("Feynman x") is equal 
(modulo target mass corrections) to the empirical Bjorken variable xgj, 
which is defined in a frame independent way as x^j = /2P ■ q. In this 
definition, = — g'^q^, with the (space- like) 4-momentum of the 
exchanged photon. The Bjorken variable scales like x^j ~ Q'^ /s, with 
s — the invariant energy squared. Thus, in deep inelastic scattering at 
high energy (large s at fixed Q^) we measure the quark distributions 
dNquark/dy^ at small x (x ^ 1). 

It is again useful to think about these distributions as a function of 
rapidity. We define this for deep inelastic scattering as y = y hadron ~ 
ln(l/x), and the invariant rapidity distribution as dN/dj = xdN/dx. At 
high Q^, the measured quark and gluon distribution functions are simply 
related (at least, in a leading order approximation) to the number of 
partons per unit rapidity in the hadronic wavefunction. 

The typical dN/dy distribution for constituent giuons of a hadron as 
measured in DIS is similar to the rapidity distribution of produced parti- 
cles in hadron-hadron collisions (see Fig. |5J . This suggests a relationship 
between the structure functions in DIS and the rapidity distributions 
for particles produced in hadronic collisions. We expect the gluon dis- 
tribution function to be proportional to the pion rapidity distribution. 
This relation is plausible (since the degrees of freedom of the giuons 
should not be lost) and is incorporated in many models of multi-particle 
production. 

The small x problem is that in experiments at HERA, the rapidity 
distribution functions for quarks and giuons grow rapidly as the rapidity 
difference 

r = ln(l/x) = Yhadron-y (1-7) 
between the quark and the hadron increases In Fig. |K1 the ZEUS 
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The Color Glass Condensate and High Energy Scattering in QCD 

data for the gluon distribution are plotted for = 5 GeV^, 20 GeV^ 
9 7 

and 200 GeV . The gluon distribution is the number of gluons per 
unit rapidity in the hadron wavefunction, xG(x, Q^) = dNgi^Qng/dj . 
Experimentally, it is extracted from the data for the quark structure 
functions, by analyzing the dependence of the latter upon the resolution 
of the probe. 

The growth seen in Fig. |S] appears to be more rapid than r or r^. 
Perturbative considerations of the high energy limit in QCD by Lipatov 
and colleagues lead to an evolution equation commonly called the BFKL 
equ ation ^ which suggests that distributions may grow as an exponential 
in T^^. Alternatively, the double logarithmic DGLAP evolution equa- 
tion ^1 predicts a less rapid growth, like an exponential in Both of 
these evolution equations would predict asymptotically a growth of the 
distributions which would exceed the Froissart unitarity bound discussed 
previously. 




10"^ 10"^ ur 10" 



Fig. 5. The Zeus data for the gluon structure functions. 



How do we understand in QCD the problem of the rapid rise of gluon 
distributions at small x? Consider Fig. El where we view the hadron head 
on. The constituents are the valence quarks, gluons and sea quarks, all 
shown as colored circles. As we add more and more constituents, the 
hadron becomes more and more densely populated. If one attempts to 
resolve these constituents with an elementary probe, as in DIS, then, at 
sufficiently small x (for a given transverse resolution), the density of the 
constituents becomes so large that one cannot neglect their mutual in- 
teractions any longer. One expects such interactions to give "shadowing" 
by which we imply a decrease of the scattering cross-section relative to 
what is expected from incoherent independent scattering. 

More precisely, we shall see later that, as a effect of these interac- 
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Low Energy 



High Energy 



Fia 



Saturation of gluons in a hadron. A view of a hadron head on as x decreases. 



tions, the parton distribution functions at fixed Q saturate, in the s ense 
of showing onl y a slow, log arithniic, increase with l/ :!^^ ' ^'^ 13,14,15 | 16 | 
(See also Refs.E^EHElIMfor recent reviews and more references.) For a 
given Q^, this saturation occurs if x is low enough, lower than some crit- 
ical value Xs{Q^). Converserly, for given x, saturation occurs for trans- 
verse momenta below some critical value Q^{x), defined as 

Ql{x) ^ asNc^^, (1.8) 
vrit^ ay 

where dN/dj is the gl uon distribution s-t y — y hadron — 

ln(l/x). Only 

gluons are included, since, at high energy, the gluon density grows much 
faster than the quark density, and is the driving force towards saturation. 
This explains why in the following we shall focus primarily on the gluons. 
In Eq. 11.811 . ttJ?^ is the hadron area in the impact parameter space (or 
transverse plane). This is well defined provided the wavelength of the 
probe is small compared to R, which we assume throughout. Finally, 
asNc is the color charge squared of a single gluon. Thus, the "saturation 
scale" Ijl.H^ has the meaning of the average color charge squared of the 
gluons per unit transverse area per unit rapidity. 

Since the gluon distribution increases rapidly with the energy, as 
the HERA data suggests, so does the saturation scale. For high enough 
energy, or small enough x. 



Q^(x)>A^CD: 



(1.9) 
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2 The Color Glass Condensate and High Energy Scattering in QCD 



and as{Qs) <C 1. This suggests that weak coupling techniques can be 
used to study the high energy regime in QCD. 

However, weak coupling does not necessarily mean that the physics 
is perturbative. There are many examples of nonperturbative phenom- 
ena at weak coupling. An example is instantons in electroweak theory, 
which lead to the violation of baryon number. Another example is the 
atomic physics of highly charged nuclei. Although the electromagnetic 
coupling constant is very weak, acm ^ 1, the large charge Z of the 
atoms introduces a scale Zacm ~ 1 producing a strong nuclear Coulomb 
field in which the electron propagates. Nevertheless, there is a systematic 
technique which enables one to compute the non-perturbative proper- 
ties of high Z atoms. Yet another example is QCD at high temperature 
T S> f^QCD\ this is a weakly coupled quark-gluon plasma, but exhibits 
nonperturbative phenomena on large distances r 3> l/T due to the col- 
lective behaviour of many quanta' — 

Similarly, the small-x gluons with transverse momenta < 
make a high density system, in which the interaction probability a x n 
(where rr ^ as/ is the typical parto n cross-se ction and n is the gluon 
density, n = ^/nR^) is of order one I^^IMSI (cf. Eq. (jTHJ). That is, 

although the coupling is small, as(Qs) <^ 1, the effects of the interac- 
tions are amplified by the large gluon density and ordinary perturbation 
theory breaks down. A resummation of the high density effects is there- 
fore necessary. Our strategy, to be developed in the following sections, 
■yilli ^jfi, effective theory — the Color Glass Condensate 
| 14 | 25 | 26|23 | 20| — which the small-x gluons are described as the classi- 
cal color fields radiated by "color sources" at higher rapidity. Physically, 
these sources are the "fast" partons, i.e., the hadron constituents with 
larger longitudinal momenta 3> xP"*". The properties of the color 
sources will be obtained via a renormalization group analysis, in which 
the "fast" partons are integrated out in steps of rapidity and i n th e 
background of the classical field generated at the previous stepl^^ESl 



1.6. Nucleus-Nucleus and Proton-Nucleus Collisions 

In Fig. 171 we plot a cartoon of the space-time evolution of a heavy ion 
collision' — '. Imagine we have two Lorentz contracted nuclei approaching 
one another at the speed of light. We choose coordinates such as the 
collision takes place at 2; = t = 0, or a;"*" = x~ = 0. Since the two nuclei 
are well localized in the longitudinal direction, they can be thought of as 
sitting at z ~t {x^ = 0) for the right mover, respectively at z ~ — t (or 
— 0) for the left mover. To analyze this problem for i > 0, namely, 
after the collision takes place, it is convenient to introduce a time variable 
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z 



Fig. 7. A space-time figure for ultrarelativistic heavy ion collisions. 



which is Lorentz invariant under longitudinal boosts'^ t = — (the 
"proper time") and a space-time rapidity variable 



For free streaming particles with velocity Vz (z = Vzt = ^t), the space- 
time rapidity equals the momentum space rapidity (11.411 : rj — y. We 
shall see later that this identification remains approximately true also 
for the off-shell quantum fiuctuations (the partons) with relatively large 
longitudinal momenta. 

At high energies, in the central rapidity region, particle distributions 
vary slowly and it should be a good approximation to take them to be 
rapidity invariant and therefore also independent of rj. Therefore distri- 
butions are the same on the lines of constant proper time r, which are 
shown in Fig.[7| 

An outstanding problem is to formulate the initial conditions for a 
heavy ion collision and to study the subsequent evolution of the produced 
partons. Can one argue from first principles that the partonic matter will 
thermalize into a quark-gluon plasma ? 

There are two separate classes of problems one has to understand 
for the initial conditions. Firstly, the two nuclei which are colliding are 
coherent quantum mechanical wavepackets. Therefore, for some early 
time, the degrees of freedom must be quantum mechanical. This means, 
in particular. 




(1.10) 



AzApz > 1 , 



(1.11) 



^This should not be confused with the rapidity variable introduced previously, in 
eq. il.7i . and which will not appear in this subsection. 
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which particularly constrains the small-x gluons, which are delocalized 
over large longitudinal distances Az ~ 1/pz, and thus overlap with each 
other. Such degrees of freedom cannot be described by semi-classical 
transport theory for particles. (Classical particles are characterized by a 
distribution function f{p,x,t), which is a simultaneous function of mo- 
menta and coordinates.) However, fortuitously, quantum coherent states 
can be described as classical fields because they have large occupation 
numbers ~ 1/qs 3> 1. Heisenberg commutators between particle cre- 
ation and annihilation operators become negligible in this limit: 

[ak, al] = 1 < alak = Nk. (1.12) 

Classical field theory is also the appropriate language to describe another 
important feature of the initial conditions, namely the classical charge 
coherence. At very early times, we have a tremendously large number 
of particles packed into a longitudinal size scale of less than a fermi due 
to the Lorentz contraction of the nuclei. We know that such particles 
cannot interact incoherently. For example, if we measure the field due 
to two opposite charges (a dipole) at a distance scale r large compared 
to their separation, the field falls off as not 1/r. On the other 

hand, in parton cascade models, interactions are taken into account by 
cross-sections which involve matrix elements squared leaving no room for 
classical charge coherence. These models should therefore not be applied 
at very early times. 

As an effective theory at small x, the Color Glass Condensate can be 
applied to study the initial stages of a heavy ion collision' ^ ■ ^" ■ ■-''-' l The 
only scales in the problem are the saturation scale Qs and the transverse 
size of the system R. On a time scale r ~ 1/Qs, the initial energy and 
number distributions of gluons can be computed by solving classical field 
equations. At later times r 3> 1/Qs, the system becomes dilute and the 
classical field approximation breaks down, but one in terms of transport 
equations becomes valid. There is presum ably an overlap region where 
the two descriptions can be both correct At early times, the dis- 
tribution of gluons in momentum space is primarily transverse. As the 
system becomes dilute, the gluons begin to scatter "off the transverse 
plane" . Baier, Mueller, Schiff and Son 34 have argued that 2~* 3 pro- 
cesses are those which lead to the most efficient thermalization. Whether 
on not the initial non-equilibrium gluon distributions thermalize to form 
a quark gluon plasma at the energies of interest is of outstanding phe- 
nomenological interest at RHIC and LHC energies. It is also of interest 
to understand how mu ch of the RH IC data can be understood purely 
from initial state effects E ^ l '^7 | 38 | opposed to the final state effects 
which are important for thermalization. Proton (or deuterium)-nucleus 
collisions should be very helpful in helping to isolate initial state effects 
from final state rescattering effects. In a proton-nucleus collision, one 
does not expect final state rescattering to dominate the measured parti- 
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cle spectrum. Initial state effects, on the other hand, should be especially 
important in the proton fragmentation region and have been computed 
recently within the CGC framework ' * l ^v l ^^ l ^^l The upcoming exper- 
iments at RHIC on deuterium-gold collisions will hopefully clarify the 
role of initial/final state effects in high energy scattering. 

1.7. Universality of High Energy Scattering 

In pion production, it is observed that, with the exception of globally 
conserved quantities like the energy and the total charge, the rapidity 
correlations are of short range. If the theory is local in rapidity, the 
only dimensionful parameter which can determine the physics at a given 
rapidity is Q^{x). For an approximately scale invariant theory such as 
QCD, a typical transverse momentum of a constituent will therefore be 
of order Qg. If S> l/R^, where -R is the radius of the hadron, the 
finite size of the hadron becomes irrelevant. Thus at small enough x, all 
hadrons become the same — specific properties of the hadrons (like their 
size or atomic number A) enter only via the saturation scale Q'^(x,A). 

Hence there should be some equivalence between nuclei and protons: 
when their values are the same, their physics should be the same. 
Eq. 11.811 suggests the following empirical parametrization of the satura- 
tion momentum: 

41/3 

qI{x,A) ~ (1.13) 
x-^ 

wher e the value A ~ 0.2 — 0.3 seems to be preferred by both the data 
l ^ l '^'^ l and the most recent theoretical calculations There should thus 
be the following correspondence: 

• RHIC with nuclei ~ HERA with protons 

• LHC with nuclei ~ HERA with nuclei 

Estimates of the saturation scale for nuclei at RHIC energies give 
Qs^ 1-2 GeV, and at LHC 2-3 GeV. 

This further suggests that for relatively simple processes like deep 
inelastic scattering (at least), the observables should be universal func- 
tions of the ratio between the transferred momentum and the satu- 
ration scale Q^{x,A). This feature, called "geometric scaling", is rather 
well satisfied by the proton structure functions measured at HERA, for 
all values of Bjorken x smaller tha n 0.01 and in a broad region of 
(between 0.045 and 450 GeV^ ) Interestingly, the observed scaling 
extends up to relatively large values of Q^, well above the saturation 
scale Qs(x), which suggests that the phenomenon of gluon saturation at 

< Qs (x) has a rather strong infiuence also on the physics at much 
larger Q^m 

The previous considerations suggest how to reconcile unitarity with 
the growth of the gluon distribution function at small x. The point is 
that the smaller x is the larger Qs (x) is and thus and the typical partons 
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are smaller. Therefore, when decreasing x, although we are increasing 
the number of gluons, we do it by adding in more gluons of smaller and 
smaller size. A probe of transverse size resolution Axj_ ~ 1/Q will not 
see partons smaller than this resolution size. Thus when Q < Qs, newly 
created partons will not contribute to the cross-section. 

2. The Effective Theory for the Color Glass Condensate 
2.1. The hadron wavef unction at small x 

As discussed in Sect. 1.2, a relativistic quantum field theory quantized 
on the lightcone has a simple structure in terms of the bare quanta or 
"partons" of the theory. For instance, in QCD, the eigenstates of the 
light-cone Hamiltonian Pqqjj can be expressed as a linear superposi- 
tion of the eigenstates of the non-interacting part of the Hamiltonian 
Pqqu = ^QCD ~ ^QCD, where Vqcd is the LC potential. The proton 
wavefunction for instance can be written in this parton basis as 

\tp >= ci\qqq > +C2\qqqg > H h Cn\qqqgggg...qqgg > H . (2.1) 

It is a priori not obvious what the advantage of such a decomposition is. 
The magic of the parton model is that, at high energies, it is apparent 
from the structure of Vqcd that the interactions of the partons with 
each other are time dilated. The very complicated picture of the scat- 
tering of the proton off an external potential can be replaced with the 
simple picture of individual partons scattering off the external potential. 
(In the eikonal approximation, each of these partons acquires a simple 
phase in the scattering At lower energies, the Fock states involving 
large numbers of partons are not very important but at higher ener- 
gies they are increasingly important, and involve predominantly gluons. 
Since the total LC momentum of the proton is divided among a large 
number of partons, a typical parton will carry a momentum fc^ << P^ , 
i.e., only a small fraction x = / P^ of the proton total momentum. 
Understanding the physics of QCD at high energies or small x thus re- 
quires that we understand the properties of the n-gluon components of 
the hadron's LC wavefunction. 

A very interesting approach to computing the properties of the n- 
gluon component of the hadronic wavefunction is Mueller's color dipole 

AO 

approach for heavy quarkonia . This is valid in the limit where the 
number of colors Nc is large, so the gluons can be effectively replaced 
by qq pairs ("color dipoles"). It will be shown later that this approach 
and the CGC formalism give identical results for the evolution of distri- 
butions in X, in this large- A'^c limit. 

The CGC approach that we shall follow here is to construct a coarse 
grained effective theory for the small-x component of the hadron LC 
wavefunction. We shall first consider a large nucleus, for which this con- 
struction is most intuitive. Then, we shall argue that, at sufficiently large 
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energy, or small enough x, a similar theory can be constructed for any 
hadron, via weak coupling calculations in QCD. 

2.2. The McLerran-Venugopalan model for a large nucleus 

Consider a nucleus in the infinite momentum frame (IMF) with momen- 
tum P+ ^ oo. We will assume that the nucleus is of nearly infinite 
transverse extent with a uniform nuclear matter distribution. As wo will 
discuss later, the model can be extended to include realistic nuclear 
density profiles. In the IMF, partons which carry large fractions of the 
nuclear momentum ("valence" partons), are Lorentz contracted to a dis- 
tance ~ 2i?/7, with 7 = /mp and mp the mass of the proton. The 
"wee" partons with momentum fractions a; ^ 1 arc dclocalized in the 
x~ direction over much larger distances: a simple kinematic estimate 
suggests that partons with momentum fraction x <^ A~^^^ are dclocal- 
ized over distances larger than the nuclear diameter. These partons then 
"see" the partons at large x as infinitely thin sources of color charge. The 
model assumes a simple kinematic distinction between wee and valence 
partons. The reality is more complex as we will discuss later but certain 
key features of this simple model will survive. 

Another salient feature of wee partons is their very short lifetimes. 
Their lifetime on the liglitconc, Aa;^, is, from the uncertainity prin- 
ciple, conjugate to their light cone energy: Ax~^ ~ l/k~. From the 
light cone dispersion relation, k~ = rn]_/2k'^ = m^/2xP"'". Thus 
Aa;"*" ~ xP"'"/m5^ is proportional to x. These estimates imply that wee 
parton lifetime is much shorter than that of the valence partons: on 
the scale of the lifetime of the former, the latter appear to live forever. 
The valence parton sources are thus static sources of color charge. Since 
their momenta are large, they are unaffected by absorbing or emitting 
soft quanta: they are recoilless sources of color charge. In this "eikonal" 
approximation, the wee parton cloud couples only to the "plus" compo- 
nent of the LC current, which, from the discussion here, can be written 
as (see also Sect. 2.3 below) : 

^ s^'+5{x-)p''{x_l) , (2.2) 

where p''(x^) is the valence quark color charge density in the trans- 
verse plane. The J-funtion in x~ assumes an infinitely thin sheet of color 
charge. The assumption can and must be relax;ed; namely, p'^{x±^) — > 
p°'(x^,x~) as wc will discuss later. Note that is static, i.e., indepen- 
dent of the LC time x^ , for the reasons explained previously. 

We now turn to the color charge density p'^{x±^) and how it is gener- 
ated. We assume that the nucleus is interacting with an external probe 
which can resolve distances of size Ax± in the transverse plane that 
are much smaller than the nucleon size ~ ^QCD- Now, in the longitu- 
dinal direction, the small probe which has x <^ A~^^^ simultaneously 
couples to partons from nucleons all along the nuclear diameter. Since 
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its transverse size is much smaller than the nucleon size, it sees them 
as sources of color charge. If the density n = NcA/ttR^ ~ Aqqj^A^^^ 

{Ra = RqA^^"^ is the radius of the nucleus) of the valence quarks in the 
transverse plane is large, n S> ^qcD' ^^'^ if <^ AS± <^ ^/^QCD, 
with AS± ~ {Ax±)'^, then the number AA*' = nAS± of color charges 
within the tube of transverse area ASj_ crossing the nucleus is large: 
AN S> 1. These sources belong typically to different nucleons, so they 
are uncorrelated with each other, because of confinement. They are thus 
random sources of color charge and the total color charge Q° in the tube 
is the incoherent sum of the color charges of the individual partons. 
Thus, 

(Q") = 0, (Q"Q") = g^CfAN = AS^_ (2.3) 

where we have used the fact that the color charge squared of a single 
quark is g^t°'t'^ = g^Cj. One can treat this charge as classical since, 
when AA^ is large enough, we can ignore commutators of charges: 
I [S", Q'] 1 = 1 ir'"Q" I « Q". 

Let us introduce the color charge density p°'{x~,x±^) in such a way 
that: p°'{x±) = J dx~ p°-{x" ,x^). Then, 

Q'' = cfx^p'^{x^) = / (fx^ / dx~ p'^{x~ ,x^), (2.4) 

Jas± jas± J 

and eqs. llO imply (recaU that Cj = (A^c - l)/2Afc) : 

^ A 

{paix±)pi,{y±))A = SabS''^^ {xj_ -yi_)p\, p\- 



2%R\ 



{pa{x ,x±^)ph{y ,y±))A = Sab3^ {3:±-y±)5{x -y )A^(x ), 

j dx^ \a{x^) = p\- (2.5) 

Here, /i^ ~ A^^'^ is the average color charge squared of the valence 
quarks per unit transverse area and per color, and \a{x~) is the cor- 
responding density per unit volume. The latter has some dependence 
upon x~ , whose precise form is, however, not important since the final 
formulae will involve only the integrated density p^A ■ There is no explicit 
dependence upon x±^ in p\ or Ayi(a::~) since we assume transverse ho- 
mogeneity within the nuclear disk of radius Ra- Finally, the correlations 
are local in x~ since, as argued before, color sources at difi'erent val- 
ues of x~ belong to different nucleons, so they are uncorrelated. All the 
higher-point, connected, correlation functions of pa(x) are assumed to 
vanish. The no n-trivia l correlators 12.511 are generated by the following 
weight function (with x — {x~ ,x±)) : 
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which is a Gaussian in pa, with a local kernel. This is gauge-invariant 
(since local) , so the variable pa in this expression can be the color source 
in any gauge. The choice of a gauge will however soon become an issue 
when we shall study the dynamics of the gluons radiated by this random 
distribution of color charges. 

The local Gaussian form of the weight function in Eq. 1)2. 6|l is valid, 
by construction, for a large nucleus, and within some restricted kine- 
matical range that we spell here again, for more clarity. As already dis- 
cussed, this is correct for a transverse resolution = 1/AS± within 
the range ^qcd ^ ^QCD-^ ■ But, clearly, the assumption that 

the valence quarks are uncorrelated must fail for transverse separations 
of order Rq ~ l/^QCD or larger, since the Nc = 3 valence quarks within 
the same nucleon are confined in a color singlet state. Thus, the total 
color charge, together with its higher multipolar moments, must vanish 
when measured over distances of the order of the nucleon size iig, or 
larger. As emphasized by Lam and Mahlon EHl (gee also Ref. '^^^ for an 
earlier discussion), the requirement of color neutrality can be included 
in the Gaussian weight function by replacing the (5-function in eq. 12.51 
with (p"" {x~ , X j_) p^ (0)) = \{x~ , Xj_)S°'^, where A(a:~, a;_L) is such that its 
Fourier transform X{x~,kj_) vanishes rapidly at momenta kj_ ^ ^qcd- 

Consider also the validity range of Eq. 12.611 in longitudinal momenta. 
As explicit in the previous analysis, the color fields that we are comput- 
ing have small values of a; <C A~^^^, that is, they carry longitudinal 
momenta k'^ much lower than those of their sources, the valence quarks. 
Thus, the weight function 12.611 can be seen as part of an effective theory 
for gluon correlations at momenta k'^ smaller than some upper cutoff 
A"*", of the order of the typical longitudinal momentum of the valence 
quarks. (This theory will be completely specified in Sects. 2.3-2.4.) Note 
however that, for the classical approximations underlying Eq. 12.611 to 
be valid, the value k~^ of interest should be not much smaller than A"*". 
Indeed, as we shall see, new color sources with momenta < A^ are 
produced by radiation from the original sources at p"*" > A'^. If the gap 
between fc"*" and A"*" is relatively large — the precise condition is that 
asNcln{A+/k+) > 1 — these new sources, which are mostly gluons, will 
completely dominate the physics at the scale k~^ of interest. As we shall 
in Sect. 3, these new sources can be explicitly constructed by integrating 
out layers of quantum fluctuations in a renormalization group analysis, 
but the ensuing weight function is generally not a Gaussian. Still, the 
Gaussian 12.611 may be a good initial condition for this quantum evolu- 
tion. Moreover, rather remarkably, it turns out that this is also a good 
approximation in the quantum theory, but with a non-trivial transverse 
momentum dependence for the 2-point function A that will be specified 
in Sect. 3. 
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2.3. The Color Glass 

Once the weight function for the classical color charge configurations 
associated with the large-x partons is known, it is possible to write down 
an effective theory for the small-x gluons. The gener atin g functional for 
the correlation functions of the small-x gluons reads 



f^^ VA5{A+) e'SlA,p] 



z[j] = I vp m/a+ [p] { r , } , (2.7) 



where the external current ja is a formal device to generate Green's 
functions via differentiations, and S[A, p] is the action that describes 
the dynamics of the wee gluons in the presence of the classical color 
charge p (see below). The path integral over the gluon fields is written 
in the light cone (LC) gauge Aa = 0, since this is the gauge which allows 
for the most direct partonic interpretation. Correspondingly, pa is the 
color charge density in the LC gauge. 

Note the dependence upon the intermediate longitudinal momentum 
cutoff A"^ in the integrals in Eq. (12.711 . As explained at the end of the 
previous subsection, A^ is the scale which separates the 'fast' partons 
(p^ > A^"), which have been 'integrated out' and replaced by the clas- 
sical color charge pa, from the 'wee' gluons (fc^ < A"*"), for which the 
effective theory is meant, and which are still explicit in the path inte- 
gral, as the gauge fields Aa- Since obtained after integrating out the 
modes with > A^ (see Sect.|3J, the weight function W/\^+ [p] depends 
upon the separation scale A^. Clearly, A^ must be chosen such that 
A^/P+ > X, with X the longitudinal fraction of interest. 

For instance, for a large nucleus at not so high energies (e.g., a gold 
nucleus at RHIC) we have x ~ 10~^ — 10"'^, so we can neglect the 
quantum evolution in a first approximation, and identify p with the 
color charge of the valence quarks. Then, we can use the MV model 
described previously, in which case Wjy+ [p] = Wa [p] is the Gaussian 
weight function given by Eqs. 12.61 and 12.51 . 

But for smaller values of x (say x < 10~^), as relevant for DIS at 
HERA, and also for nucleus-nucleus scattering at LHG, quantum effects 
are essential, so p must include the color sources generated by quantum 
evolution down to A^. As we shall see in Sect.|3 these sources are pre- 
dominantly gluons. Thus, in what follows, we shall restrict the quantum 
evolution to the gluonic sector; that is, the only fermions to be included 
among the color sources are the valence quarks in the initial condition. 

If the scale — xP^ of interest is of order A"*", or slightly below it, 
the correlation functions at that scale can be computed in the classical 
approximation, i.e., by evaluating the path integral over ^4'' in Eq. 12.71 
in the saddle point approximation. There are two reasons for that: a) 
The quantum corrections due to gluons in the intermediate range fe^ < 

< A"*" are, at most, of order as ln(A^/fc^), and thus are truly higher 
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order effects as long as ln(A~'"/fc~'') ^ 1. b) The small-x gluons liave large 
occupation numbers ^ 1 (as we shall see, at saturation the occupation 
numbers are parametrically of order 1/qs), so their mutual interactions 
can be treated in the classical approximation, i.e., by solving the classical 
field equations SS/5A^ = 0, or, more explicitly: 



where Du = du - igA^T" with iT'')bc = -ifabc- 

Note however that, precisely because the occupation numbers are so 
large, the corresponding classical fields are strong (A' ~ 1/g at satura- 
tion), and thus the classical non-linearities must be treated exactly. In 
particular, we need the exact solution to the classical equations of mo- 
tion (EOM) 12. 81 . that we shall construct in the next subsection. Once 
this solution is known as an explicit functional of p, the correlation func- 
tions of interest are obtained by averaging over p, with weight function 
For instance, the 2-point function is computed as; 



where Aa = Aa [p] is the solution to Eq. 12.81 , and is independent of the 
LC time (because so is the source pa(x)). This means that only equal- 
time correlators can be computed in this way; but these are precisely the 
correlators of interest for small-x scattering. 

The remaining question is: What is the weight function Wf^+ [p] for 
<C ? As we shall see in Sect.|3| this can be answered via a quan- 
tum calculation, the result of which will allow us to compute in 
terms of the initial condition + , by integrating out the gluons with 

momenta in the range A^ < < A^ . For the purpose of this calcu- 
lation, we need to also specify the action S[A, p] in Eq. 112.811 . The simple 
guess S[A, p] = S'YMi^] + / paAa (which would generate the classical 
EOM 12.81 ') cannot be right since the second term, involving p, is not 
gauge-invariant. This reflects the fact that, written as in Eq. lHISJ, the 
classical EOM is not manifestly gauge-covariant either. In fact, Eq. 12.81 
is correct as written only for field configurations having Aa ~ 0. This 
is not a limitation for classical calculations, since it is always possible 
to construct a classical solution having this property. But Eq. 12.81 is 
not sufficient to determine S[A, p], which is explicitly needed for the 
quantum calculation. 

To find the general equation which replaces Eq. 12.81 in some ar- 
bitrary gauge (where Aa 7^ 0), notice that, in general, the current Ja 
in the r.h.s. of the Yang-Mills equations must satisfy the covariant con- 
servation law D^J^ — (since we also have D^DvF^^ — 0). For the 
eikonal current 12.21 . this implies D~ J'^ = {-^pp —igA^)J'^ — 0, which 
reduces indeed to d p = when A~^ — 0. But in general, this is satisfied 



Pa{x) , 



(2.8) 
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by J^{x^ ,x) We have introduced here the 

temporal Wilson line: 

W[A~]{x+,x) = T expi^ig J dz+A" (2+, f ) | , (2.10) 

with T denoting the ordering of the color matrices in the exponential 
w.r.t. their a; "^ arg uments. The action generating the EOM with current 
ijmTjl reads 1^3211 

S[A, p]^- j d^x\ F^,Fr + J^J d'f Tr {p(f) W[A-]{x)], (2. 

where M^[yl~](a;) is gi ven by eq. I2.1UII with x'^ 00. This action is 
gauge-invariant indeed (Another gauge invariant generalization for 
pA~ , namely Tr|p(f) ln(W[yl~])(f )| , has been checked in Ref.l^to 

give equivalent results.) 

The mathematical struct ure of the average over p in Eqs. 12.711 and 
12.91 is that of a Color Glas^^^^. Note indeed the special structure of 
the 2-point function that follows from Eq. 12.71 : 



11) 



lTA>^(^A-(^\ f-n W r J ./^^P-4 A^(x)A'-(y)e'^[^-Pl 
{TA^ix)A iy))^.=JVpW^Ap]^ jA+p^^.,5[A,p] 

This is not the same as : 

JVp Wa[p] J^VAA''{x)A''(y) e'^^^^^P^ 
JVpWAlp]J^'DAe^S[A,p] 



(2.12) 



(2.13) 



The physical reason for this is the fundamental separation in time be- 
tween the rapidly varying wee gluons and the comparatively 'frozen' 
large-x partons. One thus solves the dynamics of the wee gluons at a 
fixed distribution of color charges, and only then averages over the latter. 
There is no feedback from the evolution of the sources on the wee gluon 
fields. And there is no interference between successive configurations of 
the color sources. These features, together with the large fluctuations in 
the color charge density, are the ultimate reasons for treating the large-x 
partons as forming a classical random distribution. 

The prototype of a glass is the "spin glass" -a collection of magnetic 
impurities randomly distributed on a non-magnetic lattice. The dynam- 
ical degrees of freedom, which are rapidly varying, are the magnetic 
moments of the impurities (the "spins"), while the slowly varying disor- 
der refers to the positions of these spins in the host lattice. To study the 
thermodynamics of such a system, one first computes the free-energy (= 
the logarithm of the partition function) of the spin system for a fixed 
disorder (namely, for a given spatial configuration of the impurities), 
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and subsequently makes an average over all such configurations, with 
some weight function. The final average over the configurations is not 
a thermal one: what is averaged is the free-energy computed separately 
for each configuration. 

Similarly, the connected correlation functions of the small-x gluons in 
the present effective theory are obtained from the following generating 
functional: 

where the logarithm is taken inside the integral over p. That is, the free- 
energy reaches its extremum as a function of the external source j for 
a fixed distribution of the color sources. The measured free-energy (or 
correlation function) is finally obtained by also averaging over p. 

Note that the presence of a non-trivial color charge background 
breaks gauge symmetry explicitly. But this symmetry is restored in the 
process of averaging over p provided the weight function is gauge invari- 
ant, which we shall assume in what follows. 



2.4. The classical color field 

In this subsection, we shall construct the solution to the classical EOM 
12.81 1. We note first that, for a large class of gauges, one can always find 

2(1 

a solution with the following properties : 

F'a^ =0, A- ^ 0, A+, A'a : static , (2.15) 

where "static" means independent of . This follows from the specific 
structure of the color source which has just a component, and is 
static. Since F"^^ = 0, the transverse fields A^ form a two-dimensional 
pure gauge; that is, there exists a gauge rotation U{x^,x^_) £ SU(A'^) 
such that: 

A\x~,xj^) = ^U{x~,xj_)d'^u''{x^,xj^). (2.16) 

(in matrix notations appropriate for the adjoint representation: A^ = 
A^T°, etc). Thus, the requirements 12.151 leave just two independent 
field degrees of freedom, A~^{x) and U(x), which are further reduced to 
one (either A'^ or U) by imposing a gauge-fixing condition. 

We consider first the covariant gauge d^A^ = 0. By eqs. 12.1511 and 
12.16II . this implies diA^ = 0, or 17 = 0. Thus, in this gauge, Aa(x) = 
5^^ aa{x~ , xx), with aa{x) linearly related to the color source pa in the 
COV-gauge : 



— V\aa{x) = pa{x) 



(2.17) 
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The only non-trivial field strength is the electric field f^a^ = —d^aa- 
Eq. 12.1711 has the solution : 

aa{x^,Xj_) = J <fy^ {xj_\ — \yi_} Pa{x^ ,yi_) 

'^'^■'1^7 ^-^T^Paix-,y^), (2.18) 



47r (x^ - ?/j_)2^2 

where the infrared cutofi^ fi is necessary to invert the Laplacian operator 
in two dimensions, but it will eventually disappear from (or get replaced 
by the confinement scale Aqqjj in) our subsequent formulae. 

We shall need later the classical solution in the LC-gauge A'^ = 0. 
This is of the form Aa ~ 5^''Aa with Aa{x~,x±^) a "pure gauge", cf. 
eq. 12.161 . The gauge rotation U{x) can be most simply obtained by a 
gauge rotation of the solution in the COV-gauge: 

Af" = U(^A'' + ^d^'^U^ (2.19) 
where the gauge rotation U{x) is chosen such that A'^ = 0, i.e., 

U\x~',xj_) ^Pexplig [ dz^ aa{z~ , x±)T'' \ . (2.20) 



The lower limit Xq —> — oo in the integral over x in eq. 12.2UII has been 
chosen such as to impose the "retarded" boundary condition: 

Aa{x) ^0 as x^ ^ -oo, (2.21) 

which will be useful in what follows. (Note that the "retardation" prop- 
erty refers here to a;~, and not to time.) 

Together, eqs. 12.16L 12.1811 and 12.2011 provide an explicit expression 
for the LC-gauge solution A^ in terms of the color source p in the COV- 
gauge. This is sufficient for the purpose of computing observables since 
the average in Eq. (12.911 can be re-expressed as a functional integral over 
the covariant gauge color source /5 by a change of variables: 

{A\x+,x)A\x+,y)---)^^ ^ JvpW^+[p]Al[p]Ailp]--- . (2.22) 

Up to now, the longitudinal structure of the source has been arbi- 
trary: the solutions written above hold for any function p°'{x~). For 
what follows, however, it is useful to recall, from Sect. 12.21 that p 
has is localized near x~ =0. More precisely, the quantum analysis in 
Sect. 121 will demonstrate that the classical source at the longitudinal 
scale A"*" has support at x~ within the range < x^ < 1/A^. From 
Eq. 12.181 . it is clear that this is also the longitudinal support of the 
"Coulomb" -field a. Thus, integrals over x~ as that in eq. 12.201 receive 
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a' 




Fig. 8. The longitudinal structure of the color source p and of the classical field solution 
A' for the effective theory at the scale A;+ . As functions of x~ , a and J^^^ are as localized 
as p. 



their whole contribution from x~ in this limited range. Any probe with 
momenta q'^ <^ A"*", and therefore a much lower longitudinal resolution, 
will not be able to discriminate the internal structure of the source. 
Rather, it will see a source/field structure which is singular at — : 
pa{x~ ,xj_) ~ S{x~)pa{x±) (see Fig.|SJ. In particular: 



{x-)'-v{dW^)(x^), 



(2.23) 



A {x~ ,x±) 

where V and are obtained by letting x^ ^ oo in Eq. Il2.20t : 

V\x^) = Pexpjig J dz~ a{z~ , x . (2.24) 



2.5. The gluon distribution 

We denote by G(x, Q^)dx the number of gluons in the hadron wavefunc- 
tion having longitudinal momenta between xP"*" and (x + dx)P'^ , and 
a transverse size Axj_ ~ l/Q- In other terms, the gluon distribution 
xG(x, Q^) is the number of gluons with transverse momenta k± Q per 



unit rapidity (see Refs 



xG(x,(3^ 



[19 20 



for more details) 



dN 



k+=xP+ 



dk+d'^kj^ 

= y d^fee(Q^ - fci)x(5(x- fc+/P+) 
where B(a;) is the step function, k = (fc^,k^) and 



dN 
Wk 



dN 

IFk 



{a',^{x+,k) al{x+,k)) = ^{AUx+,k)Alix+,-k)) 



(2.25) 



(2.26) 
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6 The Color Glass Condensate and High Energy Scattering in QCD 

is the Fock space gluon density, namely, the number of gluons per unit 
of volume in momentum space. 

The difficulty is, however, that this number depends upon the gauge, 
so in general it is not a physical observable. Still, as we will shortly argue, 
this quantity can be given a gauge-invariant meaning when computed 
in the light-cone gauge Aa = 0. Using the fact that, in this gauge, 
Fa^{k) = ik'^ Aa{k), one obtains (with k~^ = xP+ from now on): 

xG(x,g2)^i I ^Q^Q^_kl){Ft{k)Ft{-k)), (2.27) 

which so far does not look gauge invariant. A manifestly gauge invari- 
ant op erator can be constructed by appropriately inserting Wilson lines 
II I In LC gauge, this gauge invariant expression reduces to Eq. (12.271 1 
once the residual gauge freedom of the transverse components of the 
gauge field is fixed by imposing the "retarded" boundary condition 12.211 1 
This particular gauge fixing in the classical field problem has impor- 
tant consequences for the quantum calculation in Sect.|3 in that it fixes 
the ie pres cript ion to be used for the 'axial pole' in the LC-gauge gluon 
propagator i-^. 

We shall need later also the gluon density in the transverse phase- 
space (also referred to as the "unintegrated gluon distribution", or the 
"gluon occupation number"). This is defined as: 

^ (k,)- 1 _ 1 {F^(fc)F^+(^fc)) 

where r = ln(l/x) = ln(P^/A:^). Up to the factor Atv'^ , this is the 
number of gluons of each color per unit rapidity per unit of transverse 
phase-space. (As before, we assume a homogeneous distribution in the 
transverse plane, for simplicity.) 

For illustration, let us compute the gluon distribution of a nucleus 
in the MV model. We start with the low density regime, valid when the 
atomic number A is not too high, so the corresponding classical field 
is weak and can be computed in the linear approximation. By expand- 
ing the general solution 12.161 to linear order in p, or, equivalently, by 
directly solving the linearized version of Eq. 12. St . one easily obtains: 



X(fc) ^ -T+^^^^^^i^, T+\k) C.i^pa{k), (2.29) 



J. 

which together with Eq. 12.51 implies: 

{K^{k)K^{-k))A ^ {pa{k)pa{-k))A = ^R\{N^c " ' (2-30) 
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By inserting this approximation in Eqs 12.2811 and 12.2711 . one obtains 
the following estimates for the gluon density and distribution function: 

^^(fc^)~^, (2.31) 
r ( ^2, {Nj-mA 2 ["^^ dkl asANcCf , 

xG^(x,o ) ^ — ^ = , 

'^QCD J- QCD 

(with as = g'^/in). The integral over in the second line has a log- 
arithmic infrared divergence which has been cut off at the scale ^qcd 
since we know that, because of confinement, there is color neutrality on 
the nucleon size _Ro ~ 1/Aqcd ^'^^ We will argue later that, after taking 
into account quantum evolution, the actual scale for the screening of the 
infrared physics is not Aqcd but the saturation scale Qs. 

Eqs. i^TT^i are in fact the expected results, which could have been 
obtained also by a direct analysis of the gluon radiation by a single quark, 
together with the assumption that gluons radiated by different quarks 
do not interact with each other, so that the total gluon distribution is 
simply the sum of independent contributions from the A x Nc valence 
quarks. This is the Weizsacker-Williams approximation for radiation off 
independent quarks. 

2.6. Gluon saturation in a large nucleus 

According to eq. 12.31II . the gluon density in the transverse phase-space 
is proportional to A^^^, and becomes arbitrarily large when A increases. 
This is however an artifact of our previous approximations which have 
neglected the interactions among the radiated gluons, i.e., the non- linear 
effects in the classical field equations. To see this, one needs to recompute 
the gluon distribution by using the exact, non-linear solution for the 
classical field, as obtained in Sect. 12.41 By using !Fa^{x) = ?7^j,(— 9*a''), 
one can express the relevant LC-gauge field-field correlator in terms of 
the color field in the COV-gauge: 

{Ti\S)Tt\y))A = ( [Uld'^'')_^ [Uld'cf)})^. (2.32) 

One can show that the RHS of this expression can be written as 

((t/^V'"') . [uL^''-")}l - {d'a\£)d^a^^{y)) (c/,\(£)f/ea(y)} 

= S{x- - y-){TT uUmm ( - ^hAix-,x^ - y^)) , (2.33) 

where we have used Uac ~ Uca in the adjoint representation. Here we 
have made use of the following correlation function 

{aaix)ai,{y}) A = S^bSix^ - y^)"fAix^ ,xi_ - y±), 

jA{x-,k^) = -^Xa{x-), (2.34) 
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28 The Color Glass Condensate and High Energy Scattering in QCD 

which follows from p°'{x~ , fc i ) = fc '^ a°'(x~ , kj_) together with Eq. 12.51 . 

Eq. 12.3311 can be proven 1^ using rotational symmetry, the path 
ordering of the Wilson lines in x~ , and the fact that the 2-point function 
of the color fields, Eq. 12.341 . is local in x~ . The trace in Eq. 12.331 . 

SA{x-,x^-y^) = ^—{TrUHx~,x±)U{x-,y^))A , (2.35) 



2.36) 



can be explicitely computed as 

SA{x~,rj^) =exp(-g^Nc ( dz" [7^(2" , 0^) ~ ^a{z^ ,r 

where (cf. Eq. (lOH l 

lA{x-,0±)--,A{x-,r^) = \a[x-) I |^^[l-e*^--^-]. (2.37) 

The above integral over kj^ is dominated by soft momenta, and has 
a logarithmic infrared divergence which, in this classical context, can 
be screened only by confinement at the scale ^qcd- To leading- log 
accuracy, i.e., by keeping only terms enhanced by the large logarithm 
\n{l / r\h^(j j-i) , the precise value of the infrared cutoff is not important, 
and we can also expand the integrand as: 



/ 



1/rl 



^2 



r 



± 



In- 



(27r)2 fc4 J (27r)2 fc4 2 Wtv r^^A^^^ ' 

This gives, with Ha{x~) = JH^ dz^\A{z^), 

SA{x~,r±) ~ expi - ^^^ri/i^(a::~) In -j- ^ > , (2.38) 

which together with Eq. H2.33|l can be used to finally evaluate the glu on 
density in Eq. 12.281 . After simple manipulations, one obtains 1^^^" 



l-exp| - I riQ^ In p^^^ 1 

n \ / j2 —ik±-r± I. ' ±'^qcd J /r, on\ 

VA{k±) = d r^c ^ ^ — ^ ^ , 2.39) 



where 



asNcA = asNc j dx- \a{x-) ~ A^/^. (2.40) 



To study the fej_ -dependence of Eq. 12.391 . one must still perform the 
Fourier transform, but the result can be easily anticipated: 
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Let us first introduce tlie saturation momentum Qs{A) whidi, as 
we shall see, is tiie scale separating between linear and non-linear be- 
haviours. This is defined by the condition that, for r± = 2/Qs{A), the 
exponent in Eq. 12.3911 becomes of order one, which gives: 

q2(A) ~ a,Ncfi\ In^^ ~ A^/HnA. (2.41) 
^QCD 

Note that this is larger than Q^, eq. 12.401 . since we work in the hy- 
pothesis that Qyi ^ ^QCD- Then we distinguish between two regimes: 
i) At high momenta kj^^Qs{A), the integral is dominated by small 
rj_ ^ l/Qs{A), and can be evaluated by expanding out the exponential. 
To lowest non-trivial order (which corresponds to the linear approxima- 
tion), one recovers the bremsstrahlung spectrum of eq. 12.311 1: 

^Aik.) - ^ f - f for » Q,. (2.42) 

ii) At small momenta, ^ Qs{A), the dominant contribution comes 
from large distances 3> l/Qs(A), where one can simply neglect the 
exponential in the numerator and recognize 1/r^ as the Fourier trans- 
form'' of In kj_ : 

^A(fc±) « ^In^fl^ for fc^«g^. (2.43) 

asNc k^ 

Unlike the linear distribution 12.4211 . which grows like A^/^, and IS 
strongly infrared dominated (as it goes like 1/A;^), the distribution in 
Eq. 12.4311 , which takes into account the non- linear effects in the classical 
Yang-Mills equations, rises only logarithmically as a function of both A 
and 1/fc^. This is saturation. At saturation, the gluon occupation factor 
is parametrically of order 1/as, which corresponds to a Base conden- 
sate, and is the maximum density allowed by the repulsive interactions 
between the strong color fields A^ = \/7A\A^ ~ 1/g. When increasing 
the atomic number A, the new gluons are produced preponderently at 
large transverse momenta Qs{A). where this repulsion is less impor- 
tant. This is illustrated in Fig. El 

To clarify the physical interpretation of the saturation scale, note 
that, at short-distances rj_<^l/Qs{A), 



''The saturation scale provides the ultraviolet cutofi for the logarithm in eq. 12.431 since 
the short distances r± ^ l/Qs{A) are cut off by the exponential in eq. 
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Fig. 9. The gluon phase-space density ifiA{f^±) of a large nucleus (as described by the 
MV model) plotted as a function of kj_ for two values of A. Notice the change from a 
l/fc^^ behaviour at large momenta kj_ > Qs to a logarithmic behaviour at small momenta 
k± < Qs- 

is the number of gluons (of each color) having tranverse size r_L per 
unit of transverse area (cf. eq. H2.81|l ). Since each such a gluon car- 
ries a color charge squared {gT°'){gT°') = g^Nc, we deduce that 
asNc^\ In a A is the average color charge squared of the gluons 

^±^QCO ___ 

having tranverse size rj_ per unit area and per color. Then, eq. 12.41|l 
is the condition that the total color charge squared within the area oc- 
cupied by each gluon is of order one. T his is the original criterion of 
saturation by Gribov, Levin and Ryskin^^ for which the MV model 
offers an explicit realization. 

Let us finally compute the (integrated) gluon distribution in the sat- 
uration regime, i.e., for <^ Q'^{A), and compare with the correspond- 
ing result in the linear regime at high Q^, Eq. 12.311 . By using 12.431 in 
Eqs. 12.271 - 12.281 . one immediately finds: 

xGa{x,Q ) ~ —r- RaQ In ^ . (2.45) 

47rA'c as Q'^ 

This shows strong nuclear shadowing : it scales like A^^"^ In A, unlike the 
linear result H2.31^ . which scales like A. 

2.7. Dipole-hadron scattering at high energy 

Although the previous definition of the gluon distribution in terms of 
the Fock space gluon density is useful for a conceptual discussion of 
saturation, it is on the other hand less clear whether it corresponds to 
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something that could be directly measured in experiments. Recall that, 
in the lowest-order analysis of deep inelastic scattering where one ne- 
glects non-linear effects in the hadron wavefunction, the gluon distribu- 
tion is related to the scaling violation in the hadron structure function 
F2 : dF2/dlnQ^ oc asxG{x,Q^). It is therefore interesting to compute 
this quantity also in the presence of non-linear effects, and identify some 
measurable consequences of saturation. This is what we shall do start- 
ing with this section, first in the framework of the MV model, then by 
including the effects of the quantum evolution in going towards smaller 
values of x (in Sects. 3 and 4). As a general conclusion, we shall find 
that saturation effects in the hadron wavefunction correspond to uni- 
tarity effects in the high-energy virtual-photon-hadron collision, whose 
total cross-section is related to F2 via; 

^r'H{^,Q^) = ^^^^F2{^,Q^). (2.46) 

As we have seen in Sect. 11.51 when viewed in the infinite momentum 
frame (IMF) of the hadron, DIS appears as the scattering of the virtual 
photon off a quark (with longitudinal fraction equal to the Bjorken x of 
the collision) in the hadron wavefunction. At very small x, this quark is 
typically not a valence quark, but rather a see quark which is emitted, 
most probably, off the small-x gluons. It is then convenient to disentangle 
this final quark emission from the quantum evolution which involves 
mostly gluons. This can be done by performing a Lorentz boost in such 
a way to pull the 'y*qq vertex out of the hadron. That is, inste ad o f the 
hadron IMF, it is preferable to use the so-called ...o;e /.a..eIIi (and 
references therein) in which most of the energy is still carried by the 
hadron (so that the high density effects are again associated with the 
hadron wavefunction), but the virtual photon moves in the negative z (or 
positive x~) direction with enough energy to dissociate before scattering 
into a quark-antiquark pair in a color singlet state (a color dipole), which 
then scatters off the hadron. This sequential picture of DIS is appropriate 
at high energy, since the lifetime of the qq pair is much larger than the 
interaction time between this pair and the hadron. 

More precisely, if r = ln(l/a;) = — y^* is the (boost-invariant) 
rapidity gap, with r 2> 1 at small x, then the dipole frame corresponds 
to choosing y^* ^ yjj and such that asy-y* <C 1, which ensures that one 
can neglect additional gluon radiation from the quark and the antiquark 
in the dipole. 

This physi cal pictur e translates into the following factorization for- 
mula for cr..y.jyE2ISHISl^ that we shaU derive in Sect. 4.1: 

<^j'h{t-,Q'^) = dz j d^r^\^{z,r^;Q'^)\'^ ad^po\cij,r±_). (2.47) 

Here, '^{z,r^;Q^) is the light-cone wavefunction for the photon split- 
ting into a qq pair with transverse size r_L and a fraction z of the pho- 
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ton's longitudinal momentum carried by the Furthermore, 
'^dipolo('''i is the dipole-hadron total cross section for a dipole of 
transverse size = — (with the quark located at xj^ and the 
antiquark at yj_), and is obtained by integrating the (imaginary part) of 
the scattering amplitude J^T{x±,yj_) = A/V (f'_L, 6_l) over all the impact 
parameters fej^ = {x^ + y±)/2 : 

'^dipole (^^'^i) = 2 J dh^Ur{r^,h^). (2.48) 

At high energ y, the dipo le-hadron scattering can be treated in the eikonal 
approximation ^ ' ' < I v v l ^V^ I This amounts to neglecting the recoil of the 
quark (or the antiquark) during its scattering off the color field in the 
target: the whole effect of the scattering consists in a color precession. 
Then, the scattering amplitude reads Nt {x^,y^^ = 1 — (a;^ , ) , with 
the following S-matrix element: 

Sr{x^,y^) = l-^(tr{v\x^)V{y^)))^, (2.49) 

where and V are the Wilson lines describing the eikonal interaction 
between the quark (or the antiquark) and the color field at rapidity r 
due to color sources within the hadron: 

V\x^;) = Pexp^iff j dx^ At{x^ ,xi;)t°^ . (2.50) 

That is, this is the same as Eq. 12.2411 . but rewritten in the fundamental 
representation. The average in Eq. 12.4911 is over all the configurations of 
the color fields in the hadron. In the CGC formalism, where = aa, 
cf. Eq. 12.1711 . this average is computed as in Eq. 12.221 . 

In what follows we shall focus on the computation of the S-matrix 
element 12.4911 . which encodes all the information about the hadronic 
scattering, and thus about the non-linear and quantum effects in the 
hadron wavefunction. [Once this is known, F2 can be immediately ob- 
tained by using Eqs. 12.4711 and 12.461 .1 

In the MV model, to which we shall restrict in the remaining part 
of this subsection, this S-matrix element is already known, as obvious 
when comparing Eqs. 12.4911 and 12.351 . By translating Eq. 12.381 to the 
fundamental representation {Nc Cp = {N^ — l)/2Nc) and letting 
x~ ~* <x there, one obtains: 

SA{ri_) ^ eJ-dOlin-^-^], (2.51) 

where Q\ = asCpfi\ differs only via a color factor from Eq. 12.401 . 
As in the previous discussion of the gluon distribution, we distinguish 
between a small-r_L and a large-r_L regime, with the separation between 
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the two regimes given by the saturation scale Q^{A), defined by analogy 
with Eq. 

i) A small dipole, with <C l/Qs{A), is only weakly interacting 
with the hadron: 

^A{r±) = 1 - SA{ri_) « i rlQl In ^-i « 1, (2.52) 

a phenomenon usually referred to as "color transparency" 

ii) A relatively large dipole, with 2> (but rj_ <^ ^/^qcd 
though, for the present perturbative treatment to apply) is strongly ab- 
sorbed : Sa ^ 1 or A/'yi(rj_) « 1 which corresponds to the unitarity (or 
"black disk") limit. This is in sharp contrast with the single-scattering 
approximation in Eq. H2.52|l . which would predict a scattering ampli- 
tude growing indefinitely with the area spanned by the dipole, and 
also with the atomic number A of the target. Eq. 12.511 shows that, 
when rj_ and/or A are large enough, the effects of multiple scattering 
become non-negligible, and eventually ensure the unitarization of the 
scattering amplitude at fixed impact parameter'^. Note that, according 
to Eq. H2.51^ . the multiple scatterings are higher-twist effects, i.e., their 
contributions are of higher order in . 

We see that, remarkably, it is the same scale — the saturation mo- 
mentum — which sets the critical transverse size for both gluon satura- 
tion and the unitarization of the dipole-hadron scattering. This conclu- 
sion, that we have found here in the framework of the MV model, will 
be seen in Sect. 3 to remain valid after including the quantum evolution. 
Physically, this can be understood as follows: Both saturation and unita- 
rization (when the scattering is seen in the dipole frame) require strong 
color fields in the hadron wavefunction, such that g J dx~ A'^ ~ 1. Qs 
is the critical transverse scale at which this strong field condition begins 
to be satisfied. 

What is specific to the present MV model (and, more generally, to 
any approximation in which the color sources are only weakly correlated 
with each other, like the gluonic sources in the DGLAP approximation 
l ^ l l'^ l) ig that the dipole scatters independently off the color sources in 
the hadron (here, the valence quarks). This is best seen by noticing that 
Eq. 12.5111 can be rewritten as a Glauber formula : 

^ , . j 2 tv^Cf AxGN{x,l/rl) \ 

SAir^) - expj - a.r^ -^^ 1, (2.53) 

where xGjq{x,Q^) is the gluon distribution of a nucleon, and is given in 
the present approximation by the second line in Eq. 12.3111 with A ^ 1. 



'^The impact parameter dependence was omitted in writing Eq. 12.511 since trivial for 
the case of a homogeneous target. This will be reintroduced later, when needed. 
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As we shall see in the next section, the previous picture changes 
quite substantially after including quantum corrections, due to the fact 
that the evolution towards small x induces correlations among the color 
sources. As a result, not only the general Glauber formula 12.531 becomes 
inapplicable (the successive scatterings are not independent any longer), 
but even its linearized 'leading-twist' approximation, corresponding to 
a single scattering, fails to apply when l/rj_ is close enough to Qs, 
while still above it. This is the BFKL regime where 'higher-twists' effects 
appear already in the linear evolution. 

Most interestingly, we shall see that gluon saturation at small x holds 
independently of the non-linear effects in the classical Yang-Mills equa- 
tions. Rather, this is the consequence of the correlations among the color 
sources induced by non-linear effects m the quantum evolution. These 
same correlations will be shown to ensure color screening already over the 
perturbative scale 1/Qs ^ ^/^QCDi which thus eliminates the infrared 
sensitivity of the classical MV model to the non-perturbative physics of 
confinement (see, e.g., Eqs. 12.3911 or 12.51ll l. 

3. The Quantum Evolution of the Color Glass Condensate 

In this section, we shall explain how to construct the CGC effective 
theory at small x by integrating out the gluons with x' > x in pertur- 
bation theory, in the presence of high density effects. The central result 
of this analysis will be a renormalization group equation for the weight 
function Wf^+ [p] in Eq. 12.71 , which generalizes the BFKL equation by 
including non-linear effects, and has important physical consequences 
among which gluon saturation. 

3.1. The BFKL evolution and its small-x problem 

Within perturbative QCD, the enhancement of the gluon distribution at 
small X proceeds via the gluon cascades depicted in Fig. 1101 Fig. llOl a 
shows the direct emission of a soft gluon with longitudinal momentum 
fc+ = xP+ < P+ by a fast moving parton (say, a valence quark) with 
= xqP^ and 1 > xq ^ x. Fig. llUl b displays the lowest-order radiative 
correction"^ which is of the order (with Qs = asNc/ir) 



relative to the tree-level process in Fig. llOl a. This correction is enhanced 
by the large rapidity interval At = ln(xQ/x) available for the emission 
of the additional gluon. 




(3.1) 



"At the same level of accuracy, a complete calculation must include also the appro- 
priate virtual corrections (self-energy and vertex renormalization); but for the present, 
qualitative purposes, it is sufficient to consider the real gluon emission. 
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P' 




Fig. 10. a) Small— X gluon emission by a fast parton; b) the lowest-order radiative cor- 
rection; c) a gluon cascade. 

A similar enhancement holds for the gluon cascade in Fig. llUl c. in 
which the succesive gluons are strongly ordered in longitudinal momenta: 

^ 3> P2 ^ ■ ■ • 3> Pn 3> fc^ . This gives a contribution of relative 
order 

l(..ln|l)", (3.2) 

where the factorial comes from the ordering in p"*". Clearly, when x is so 
small that ln(xo/x) ~ l/Ss, all such quantum "corrections" become of 
order one, and must be resummed for consistency. A calculation which 
includes effects of order {as ln(l/x))"' to all orders in n is said to be valid 
to "leading logarithmic accuracy" (LLA). 

The gluon cascades in Fig. IIUI contribute all to the production of 
(virtual) gluons with longitudinal fraction x. Thus, by resumming these 
cascades, one can compute the number of such gluons per unit rapidity, 
i.e., the gluon distribution (12.2511 . One can recognize in Eqs. 13.111 - 13.21 
the expansion of an exponential. Therefore: 

= xG{x,Q^) ~ g"""'^ — ^ (3.3) 

dr 

with uj a pure number. We have tacitly assumed that all the gluons in 
the cascade have transverse momenta of the same order, namely of order 
Q. A more refined treatment, based on the BFKL equation^-, allows one 
to compute uj and specifies the Q^-dependence, and also the subleading 
r-dependence (beyond the exponential behaviour shown in Eq. 13.311 ') of 
the gluon distribution. 

To describe the effects of the BFKL evolution in more detail, it is 
instructive to consider the dipole-hadron scattering introduced in Sect. 
12.71 With increasing energy, the gluon fields change in the hadron wave- 
function, and therefore so does also the cross-section for the dipole which 
couples to these fields (cf. Eq. 12.491 '). Specifically, to LLA, and in the 
linear regime where one can neglect multiple scattering, the amplitude 
■^rix±,y±) = 1 - 5'r(a::_L,y_L) = Hxy obeys to: 

^Mxy =as -2 {j^ccz + M.y ~Mxy], (3.4) 

dr " J 2tt ixj_- z^)^{yj_- z±y i tfj. v / 
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6 The Color Glass Condensate and High Energy Scattering in QCD 

which is the coordinate form of the BFKL equation!^ 



Fast hadron \ ^ /\/\/^ \ fast partoiis (p) 




small-x gluoiis (A[p]) 

Dipole 

Fig. 11. Deep inelastic scattering in the dipole frame 

The physical interpretation of this equation depends upon the 
Lorentz frame that we choose to visualize the process. When using the 
dipole frame of Sect. 12.71 the quantum evolution is put entirely in the 
wavefunction of the hadron (which is boosted to higher and higher ener- 
gies with increasing r), while the dipole remains a simple qq pair. In this 
frame, Eq. (13.411 describes the dipole-hadron scattering as the exchange 
of a BFKL ladder; the dipole couples to the last gluon (with the smallest 
value of x) in a gluon cascade which develops fully inside the hadron. 
(See also Fig.lTTI') 

Alternatively, by a change of frame, one can use the increase in the 
total energy to accelerate the dipole, and study the evolution of its wave- 
function with T. Under an increment dr such that asdr ~ 1, the dipole 
evolves by emitting one gluon (from either the quark or the antiquark), 
and the ensuing qqg state scatters off the hadronic target. It is conve- 
nient (although not necessary) to view this evolved state in the large- A^c 
limit, in which the radiated gluon is effectively replaced by a qq pair in 
a color octet state. Then, the evolution looks like the splitting of the 
original dipole into two new dipoles, each of them made of a quark (or 
antiquark) from the initial dipole and an antiquark (or a quark) from 
the emitted gluon. From this perspective, the various terms in Eq. 13.41 
have a simple interpretation: the quantity 

is the differential probability for the initial dipole {x_i,y±^ to decay 
into a pair of dipoles {x±,zj_) and {xj_,zj_), while Mxz and A/'zy are 
the amplitudes for the scattering between any one of these final dipoles 
and the target. Finally, the negative contribution proportional to —Nxy 
represents the decrease in the scattering amplitude of the original dipole 
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due to its dissociation (this term is necessary for the conservation of the 
probability) . 

One should mention here that this different perspective, in which the 
quantum evolution is put in the dipole and studied in the large-JVj limit, 
lies at the basis of an approach originally developed by Mueller'-*' — the 
'Color Dipole approach' — , in which the wave function of a very energetic 
color dipole (an "onium") is constructed in the BFKL approximation. 
In this approach, the scattering between two "onia" (physically, this 
corresponds, e.g., to the 7*7* scattering) can be treated as the product 
of the number of dipoles in each onium times the dipole-dipole cross- 
section. The ensuing scattering amplitude has been shown to satisfy 
equation 13.41 . 

The BFKL equation (EH 

can be solved by standard techniques'^. 
At high energy, otsT ^ ln(l/r^AQp^), and for a homogeneous target 
(e.g., for a large nucleus, and impact parameters near the center of the 
target), the solution reads (with r±_ = x±_ — y±^) 

Mr{r^) ^ I exp —-^^ Ly (3.6) 

where the reference scale Qq is introduced by the initial conditions at 
low energy, and thus is of order Aqqjj (for a nucleus, this carries the 
dependence on A). Furthermore, = 41n2 « 2.77 and (3 = 28^(3) « 
33.67. In writing Eq. we have assumed a fixed coupling as, as 

appropriate at leading-order BFKL accuracy. The modifications due to 
the running of the coupling will be discussed in Sect. ??. 

Eq. 13.61 exhibits two essential features of the BFKL approximation, 
which eventually provoke its failure in the high energy limit : 

(a) Violation of the unitarity bound : The solution 13.61 increases 
exponentially with r, that is, as a power of the energy. At high energy, 
such a behaviour violates both the unitarity bound A/V(r^,&j_) < 1 on 
the scattering amplitude at fixed impact parameter, and the Froissart 
bound o-(jipoi(,(s) < (tq In^ s'^on the total dipole cross-section 12.481 . 

(b) Infrared diffusion : When seen as a function of In {^\Q'q^, the 
solution 13.61 shows a diffusive behaviour, with diffusion 'time' propor- 
tional to r. With increasing energy, the typical transverse momenta car- 
ried by the gluons within the BFKL ladder diffuse away from the hard ex- 
ternal scale l/r\^ ^QCD' ''■'i*^ eventually enters the non-perturbative 
region at k\ ^Qcr>' which contradicts the use of perturbation theory. 

The first difficulty is similar to that with the linear approximation 
to the scattering amplitude in the MV model, Eq. 12.521 . which violates 
unitarity" for very large A. As in that case, we expect unitarity to be 



"We mean here the unitarity of the scattering amplitude at fixed impact parameter. The 
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restored by multiple scatterings, i.e., by keeping terms of all orders in 
gA+ in the Wilson lines in Eqs. JTI^ - dTtnt . Both Eqs. llH^ and 
correspond to a single scattering approximation-they are obtained 
by retaining terms which are, at most, quadratic in gA~^ in the expansion 
of the Wilson lines in Eq. 

However, unlike in the MV model, at small x we expect an additional 
source of non-linearities which arise from the interactions among the glu- 
onic sources. Such interactions lead to the fusion of gluons from different 
parton cascades ("gluon recombination"), a phenomenon which should 
tame the rapid growth of the number of partons. Thus saturation arises 
from the competition of two effects: the growth of the gluonic density 
due to radiation and its depletion due to recombination effects ' — ' — ! 
Since now the dominant color sources are themselves gluons, the satura- 
tion we speak of here refers simultaneously to the source and the fields 
radiated by them. This is because what we call "sources" and "radiated 
fields" is only relative, as it depends upon the scale — xP^ at which 
we consider the effective theory. 

Both the multiple scattering and gluon recombination mechanisms 
are illustrated in Fig. 1111 We expect both mechanisms to become impor- 
tant at the same scale, Q'^{t,A), which is the critical gluon density at 
which the non-linear effects become of order one. This saturation scale 
is also the typical transverse momentum of the gluons in the hadron 
wavefunction at small x. As anticipated in the Introduction, and will 
be verified explicitly in what follows, this scale increases rapidly with r 
and A. The emergence of such a hard intrinsic momentum scale can also 
solve the 'infrared diffusion' problem of the BFKL approximation and 
therefore restore the applicability of perturbation theory to high energy 
processes in QCD. 

In the dipole frame, gluon recombination is seen as the merging of 
two gluon cascades, as illustrated in the r.h.s. of Fig. 1111 It is interesting 
to see this process also from the boosted frame in which the quantum 
evolution proceeds via the dissociation of the incoming dipole into two 
dipoles. After the boost, the final gluon in the cascade in Fig. 111! — the 
one which couples to the dipole — gets incorporated within the dipole 
wavefunction, so the merging of two cascades now happens inside the 
dipole. Thus, from the boosted frame, the non-linear process is seen 
as the simultaneous scattering of the two final dipoles off the hadronic 
target. 

These considerations suggest the following simple equation which 
generalizes Eq. (13.411 by taking non-linear effects into account: 




X {N'xz +Mz 



zy 



xy 



MxzMzy} ■ 



(3.7) 



discussion of the total cross-section is more involved, and deferred to Sect. 4.4. 
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This is the equat ion derived by Kovchegov within Mueller's 'Color 
Dipole approach' BHl The main assumption used in its derivation was 
the fact the two final dipoles scatter independently off the target; it is 
correct only in the large- A^c limit. Deriving this equation, (tog ethe r with 
its generalizations to finite Nc originally obtained by Balitsky 1^ within 
a different formalism) within the framework of the CGC effective theory, 
will be a main objective of the forthcoming developments in this section. 

3.2. Non-linear evolution for the CGC 

The CGC provides a natural framework for the description of the non- 
linear effects in the quantum evolution towards small x, and of the 
phenomenon of saturation. The main observation is that, to LLA, all 
the quantum corrections described previously — both the exponentially 
developing BFKL cascade, and gluon recombination which tames this 
rapid growth — can be incorporated into a change of the classical color 
charge and its correlations, namely, into a renormalization of the weight 
function Wa+ [p] in Eq. 112.711 . 

To see this at an intuitive level, let us reconsider the first radiative 
correction, the one-gluon emission in Fig. llUl b. and note that, to LLA, 
the typical contributions to the integral in Eq. 13. H come from momenta 

such that p"*" ^ p'f 3> fc"*". That is, the condition of separation of 
scales is indeed satisfied for the intermediate gluon with momentum p^ 
to be treated as a 'frozen' color source for the final gluon with momentum 
fc"*". The effect of this quantum correction is therefore simply to renor- 
malize the effective color source at scale , as pictorially illustrated in 
Fig. EH 




Fig. 12. Effective color source after including the lowest-order radiative correction. 

By iterating this argument, it is quite clear that a whole BFKL 
cascade (see Fig. llOl c'l can be included in the definition of the classical 
color source at the scale A"*" — xP^ of interest. It is furthermore clear 
that the fusion between two gluon cascades, as illustrated in the l.h.s. 
of Fig. 1131 can be represented in the CGC theory as a non-linear effect 
in the classical dynamics of the color fields generated by this effective 
source (see the r.h.s. of Fig. 1131 . 

But non-linear effects are important also in the quantum evolution, 
and actually interfere with it, as illustrated in Fig. 1141 Fig. 1141 a is an 
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Fig. 13. The fusion of two gluon cascades and its interpretation in the CGC theory. 

immediate generalization of the one-gluon emission in Fig. llUl b. It is 
clear that what is renormalized by the scattering off the "semi-fast" 
(A^ ^ ^ k^) quantum fluctuation is the classical field A^[p] at 
scale A^, which in turn is non-linear in p. (The Feynman rules for eval- 
uating diagrams like those in Fig. 1141 and also the present discussion, 
are adapted to the LC gauge = 0, in which the quantum effective 
theory is written, cf. Eq. H2.7|l .'l Fig. 1141 b shows an additional source of 
non-linearity, arising from the propagation of the radiated gluon in the 
classical 'background' field .4*[p]. If A'*' = xP^ is small enough (x ^ 1), 
the classical field is very strong, ~ 1/g, and gluon rescatterings must 
be included to all orders in ■ The diagrams in both Figs. 1141 a and b 
can be taken into account as the cut of the diagram in Fig. 1141 c. The clas- 
sical field that enters the vertices is the fully non- linear solution A^[p] 
constructed in Sect. 12.41 and the propagator of the quantum gluon is 
computed to all orders in this background field, the resummation indi- 
cated here by a blob. 




a) b) c) 



Fig. 14. Some typical non-linear effects in the quantum evolution 



The diagram in Fig. 1141 c is manifestly a quantum correction to 
the 2-point function of the gauge fields at scale , and is of order 
as ln(A^/fc^"). Thus, for this to be computable in perturbation theory, 
the separation of scales between A"^ and must not be too large: 
A^ S> , but asln(A^/fc^) <^ 1. The quantum modes must be in- 
tegrated out in layers of , within a renormalization group procedure 
I— . At each step in this procedure, one has to perform a one-loop 
quantum calculation, but with the exact background field propagator 
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for the "semi-fast" gluons (the quantum gluons that are integrated out 
in that particular step) . Such an all-order inclusion of the classical field 
effects permits one to resum not only the large energy logarithms, namely 
the terms ~ (a^ In l/a:)", but also the dominant high density effects- the 
non-linear effects (like gluon recombination) which become of order one 
at saturation. 

The corresponding analysis is technicall y quite involyed and has been 
described in detail in previous publications l ^'^ l ^" ! Here we shall pre sent 
only the final results and their consequences. As shown in Ref. the 
condition that the new correlations induced by integrating out quan- 
tum fluctuations be reproduced by the CGC effective theory leads to a 
functional renorraalization group equation (RGE) for the weight function 
Wj^+ [p] = Wt [p] , which is most succinctly written as 

j;;^,Xa,i^.,y.M-^WAp], (3.8) 

in notations that we shall shortly explain. Ear ly vers ions of this equation 
can be found in the pion eerin g works of Refs.' — A formally similar, 
and physically equivalent!^, functional evolution equation has been ob- 
tained by Weigert^^, within a different formalism We shall discuss 
this latter approach in the next subsection. 

Let us now discuss the meaning and structure of the terms in 
Eq. I|8.8^ . The rapidity variable r = ln(l/a;) = ln(P+/A''") indicates 
the dependence of the effective theory upon the separation scale. This is 
convenient since (as illustrated by the BFKL evolution discussed in Sect. 
13. in . r is the natural "evolution time". The contribution of the quantum 
modes within a small layer in p"*" (say A^ > > fc^) is proportional 
to the rapidity extent Ar = ln(A~''/A;~'') of that layer. 

The kernel x[p\ is a positive definite non-linear functional of p 
(Eq. 13.81 is a diffusion equation) and is highly non-local in both longi- 
tudinal and transverse coordinates. The non-linearity in p and the non- 
locality in x~ are strongly correlated, since they have a common origin: 
x[p\ depends upon p via the Wilson lines (12.2411 . Physically, x[p\d''^ is 
the charge-charge correlator induced when integrating out quantum glu- 
ons within a rapidity interval dr and in the pr esence of a classical color 
charge distribution with density pa(x) Technically, this is com- 

puted by evaluating the diagram in Fig. 1141 c with the Feynman rules 
in Eq. (12.711 . Note that, in addition to this real-gluon emission diagram, 
there are also virtu al (self- energy and vertex) corrections which must be 
similarly computed' — Such corrections are already included in the 
RGE (13.811 ._where they correspond to the functional derivative of the 
kernel xM*^- 

The argument p°(x^) of the functional derivatives in Eq. (13.811 de- 
notes the color charge density p°'{x~ ,x^) at x~ = x^ . The color source 
generated by the quantum evolution up to rapidity r has support within 
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a limited interval in x~ , namely at < x~ < x^ , with x^ oc e'^. The 
reason this is so follows from the uncertainty principle: since the classical 
source at rapidity r is obtained by integrating out quantum modes with 
large longitudinal momenta p"*" ^ = e~^P^, it must be localized 
near x~ = 0, within a distance Ax~ ~ ^^^q (with x^ = 1/P~^). 

However, eq. I|8.8^ shows that the correlation between the quantum 
evolution in r and the longitudinal distribution of the resulting color 
source is even stronger. When the rapidity is further increased, say from 
r to T + dr, the additional contribution to the color source which is 
generated in this way has no overlap in a::~ with the original source at 
rapidity r. Instead, this new contribution makes a new layer in a::" , which 
is located between x^ and x~_^_^^. This is why the functional derivatives 

in eq. (13.811 involve just the color source pr = p{xr) in this outermost 
layer. 

This correlation is most simply formulated if one uses the space-time 
rapidity y, 

y = ln(a;"/a;Q ), Xq = , — oo < y < oo , (3.9) 

to indicate the longitudinal coordinate of a field. For example, 

Py{xj^) = x^p"'{x^,x^) for x^=Xy=XQe^, 
j dy py{xx) = j dx^ p'\x^ ,xx), (3.10) 

and similarly for a, Eq. 12.1811 . or any other field. Eq. TTM can be 

rewritten as : 



V\x^) = Pexp|igydyay(a::x)i°| 



(3.11) 



The previous discussion shows that the space-time rapidity y of a given 
layer in p is identical to the usual (momentum) rapidity of the fast gluons 
that have produced that layer. In particular, the color source created 
by the quantum evolution up to r has support at space-time rapidities 
y < r (in agreement with the simple argument based on the uncertainty 
principle). Formally: 

Wr[p\ OC 5r[p\, (3.12) 

where 5t[p\ is a (5-functional enforcing that py = for any y > r. As we 
shall see, this constraint is important because the r-dependence of the 
observables in the effective theory comes precisely from the upper limit 
on the longitudinal support of p. The Color Glass evolves by expanding 

in y. 

Since the Wilson lines 13.1111 and many interesting quantities (like 
the ^-matrix element (12.491 . or the gluon distribution 12.3211 ) are more 
directly expressed in terms of the COV-gauge field oiy{x±_), rather than 
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the color charge p, it is often preferable to use the 'a-representation', 
whose weight function Wr[a] = Wt[p = — V^^a] satisfies the following 
RGE, obtained after a change of variables in Eq. H3.8|l ; 

dWr[a] 1 5 „ 1 SWr 

We use compact notations in which repeated color indices (and coordi- 
nates) are understood to be summed (integrated) over. The relation of 
the kernel here to that in Eq. 113.811 is 

Jz±,u± Vj_ V_|_ 

The analysis in Ref.l^ yields (see alsol^ : 



ah ( 



IT J (27r)2 (a:;j_ - 2_L)2(y^ - 2^)2 

X {l + V^Vy ~ - V^VyY', (3.14) 

with Vx = V^{xj_), etc. This is real and symmetric {ri'^'^ {xj_,y^) — 
ri'"^{y±,xj^)), and also positive definite, as anticipated, since: 

{l + Vhy - V^Vz - V.^VyY'' = (1 - V^Vx)ca{l - V^Vy),t, (3.15) 

and the color matrix 1 — V^Vx is hermitian. The transverse kernel in 
Eq. 13.1411 is similar to the 'dipole kernel' in the BFKL equation 13.41 . 
Their relation will be discussed in the next subsection. 

The RGE in eq. 13.1311 [or eq. 13.81 1 has the structure of a Fokker- 
Planck equation. In the CGC formalism, the quantum evolution towards 
small X is a random walk in the space of Wilson lines The random 
process is one by which the Wilson lines are built. The physical random 
variable in the evolution is the elementary contribution a'^{x^_) to the 
classical field in the hadron arising from integrating out quantum fluc- 
tuations in the rapidity strip [r, r -I- dr] . Such a contribution changes the 
Wilson lines according to: 

Ul+,^x^) = e'3''^<(--)^" C/t(xj_) , (3.16) 

whose iteration defines a path in the space of the U fields. (This path 
is unam bigu ously defined only after discretizing the rapidity variable; 
see Ref. I^for details.) By exploiting this representation, an exact but 
formal s olut ion to Eq. 13.131 has been constructed in the form of a path 
integral This random walk can be equivalently reformulated as a 
Langevin equation!^, a formulation which is better suited for numerical 
simulations on a two-dimensional lattice. 
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3.3. The Balitsky-Kovchegov equation 

In addition to the numerical simulations on a lattice, Eq. 13.131 can be 
made tractable via two strategies. Both involve some approximations. 
The first strategy consists in using this functional equation to deduce 
ordinary differential equations for quantities of interest. Because of the 
non-hnearity of Eq. 13.131 . the ensuing equations will generally not be 
closed, but rather form an infinite hierarchy of coupled equations. Nev- 
ertheless some progress can be made in various approximations, partic- 
ularly in the large Nc limit, where we shall see that a closed equation 
emerges. The other strategy, to be developed in the next subsection, is 
to search directly for approximate solutions to the functional equation 
13.131 for the weight function. 

If {0[a] )t is any observable which can be computed as an average 
over a (cf. eq. lIT^ ) : 

{0[a]}r = [ VaO[a]Wr[a], (3.17) 



then its evolution with r is governed by the following equation: 
-{0[a])r= / VaO[a] ^ 



d 



2 5a?(xi) 5a%{y^) 

where, in writing the second line, we have used eq. 13.131 and then in- 
tegrated twice by parts within the functional integral over a. The last 
expression involves the functional derivative of r][a\, which is easily com- 
puted by using [with Sxy = ^^^•'(xj^ ~ J/±)] ■ 

gg^^.,..,T«^t(,^), ^^^~^O^.yVi.,)T^. (3.19) 

Here Oy = for y > r, cf. Eq. lISTT^ . 

The 2-point function Sr(a;x,J/±) of the Wilson lines, Eq. 12.49^ 
(which physically represents the ^-matrix element for dipole-hadron 
scattering) ca n straightforwardly be computed using repeatedly 
Eq. lITT^ (see 1^ for details): 



X l^NMV^Vy) - tT{V^Vz)tr{vhy)) ^ ■ (3.20) 

This equation was originally derived by Balitsky within a formal- 
ism based on the evolution of observables (in high-energy dipole-hadron 
scattering) which are built from Wilson lines. This is similar in spirit 
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to the 'Color Dipole approach' by Mueller but it is not restricted to 
the large-A'^c limit. It is better suited for an asymmetric collision, like 
that between an "onium" (= a high-energy dipole) and a dense hadronic 
target, like a nucleus. In this respect, Balitsky's formalism is closer to 
the CGC formalism, where the focus is fully on the target wavefunction. 

As anticipated, the above equation is not closed; It relates the 2- 
point function to the 4-point function {tr(14f Vz)tr(Vz^Vy)). Physically, it 
is so since, except at large Nc, the qqg system formed after radiating one 
gluon from the original dipole is not exactly the same as a system of two 
dipoles (recall the discussion around Eq. I|8.5^ ). On e can similarly derive 
an evolution equation for the 4-point function but this will in turn 
couple the 4-point function to a 6-point function, and so on. Eq. 13.201 
is merel y the first in an infinite hierarchy of coupled equations . In 
Ref.l^, Weigert managed to reformulate Balitsky's hierarchy as a single 
functional evolution equation for the generating functional of the n-point 
functions of the Wilson lines. As shown in Ref. , Weigert's equation 
is equivalent to the RGE 13.131 as far as the correlations of the Wilson 
lines are concerned. Mo re r ecently, Mueller used a similar approach to 
give a simple derivation!^ for Eq. 13.1311 . 

A closed equation can still be obtained in the large Nc limit in which 
the 4-point function in eq. 13.21)1 factorizes: 

{tr{V^V,)tr{V,^Vy))^ (^tr{V^Vz)'j ^ {tr{V,Wy)) ^ for Nc oo. 

Then Eq. 13.2U1 reduces to a closed equation for the 2-point func- 
tion, which, when rewritten in terms of the scattering amplitude 
■^r{x±,y±) = j^(tr(l — VxVy))T, is recognized as the Kovchegov equa- 
tion 13.71 An early v ersio n of this equation has been pr opos ed by 
Gribov, Levin and Ryskin^^, and proven by Mueller and Qiu^^in the 
'double-logarithmic approximation'. More rece ntly, Braun has rederived 
Eq. 13.71 by directly resumming 'fan' diagrams'*—'. Following the recent 
literature, we shall refer to Eq. 13.71 as the "Balitsky-Kovchegov (BK) 
equation" 

Clearly, in the weak scattering approximation A/'r(7'_L) ^ 1 (which 
corresponds to a very small dipole, or, equivalently, to a relatively low 
gluon density in the hadronic target), the non-linear term can be ne- 
glected in the r.h.s. of Eq. 13.71 . which then reduces to the BFKL equa- 
tion 13.41 . But in general, the feedback provided by this non-linear term 
ensures that the solution J^T{r±) to Eq. 13.71 respects the unitarity 
bound A/'t('"j_) < 1- Thus, Eq. 13.71 is a simple QCD-based non-linear 
equation consistent with unitarity. This explains the large interest in this 
equation in the recent literature, with important progress to wards its res- 
olution via both analytic | 53 . 65.46] numerical methods I^^I^^EEi 
The conclusions reached in this approach are equivalent to those ob- 
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tained from direct investigations of the RGE l|3.13|l l ^^ l ^^ l and will be 
described in the next sections. 



3.4. Saturation momentum and geometric scaling 




Fig. 15. The solution to the BK equation as a function of r±_ for two values of r. 



The solution to Eq. 1)3. 7|l is shown qualitatively in Fig. 1151 which 
displays Mrirj^) as a function of r_L = — yj_ for two different rapidi- 
ties. The scattering amplitude vanishes as rj_ 0, as it should from its 
definition JVT{r±) — -^{tr(l — V^Vi/))t, and the fact that the Wilson 
lines are unitary matrices. For small rj_, A/'r(fj_) remains small ("color 
transparency"), and is well approximated by the BFKL solution 13.61 1. 
For large r_L, it approaches the unitarity bound A/t('"_l) = 1. The transi- 
tion between "color transparency" at small rj_ and "blackness" at large 
r_L takes place at a characteristic value of r± that we shall identify with 
the saturation length l/Qsij). More precisely, we shall define Qsij) by 
the following convention: 

NT{ri) = 1/2 for r^^l/Qs{T). (3.21) 

As shown in Fig. 1151 this saturation length decreases with r. 

The result illustrated in Fig. ll5l clearlv shows that the non- linear BK 
equation (13.711 solves the unitarity problem of the BFKL equation (at 
least, at a fixed impact parameter; see also Sect. 4.4 below). Moreover, 
the emergence of an intrinsic "saturation scale" Qs{t), which increases 
with r, also solves the problem of "infrared diff usio n" , as convincingly 
demonstrated by the numerical analysis in Ref. In addition to the 
numerical studies in Refs. I^ i i^*^ - ^^ - the behaviour shown in Fig. 1151 
is supported also by analytic investigations focusing on qual itative fea- 
tures like the energy dependence of the saturation scale |15|4fci|94 J44| 
the "geometric scaling" behavio ur or the approach towards 

the blackness with increasing r_L We shall describe here 
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some of these analytical studies, whose results follow from general argu- 
ments, such as the validity of the BFKL dynamics at small r^. and the 
emergence (via non-linear effects) of an intrinsic momentum scale, the 
saturation momentum Qs{t). 

Consider first the calculation of the saturation scale. Even though the 
BFKL solution 1)3. 6|l is valid only at small r_L, well below the saturation 
length, it is nevertheless possible to compute the energy dependence of 
the saturation scale by extrapolating Eq. 13.611 up to r± ^ 1/Qs{r) 
and then imposing the saturation condition 13.2111 ' ' '. More precisely, 
we shall see shortly that Q'^^t) is increasing exponentially with r. The 
BFKL computation alluded to above should then correctly reproduce the 
value of the exponenl^but not necessarily the (slowly varying) prefactor 
as well (see however**-^. 

Towards this end, the solution of the BFKL equation in Eq. 113. till is 
first re-written as (p = ln(l/r^Qo))- 

AAr(r^)~expLQsr-^--|^|, (3.22) 

where we have kept only the dominant terms in the regime p 3> 1 and 
asT ^ 1, with ctsT 3> p. Since the BFKL equation is now seen as just an 
approximation to more general non-linear equations like the BK equation 
Ij3.7|l . the solution 13.22^ is acceptable only as long as r±_ ^ 1/Qs{t), or 
p > ps(r), with p.(r) = In(g2(r)/Q§). 

If nevertheless extrapolated down to p ~ Ps{t), the saturation con- 
dition 13.211 amounts to the vanishing of the exponent in Eq. 13.221 1. 
This giv es a se cond-order algebraic equation for ps(r) with the physical 
solution 

Q?(r) = qI e"""^ , c=[-p + y/(3{/3 + 8Lu)] /2 = 4.84... (3.23) 

This estimate is consistent w ith the nu merical solutions to the BK equa- 
tion, which have found c ^ 4 165 1 66 1 6 51 j^^^ with the phenomenology 
of DIS at HE RA, which suggests rather a significantly lower value for the 
exponent (see Sect. 4.3 below), namely, A ~ 0.3 instead of cqs ~ 1. 

The factor in front of the exponenti al in eq. 13.231 is not under 
control in present approximations. In Ref. a more refined treatment 
was proposed, where the BFKL equation was solved with an absorbtive 
boundary condition at rj_ ~ 1/Qs{t), and a weak dependence on r for 
this prefactor was obtained. It would be interesting to test their results 
against more accurate numerical solutions to Eq. 13.71 . 

It is also interesting to study the behaviour of the scattering ampli- 
tude 13.221 for rj^ below but relatively close to l/Qs(r)-for p slightly 
above Ps{t). Since: 

P - 1- 3^ = P^(-) + 1- ;ra(7y - + SP , (3.24) 
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a simple calculation yields: 

AA,(r^)^exp{-7<5p-;^|^|, (3.25) 

where 7 = 1/2 + c//3 ~ 0.64. Eq. H3.25|l suggests a remarkable simpli- 
fication: Assume that is sufficiently close to 1/Qs{t) (although still 
below it) for Sp/dsT ^ 1- Then, the second term in the exponent can 
be neglected compared to the first one, and we are left with 

Ur{r^) ^ (rlQl{r)y , (3.26) 



which shows geometric scaling' ' ' : it depends upon the two kinemat- 
ical variables rj_ and r only via the combination »"^(3s(t). Numerically, 
such a scaling behaviour has been seen in the solutions to the BK equa- 
tionB5 68 

Since Eq. 13.2611 is the first term in an expansion in powers of Sp/ps, 
with ps = cusT, this approximation is correct for 

^<^-^r^)«'^^sr. (3.27) 



The condition on the left, rj^ < 1/Qs{t), ensures we are still in a linear 

M 



regime. For a dipole transverse resolutio n = this condition 

translates to the following scaling window^ 



Q?(r) « « (3.28) 

Since Qq ~ Aqcd and Qsij) 3> Aqcd for sufficiently large r, the 
upper boundary of this scaling window is rather large. In particular, 
it is much larger than the saturation scale itself: Qs(t)/Qo ^ Qs{t). 
Remarkably, as a consequence of saturation, knowledge of an intrinsic 
momentum scale is propagated through the linear evolution equations 
up to relatively large values of well outside the saturation regime. 
This is especially interesting since such values of are large enough for 
perturbation theory to be fully trustworthy. A prop ertv like 113.261 can 
be and has been tested against the experimental data^^^ ^^ l ^^ l We^shall 
return to phenomenological aspects of geometric scaling in Sect. 4.3. 

The previous results are obtained from the leading order BFKL equa- 
tion. Recently, there has been some progress in inclu ding ne xt-to-leading- 
order as effects in the physics of saturation. In Refs.^ ' \ this was done 
heuristically, by simply replacing the fixed coupling as in the BFKL 
equation by the one-loop running coupling of QCD with the running 
scale set by the saturation momentum : as — > Qs(Qs(''"))j with 
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The only modification due to the running coupling is in the functional 
form of the saturation scale, whose growth with r becomes somewhat 
milder (tq is an arbitrary constant, and c is the same number as in in 
Eq. (Ei^): 

where the overall scale is now set by Aqcd rather than the initial scale 
Qo . All the previous results on geometric scaling (the scaling law H3.26^ , 
including the value of the "anomalous dimension" 7, and its range of 
validity 13.281 ') remain unchanged, except for the expression for the sat- 
uration scale entering these results. 

For nuclei, Eq. (13.301 has an intriguing consequence for the depen- 
dence of the saturation scale upon the atomic number ^1^. Assume an 
initial condition of the MV type at r = 0, Qs{t = 0, A) = Qs(^) ~ 
A^^^lnA. For fixed coupling BFKL evolution, where Eq. 13.231 applies, 
this initial condition identifies the hitherto unspecified 'initial' scale Qo 
with the MV saturation scale, 

Qs(r,^) = (3^(^)6*^"°^ (fixed coupling) , (3.31) 

which preserves the yl-dependence of the initial condition at any later 
'time' r: 0?(t,^) A^/^lnA. 

For the running coupling BFKL evolution case where Eq. 13.301 
applies (with tq fixed by the initial condition as 2feo crp — 
[ln((3s(A)/AQ(-.j-))]2), one obtains a very different yl-dependence at small 

and large r, respectively 

i) At relatively small energies, such that r <C tq ~ In^ A^^^, 

QI{t,A) « Q2(A)e'="=('3'(^»^, (3.32) 

which is the 'fixed-coupling'-like behaviour, with the as in the exponent 
being the running coupling 13.291 evaluated at the initial saturation 
scale. 

ii) At higher energies, r ^ In^ A^^"^, one obtains 



)?(., A) . A|cD exp | ( In %^ ) }, (3.33) 



^^QCD 



2 ' 



which, for very large r, is nearly independent of A. 

In Ref. ' — ', which did not consider the A-dependence, Triantafyl- 
lopoulos presented a complete computation of the NLO effects on the en- 
ergy dependence of the sat uration scale. Recall that the NLO corrections 
to the BFKL equation turn out to be anomalously large and require 
resummation to obtain sensible results. Ref. used the RG-improved 
resummation scheme of Ciafaloni, Colferai, and Salam^^and found that. 
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although <3s(t) is in general a more complicated function than the simple 
exponential H3.23|l . it can nevertheless be represented as such for a rather 
wide range of rapidities (including those of phenomenological interest). 
Specifically, if one defines A(r) = dln(Q^/A^)/dr, then A(r) turns out 
to be a very slowly decreasing function, with A(r = 5 — 9) ~ 0.30 — 0.29. 
Remarkably, this value is also favoured by the current phenomenology at 

AO 

HERA (and Sect. 4.3 below). It would be very interesting to compute 
the j4-dependence within this fully NLO formalism. 

3.5. Gluon saturation and perturbative color neutrality 

In this section, we return to the RGE 13.13II . which describes the evolu- 
tion of the hadron wavefunction as a whole, and construct approximate 
solutions to it. As usual, these approximations depend upon the trans- 
verse resolution scale at which correlations are measured. If is 
large enough (Q^ ^ Qs (''"))) we probe color sources with small trans- 
verse size, which do not overlap with each other. In this dilute regime, a 
description in terms of uncorrelated sources, such as the MV model, may 
be a good approximation. With increasing r at fixed (or, equivalently, 
with decreasing at fixed r), spatial correlations start to develop, ini- 
tially according to the linear BFKL evolution and then, once the density 
is high enough, according to the general non-linear RGE which predicts 
gluon saturation. A schematic map of the kinematical regimes for quan- 
tum evolution is shown in Fig. 1161 f see the discussion below for details). 

i) High— momentum regime (Q^ > Qs): Recovering BFKL 

In the dilute regime at 2> Q^, the color charge density is low, 
hence the classical field is weak, ga <^ 1. One can therefore keep only the 
linear term in the expansion of the Wilson lines in Eq. (13.131 1 in powers 
of ga : 



{aa{xj_) is the effective color field in the transverse plane.) The kernel 
ri[a] in Eq. (13.1411 then becomes quadratic in a{x±), and the RGE takes 
the generic form: 



where the new kernel IC is non-local in the transverse coordinates. Its 
explicit form is easily extracted from Eq. (13.141 using Eq. (IS331- Even m 
this dilute regime, the RGE is non-linear and the corresponding weight 
function Wr [a] is not a Gaussian. Nevertheless, compared to the general 



V'^{x±^) ^ I + igj dyay{xi} = 1 + iga{x±), 
1 - V^Vx ^ ig{x^ - x^)d'a{xj^). 



(3.34) 




(3.35) 
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51 



Linear 



InA^ Ink^ 



Parton Gas 




Fig. 16. A map of the quantum evolution in the r — fcj^ plane. 



RGE 13.131 . the evolution generated by Eq. 13.351 exhibits an impor- 
tant simplification: it does n ot m ix correlations {a(l)a{2) . . . a{n)) with 
different numbers n of fields Indeed, the quartic operator acting on 
Wt[q:] in the r.h.s. of Eq. 13.3511 is formally the same as a the second- 
quantized Hamiltonian for a non-relativistic many-body system, and is 
diagonal in the number of "particles" . 

Eq. 13.3511 provides a closed evolution equation for the 2-point func- 
tion, which is the BFKL equation. This equation is most commonly 
written for the charge-charge correlator i-tr(k±) oc {pa{k±)pa{—k±))T- 
In the linear regime, it is the same as the unintegrated gluon distribution 
12.281 . Specifically, if one defines p,'^{kj_) as the Fourier transform of : 



2, X _ (P°(^±)p°(j/jj)r a, ^- a, . 

fj'T{xj_,yi_) = -I ' ^ " J '^^Py^^^)' 

then Eqs. ETEl and imply, similar to Eq. lICT^ . 

<^,(fc^) ~ ^^igil for fci»Q.(r), 



(3.36) 



(3.37) 



while from the RGE 13.351 one obtains' 



dr 



d^PA 



■^tJ-T{k±) 



(3.38) 



1, 2008 20:44 



WSPC/Trim Size: 9in x 6in for Review Volume 



QGPS'finar 



2 The Color Glass Condensate and High Energy Scattering in QCD 

which is indeed the BFKL equation 

Note that fJ-rik^) corresponds to of the MV model, but, unlike 
the latter, it carries a non-trivial transverse momentum dependence, 
and a r dependence, both of which are obtained by solving Eq. I|3.38^ . 
In fact, given the formal similitude between Eqs. 13.3811 and 13. 4L it 
is clear that the corresponding solution for ^^{k±) can be obtained by 
replacing t^Qq k\/Q^ in Eq. 13.611 . This shows the expected rapid 
exponential growth with r and infrared diffusion as well. 

ii) Low— momentum regime (Q^ <^ Q^) : Saturation 

An external probe with low transverse resolution <^ couples 
mostly to the saturated gluons, which have momenta kj^ ^ Qs and oc- 
cupation numbers ~ 1/as. The corresponds classical fields are strong, 
ga{x±) ~ 1, so the Wilson lines 13.111 — which are complex exponentials 
built with these fields — oscillate around zero over a characteristic dis- 
tance ~ 1/Qs(t) in the transverse plane. This implies that Wilson lines 
which are separated by large distances ^ 1/Qs{t) are necessarily uncor- 
related (since their relative phases are random) . Thus, when studying the 
dynamics over large transverse separations rj_ 3> 1/Qs(r), it should be 
a good approximation to neglect the correlations of the Wilson lines (or, 
more generally, to treat them as small quantities ). This is the "random 
phase approximation" (RPA) introduced in Refs. I -^^ I ^^ I 

In this approximation, the RGE 13.811 simplifies drastically^^. Ne- 
glecting the Wilson lines, the kernel rj becomes independent of a, and 
the RGE reads in momentum space, 

dWr[a] _ 1 fd^k^ 1 S^Wrlp] (3 39) 



dr 2 J (27r)2 7rfc2_ Sa^{kj_)5a'i^{-kj_) ' 

Being quadratic, this equation can be immediately integrated 

[P] - AA.expl - ^ / dy Tg% ^°(^^)f(-^^) 1, (3.40) 



2 i J {2tv 



which for convenience has been written as a functional of py{kj_) — 
k]_a^{k^). Eq. ijTint is a low-momentum approximation: for a given 
rapidity y (with y < r). It is to be used only for modes kj_ < Qs(y). 

We see an interesting duality emerging at saturation: this strong field 
regime allows for a description in terms of a Gaussian weight function, 
like for a free theory. But even with such a Gaussian weight function, 
the CGC effective theory remains non-trivial, since the classical solution 
12.161 and the observables for high energy scattering, like 12.491 . are non- 
linear functionals of p. 
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Eq. (I3.4UII shows that the only non-trivial correlation of the color 
sources with small kj_ is the 2-point function, and reads: 

{Pyix±) Py'iy±))T = 5°''<5(y - yX-r - y) hixj_,yj_), 

Ay(fc^) ~ i fci , for « Qs{y). (3.41) 

This distribution is local in (space-time) rapidity y, and homogeneous 
in all the (longitudinal and transverse) coordinates. In the transverse 
plane, it is only a function of the relative coordinate x± — yj_, and in the 
longitudinal direction, it is independent of y. For a given kj_ <g; Qs{t), 
Eq. 13.411 applies only for y in the interval Ts(fcj_) < y < r, with Ts{k^) 
being the rapidity at which the saturation scale becomes equal to the 
momentum kj_ of interest : 

Ql(r) = ki for T = Tsikj^), (3.42) 

(see Fig. Il(il . It follows that the integrated quantity (cf. eq. 13.3611 ) : 

T 2 2 

Mr(fci)| = / rfy = (r-r^fc^))^, (fc^ « Q^r)), (3.43) 

I sat J TT V /TV 

which measures the density of saturated color sources (with given k^_) 
in the transverse plane, grows only linearly with t. This should be con- 
trasted with the exponential increase of the corresponding quantity at 
k±_ > Qs{t), obtained from the BFKL equation (lOSll . (In Ref. the 
result in Eq. 13.431 was obtained from a study of the BK equation.) 

We conclude that, at low momenta kj_ ^ Qs(t), the color charge 
density saturates, because of the non-linear effects in the quantum evo- 
lution In contras, in the MV model, the corresponding quantity /i^ 
increases indefinitely as a power of the atomic number A (the analog of 
the energy in the MV model), since the respective sources are uncorre- 
lated. Due to the quantum evolution, correlations arise already in the 
linear regime via the BFKL equation. However it is only after includ- 
ing the non-linear effects associated with gluon recombination that these 
correlations are such as to limit the growth of the color charge density 
with r. 

Since the color sources saturate with r, so does also the gluon dis- 
tribution radiated by these sources. This is obvious if the linear ap- 
proximation^, Eq. 13.371 . is used to compute this distribution. Namely, 
Eqs. lIT?! and imply: 

^ = ^ In ^(l) , (3.44) 



fRecall that Eq. was obtained by using the linearized solution to the classical 

EDM, i.e., Eq. VTJM . 
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where in writing the second equahty we have used Eq. 13.2311 for the 
saturation scale^, together with the definition 13.421 of Ts{kj_). A more 
careful calculation, based on the non-linear solution l|2.1()ll . shows that 
the correct answ er for ifir(k \ ) differs from Eq. 13.4411 only by a numer- 
ical factor 5^1 l ^^ l ^^ l This factor cannot be computed in the present 
approximations, as it is sensitive to the transition regime at fcj_ ~ Qs{t). 
This overall factor will not affect the salient features of Eq. 13.441 , namely 
its dependence upon the energy and the transverse momentum. 

^4 




Fig. 17. The gluon phase-density (pT{k±) in the effective theory plotted as a function 
of fex for two values of t = ln(l/x). 



In addition to t -saturation (the linear increase with r), Eq. 13.44t 
shows also k^-saturation-the kj_ spectrum is only logarithmic in at 
low momenta. Recall that, in the classical MV model, a similar spectrum 
emerged (Eq. 12.431 ) only after fully taking into account the non- linear 
effects in the classical EOM 12.81 . By contrast, in the quantum case, the 
non-linear effects responsible for saturation have already been encoded 
in the distribution of the color sources. This is natural since, as noticed 
at the end of Sect. 2, the difference between "sources" and "wee gluons" 
is a matter of convention and depends upon the resolution scale in k~^ . 
As anticipated, the gluon occupation factor at saturation, Eq. 13.441 . is 
parametrically of order 1/as, as in the MV model. 



s Recall that Eq. I3.2HI corresponds to a fixed coupling Qs. In the case of a running 
coupling, where Eq. 13.301 applies, Eq. 13.441 remains formally the same, but as in 
the denominator must be understood as the running coupling <3.29i evaluated at = 
k±Qs. 



1, 2008 20:44 



WSPC/Trim Size: 9in x 6in for Review Volume 



QGPS'finar 



The Color Glass Condensate and High Energy Scattering in QCD 55 

To conclude, the unintegrated gluon distribution predicted by the 
CGC effective theory is illustrated in Fig. 1171 and looks qualitatively 
similar to that in the MV model, cf. Fig. El At large kj_ ^ Qs{t), the 
distribution in Fig. ll7l is given by the solution to the BFKL equation, cf. 
Eqs. llOTl - llI!^ . while at low k± < Qs(r), it is given by Eq. llOIl . 
The saturation condition 113.21^ can be also formulated in terms of the 
unintegrated gluon distribution as 

^r{kj_) for fc^~Q4r). (3.45) 

as 

This condition, together with BFKL evolution at higher momenta, im- 
plies that, for kj_ within the scaling window 13.2811 . the unintegrated 
gluon distribution has a scaling form similar to Eq. 13.2611 

..(.,). ^(^)\ (3.46) 

where k is a yet undetermined numerical prefactor. We see that the 
scaling property characteristic of saturation (cf. Eq. 13.4411 ') is preserved 
by the linear BFKL evolution up to a relatively large moemtum kj_ ~ 
Qlir) / Aqqi) (cf. Eq. 13.281 ). which is well above the saturation scale. 

This "extended scaling" region where the gluon density is relatively 
low, but takes the scaling form 13.46L is represented on the diagram in 

Fig. im 

iii) Color neutrality at saturation 

The vanishing of the charge-charge correlator (see 13.431 '). as fc^ 
when ^ has important consequences for the infrared behaviour of 
the CGC effective theory. Consider, the calculation of the dipole-hadron 
S-matrix element. In the MV model, it was found to be logarithmically 
infrared divergent, cf. Eqs. or JTHl . If in Eq. fTTll one replaces 

the MV estimate for the charge correlator A^j by the corresponding quan- 
tum expression, the ensuing integral over k^ becomes infrared finite, 
due to Eq. 13. 4H . "When computed in the quantum effective theory, the 
dipole scattering amplitude is infrared safe, and therefore insensitive to 
the non-perturbative physics of confinement. A similar property holds 
for other gauge-invariant quantities like the gluon distribution. 

Furthermore, the smooth infrared behaviour in Eq. 13. 4H is responsi- 
ble for fc^-saturation in the gluon distribution. In what follows, we shall 
argue that this behaviour has a simple physical interpretation-gluonic 
color sources are correlated over long distances to ensure that color neu- 
trality is achieved over a transverse area of order 1/Q^{t) '^^ 

Consider the total color charge enclosed within a surface E, as 
given by Eq. 12.41 with AS±^ — > E. This is a random quantity with zero 
average (since (p°(s)) = at any point x), so we shall compute the 
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average of the color charge squared — Q°'Q'^. We have, 

{Q')r = {n!-1) f d^x^ f d^y^,,r(x^,y^). (3.47) 

In the MV model, where the sources are uncorr elated (cf. Eq. JTHJ), 
we have 

{Q')a = (^c - ~ — SQi , (3.48) 

as 

which increases rapidly with A, like A^/'^ . 

After including quantum evolution, the charge correlator acquires a 
non-trivial momentum dependence ^^iki^), and Eq. 1)3.47^ can 

be estimated as (up to a color factor): 

{Q2)r ~ EMr(fci ~ 1/S). (3.49) 

For a relatively small area, or moderately high energies, ^ 
(3s(r), and /^^(fcj^) is given by the BFKL equation 13.3811 . Then 
Eqs. llS371l and llob imply : 

{Q?)r —{t.qI{t))' for l/S»Q?(r), (3.50) 
as \ ) 

which shows incomplete color shielding. With increasing E, the total 
charge squared enclosed within this surface increases, but not as fast 
as the area itself. Thus, the density (Q^) /T, ~ 1/T,^~'^ vanishes in the 
limit E —> oo, which is consistent with (global) charge conservation. 
Furthermore, for fixed E, the total charge {Q^)t increases exponentially 
with r. 

For larger surfaces, of the order of the saturation disk l/(5s(r) or 
larger, one should rather us eq. 13.4311 . which gives: 

(Q^)r ~ — lnfEQ2(r)) for 1/E « Q2(r). (3.51) 

The total charge squared (13. 5H is still non-zero, but unlike Eq. (UMl 
it is only logarithmically increasing with both 1/x and E. This is com- 
plete shielding : When increasing E, the newly included color sources are 
completely screened by the other sources. If the total charge (squared) is 
still increasing with E (albeit only slowly) , it is because the longitudinal 
width r — Ts(E) ~ ln(EQs(r)) of the "saturated" piece of the hadron 
rises logarithmically with E. 

This complete shielding, together with the fact that the total charge 
13.5111 is much smaller than the total charge for a system of uncorre- 
lated color sources with surface density Qj/qs, cf. eq. 13.4811 . enables 
us to speak about color neutrality already at the relatively short scale 
l/Qsij) <^ 1/Aqcd- This interpretation is further confirmed by the fact 
that the color field created by gluon sources at la rge distances 3> 1/Qs (t) 
cannot be distinguished from a dipolar field 
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3.6. A Gaussian effective theory 

Consider the calculation of the dipole-hadron S-matrix element 12.491 
within the CGC effective theory. After expanding the Wilson lines 
in Eq. (12.4911 . one is led to evaluate n-point functions of the type 
{a{zi)a{z2) ■ ■ ■ a(£n))r, where the transverse arguments 2j_ j are ei- 
ther xj_ or . Each such a n-point function receives contributions from 
either hard {kj_ > Qs{t)) or semi-hard'' (Aqcd < k± < Qs{t)) mo- 
menta. The contributions of the modes with fcj^ ^ Qs{t) are easy to 
evaluate. In the saturation regime, the weight function is the Gaussian 
13.4011 ■ and the only non-trivial correlation is the 2-point function 13.411 1. 
For hard momenta kj_ 2> Qsij), the weight function is strictly speaking 
not a Gaussian. In this regime, fields are rather weak [ga ^ 1), so the 
(yet unknown) contributions of the higher-point correlations are less im- 
portant than that of the 2-point function. The latter is explicitly known 
as the solution to the BFKL equation. Finally, there are the contribu- 
tions of the modes ~ Qs{t), which we don't know how to evaluate 
so far. But these contributions should be relatively unimportant as long 
as — yj_ is either very small, or very large, compared to 1/Qs{t). 
To summarize, in order to compute 5't(t'_l) for a dipole which is either 
very large, or very small, it should be enough to use the previous ap- 
proximations for the 2-point function, which are valid far away from the 
saturation scale. 

These considerations can be extended to any quantity which is not 
very sensitive to the momenta in the transition regime around Qs{t). To 
systematically compute such quantities, construct a Gaussian approxi- 
mation to the weight function which encodes the limiting behavior of the 
the 2-point function at high and low momenta and which interpolates 
smoothly between these regimes. Such a Gaussian will be as simple to 
use as the original MV model, but will extend the latter by including 
the BFKL evolution at high momenta, and the physics of saturation and 
color neutrality at low momenta. 

Such a Gaussian weight function has been constructed in Ref. 
and reads: 

where the kernel ^y{xj_,yj_) — the 2-point function of the color sources 
(cf. Eq. 13. 4111 1 — is such that its Fourier transform Ay(fcj^) satisfies the 
BFKL equation at momenta k^_ ^ Qs(t) and reduces to Eq. 13.4111 for 
momenta fcj^ <^ Qs{r). We can match these limiting behaviours contin- 
uously at k±_ = Qs(t) by imposing the following boundary condition on 



'^As explained above, gauge invariant quantities computed in the effective theory are not 
sensitive to soft {k± ^ Aqcd) momenta. 
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the solution to the BFKL equation: 



-Ql{r) ~ A™(fc^=0.(r)). 



(3.53) 



Then, the following function provides a smooth interpolation, that we 
shall use in Eq. H3.52^ : 



, , _ fciA™(fc^) 



(3.54) 



For momenta within the scaling window l|3.28|l . the BFKL solution takes 
the scaling form (with 7 = 0.64) : 



(3.55) 



which allows us to write a simple explicit expression for the kernel 
Ay(A;j_), valid for all momenta kj_ ^ Q'ii'r)/ Aqcd'- 



Ay(fc^) = 0{T-y)^ 



(3.56) 




1 10 
ki/Qs(xo) 

Fig. 18. Energy and momentum dependence of ipT(k±). We have plotted ipr{k±) as a 
function of fcx/Qs(ro) (with to some value of reference) for six values of r. The lines, 
from the bottom to the top, correspond successively to r = tq, tq + 2, ■ • • , tq + 10. The 
increase with r is exponential at high momenta (giving eguidistant curves in this log— log 
plot), but only logarithmic at low momenta. From Ref. Iliiil 
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In Sect. 4.3 , Eq. I3.5(ill will be used to compute the dipole-hadron 
scattering amplitude. Here, we shall present a different application, 
namely we shall deduce a simple analytic expression for the uninte- 
grated gluon distribution. Specifically, Eqs. (13.3611 and 13.411 imply 
~ f dyXy{k±^). When further combined with Eqs. 13.371 and 
H3.56|l . it leads to the following final result: 

Mk.) - ^ / dyXAk.) - ^ m (l + (^y) . (3.57) 

(In performing the integral over y, we assumed Eq. H3.23^ for the satura- 
tion scale.) Eq. 13.571 interpolates smoothly between Eq. 13.441 deeply 
at saturation {kj_ <^ Qs) and Eq. 13.461 for momenta within the scal- 
ing window Qs ^ k^ ^ 0^1 ^QCD- In Fig- H^l we illustrate both the 
A;^ -dependence and the r-dependence of the function 13.571 . which is 
plotted as a function of fcj^ for several values of r. 

4. Deep Inelastic Scattering and the CGC 

In this section, we shall discuss some applications of the Color Glass 
Condensate picture to Deep Inelastic Scattering. We start by explicitly 
deriving a factorization formula introduced in Sect. 2.7, which allows 
one to compute the F2 structure function at high energy in terms of the 
dipole scattering. We shall then discuss a phenomenological "saturation 
model", proposed by Golec-Biernat and Wusthoff, which is based on 
similar physical premises, and compares remarably well with the HERA 
data. This model will be further reconsidered from the CGC perspective, 
in which the dipole-hadron scattering is computed within QCD. This al- 
lows us to relate the "geometric scaling" observed in the HERA data 
to properties of the quantum evolution towards saturation. In addition, 
we discuss the Froissart bound can be realized conceptually in the Color 
Glass Condensate. We discuss next the relation of saturation to shadow- 
ing in DIS. Finally, we discuss inclusive and semi-inclusive signatures of 
the Color Glass Condensate in Deep Inelastic Scattering and specifically 
how they may be observed at a future Electron Ion Collider. 

4.1. Structure Functions in the Color Glass Condensate 

We shall show here how one computes structure functions for DIS in 
the presence of non- linear effects associated with saturation . Towards 
this end, one needs the correlator of the electromagnetic current in the 
background of the strong classical gluon field representing the CGC. In 
DIS, the interaction between the hadron and the virtual photon is en- 
coded in the following tensor expressed in terms of the forward Compton 
scattering amplitude T^v^- 

W'"'{q^,P-q) = 2 Disc T^"'{q'^,P-q) 
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i- d'^xe'i-'' < P\T{Jf'{x)J''{0))\P > , (4.1) 



27r 

where "T" denotes a time-ordered product, = tp^^tl) is the hadron 
electromagnetic current and "Disc" denotes the discontinuity of T^i/ 
along its branch cuts in the variable P- q. Also, < is the transferred 
momentum squared (i.e., is the momentum of the virtual photon, 
and = ~q^) and P^ is the momentum of the target. In the IMF, 
P+ — > oo is the only large component of the momentum. Since 

< T{J^\x)r{y)) > = < T (^(a;)7^V(xWj/)7"^(2/)) > , (4.2) 

the time ordered produced of currents can be expressed, in complete 
generality, as 

<T(J^(x)r(jy))> = 
= Tr{j^GA{x))Mj-'GAiy)) + Tiil''GAix,yh''GA{y,x)) , (4.3) 

where GA{x,y) — —i < tp{x)tp(y) >a is the quark Green's function in 
the external background gauge field Aa of the hadron. 

The first term on the right hand side of Eq. 14.31 is a tadpole con- 
tribution without an imaginary part. It therefore does not contribute to 
W". We are thus left with: 

W^^''(g^P.g) = ±^Im [d'X [ /xe^"'^ (4.4) 



2tt M _ 

(tt(j''Ga{X + x/2,X- x/2)'y'' Ga{X - x/2,X + x/2)^Y 

The approximation one makes here is to replace the full background 
gauge field by the classical background field A'^ classical- other 
words, the Green's function which, in general, is computed in the full 
background field of the nucleus, is now computed in the saddle-point 
approximation where — > A'^ classical- Note that this expression makes 

no reference to the operator product expansion of DIS ■ -. Thus, it is also 
valid at small values of x and moderate Q^, where the operator product 
70 

expansion is not reliable but where the classical approximation is 
sensible. 

Since the current-current correlator is gauge invariant, one can com- 
pute it by using background field propagator in an arbitrary gauge. For 
the reasons explained in Sect. 2.4, it is most convenient to use the co- 

7l 'i2ii'iQi 

variant gauge, in which the propagator reads as follows - ' ' ■ 

GA(x,y) = Gq{x ~y)-i J d'^zGoix - z)-/^ 5{z^)Go{x ~ z) 

{e{x~)e{-y-){v\z^) - 1) - e{-x-)e(y-){v(z^) - i)} , (4.5) 
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where Gq is the free propagator and V and are the Wilson Unes of 
Eq. 12.241 . This is obtained by assuming the color source to be a 5- 
function in a;^, which is appropriate since this source is due to relatively 
fast partons with longitudinal momentum fractions much larger than the 
Bjorken x = /2P ■ q of the collision. 

Inserting the fermion propagator 14.511 in Eq. 14.411 and performing 
the integrations there, one obtains the final result for W", which is 
conventionally expressed (for an unpolarized target and for « My^) 
in terms of two structure functions Fi and defined by^ 



,(p^ _ q'^jP-q) ) (pu _ q-{P-q) \ ^2 
^\ g2 J ql ){p.q) 



(4.6) 



Consider first F2: the corresponding result takes the factorized struc- 
ture anticipated in Eqs. 12. 461 - 12. 4811 . where |vl>(2, r-j^)|^, which signifies 
the probability that the virtual photon splits into a gg-pair can be ex- 
pressed as [^{z, r±_)\^ = \^x{z, r^)\^ + \^ i^{z, r\_)\^, and where each of 
the terms has the explicit form 



l*T(^,ri)i' = ^ E ^/ {(^' + (1 - ^f)Q}Kl{Qfr)+m)Kl{Qfr)] , 
f 

l*L(^,^x)l' = ^E4{40'^'(l-^)'^o(Q/r-)} . (4.7) 



Above, the sum runs over the quark flavours, = z{l — z)Q'^ + m?p 
is the quark mass, and Kq and Ki are modified Bessel functions. 
I^t(2)^±)P '"i)!^) denotes the probability that a transversely 

(longitudinally) polarized photon splits into a qq-palv. This decomposi- 
tion implies a similar decomposition for F2, namely, F2 — Fx + -Fl- Fij 
this is proportional to Fy-explicitly, Fx = 2xFi and the longitudinal 
structure function is Fx = F2 — 2xFi . In the parton model, Fx = 0-the 
Callan-Gross relation-but is non-zero in QCD and is directly propor- 
tional to the gluon distribution. An independent measurement of Fx is 
therefore of great phenomenological interest as will be discussed in Sect. 
4.6. One can similarly derive expr essions for the diffractive structure 
functions as discussed in Ref. f-^^. Some of the phenomenological im- 
plications of the CGC picture for measurements of structure functions 
were discussed in Sects. 2.6 and 2.7. 



4.2. The Golec-Biemat-Wusthoff model 

In Ref.l^SI Golec-Biernat and Wiisthoff introduced a simple phenomeno- 
logical model for the dipole-hadron cross-section, Eq. I|2.48|l . which is 
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generally referred to as the "Saturation Model" : 



'^dipole(x, rj.) = ao (l - e-'--'5=('')/4) (4.8) 



with the parametrization Qs(x) = Qq{xq/x) . This shows color trans- 
parency at low rj^ <^ 1/Qs(x) and respects the unitarity bound, as it 
approaches a constant value (Tq when r_L ^ l/(5s(x). If one interpretes 
(Tdipoio (x, ) /ao as a scattering amplitude (compare to Eq. I|2.48^ ) , then 
this has no dependence upon the impact parameter bj_. Rather, this 
should be seen as an average of the true scattering amplitude over all 
impact parameters. This amplitude is consistent with the physics of sat- 
uration (compare to Fig. IIS! , so it is natural to interprete the scale 
Qs(x) in Eq. 14.81 as a "saturation momentum". Note however that the 
approach towards the unitarity bound J^T{r±) = 1 for large rj_ is much 
faster then predicted by the CGC approach, or the solution to the BK 
equation (see Sect. 4.3 below). Besides, at small rj_, Eq. (14.811 fails to 
reproduce, as it should, the leading-twist approximation-it misses the 
logarithmic factor ln(l/r^A^) in the exponent of Eqs. 12.5111 . Neither 
does it recover the BFKL prediction 13.61 . 

By using the dipole cross-section 14.81 and the factorization formula 
iTTfll for DIS, Golec-Biernat and Wiisthoff were able to fit the HERA 
data for a^^p for x < 10~^ and a wide range in with only three 
parameters, namely ctq ~ 23 millibarns, A = 0.288 and xg = 3.04 ■ 10~^. 
(The reference scale Qq has been fixed as Qq = 1 GeV^.) These fits 
were performed for three light quark fiavors; with the addition of charm 
quarks, the best fit was obtained with slightly changed parameters: (tq = 
29.1mb, A = 0.277, and xq = 0.4 • lO"'*. 

The fits are reasonable for up to ~ 20 GeV^ but are less suc- 
cessful beyond. This is related to the above observation that Eq. 14.81 
does not have the right perturbative behaviour at small r_L . By replacing 
Eq. 14.81 with a Glauber-type formula as shown in Eq. 12.531 . the situ- 
ation at high improves considerably an d a wider range in can be 
fit within the framework of this model For other phenomenological 
models of the HERA data based on ideas of saturation see Refs. ' ' '. 

The diffractive structure function F^{x-p,Q^,P) (where x-p = 
{M'^ + Q'^)/{W'^ + Q'^) and /3 = /{M'^ + 0^) where M is the diffractive 
mass and W is the total energy of the virtual photon-proton process) 
was measured at HERA and several striking properties of the diffractive 
structure function were observed. For instance, the ratio of a'^^^^/cj is 
large and is nearly independent oi W , a feature that was not anticipated 
in pQCD based models. In Ref. ^31 the simple model that describes the 
inclusive scattering data also describes the diffractive structure function 
data. Interestingly, the form of the diffractive cross-section in this model 
is similar to the inclusive one except the dipole-hadron cross-section ap- 



February 1, 2008 20:44 



WSPC/Trim Size: 9in x 6in for Review Volume 



QGPS'finar 



The Color Glass Condensate and High Energy Scattering in QCD 63 

pearing in the latter is replaced by the square of this cross-section The 
parameters appearing in the fit of the inclusive cross-section to the data 
are therefore the same as those used in the diffractive fit. The agree- 
ment with data is quite impressive given these constraints. A version of 
the model has als o been a pplied to study vector meson production with 
reasonable results ^^^l^^l 

Whereas the behaviour at high can be improved rather easily, 
by replacing Eq. 14.81 by Eq. 12.5311 as mentioned above, the impact 
parameter dependence is problematic. This is expected to be especially 
important in a proton. Data on the f-dependence of differential cross- 
sections, particularly for vector meson photoproduction, may be useful 
in understanding parton distributions in impact parameter space ESilZI 



4.3. Geometric Scaling in DIS 

The Golec-Biernat-Wiisthoff dipole cross-section, Eq. 14.81 1. has the 
remarkable feature to depends upon the two kinematical variables x 
and r±_ only via the dimensionless combination (the "scaling variable") 
T = rJ^Qs(x). Via the factorization formula 12.471 . this scaling prop- 
erty transmits to the virtual photon total cross-section a~^*p which, in 
the limit where the quark masses are negligible, becomes a function of 
the ratio (3^/(3s(x) alone (a property usually referred to as "geomet- 
ric scaling"). At a first sight, this may appear an artifact of the simple 
parametrization 14.81 specific to the saturation model. Indeed, while the 
scaling looks natural at saturation one would naively expect this scaling 
to be broken after generalizing Eq. 14.81 to reproduce the perturbative 
behaviour at high (e.g., the inclusion of the logarithm ln(l/r^A^) 
as in Eq. 12.511 would clearly violate scaling). 

It thus appeared as a surprise when Stasto, Golec-Biernat and 
Kwiecihski showed that, to a rather good accuracy, the HERA data 
on cr-y*p do show scaling for small enough x (x < 0.01) and all up 
to 450GeV2 (see Fig.[T^. Such are significantly higher than the es- 
timated value of the saturat ion scale at HERA, as extracted from the 
"saturation model" fits to ^2'^ : Qs - 1 • • • 2 GeV^. On the other hand, 
the data show no scaling for larger values of x. 

Subsequently, some indicatio ns o f geometric scaling have been found 
also in t he d ata for DIS off nucleil^, and even in the particle production 
at RHIC^^, although, in these cases, the experimental evidence is more 
uncertain. 

In the light of the previous discussion in Sect. 3, such a scaling be- 
haviour should not look surprising any more. As explained in Sects. 
3.4 (for the scattering amplitude) and 3.5 (for the gluon distribution). 



'This feature of the model can be understood very simply in the CGC pictureESni 
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Fig. 19. HERA data on the cross section for 7*p DIS from the region x < .01 and 
< 400GeV^ plotted versus the scaling variable T = Q^/Qg(x) (from Ref.E^. 



the scaling property at saturation is preserved by the BFKL evolu- 
tion up to relatively large Q^, of order Qt{''')/^QCT> i'^^- -^^l- 
If one uses the phenomenological values of Qs alluded to before, and 
Aqcd = 0.2 ■ • ■ O.SGeV, one finds that the maximum up to which 
scaling is expected is indeed of the order of a few hundred GeV^, as 
observed at HERaI^. 

Within the CGC formalism, the scaling properties of the dipole- 
hadron scattering amplitude can be explicitly studied by using the Gaus- 
sian approximation to the weight function introduced in Sect. 3.6. The 
kernel of the Gaussian shows explicit scaling for all momenta kj_ ^ 
Q^{t)/ Aqqd, cf. Eq. Il3.5till . and therefore so does also the S'-matrix 
element compute d in this approximation. Specifically, a straightforward 
calculation yields (in the case of a fixed coupling, cf. Eq. I|8.23^ , for 
definiteness, and with n = AttCp/'jcNc) : 




exp 



(2^): 



kl 



In 



Unlike the corresponding prediction of the MV model, Eq. H2.51|l . which 
was infrared sensitive, and thus dependent upon the non-perturbative 
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scale Aqcd, the integral in the equation above is well behaved both 
in the infrared and in the ultraviolet, so, clearly, the result is a scaling 
function: Srir^) = S{rj_Qs{T)). 

The limiting behaviour of Eq. (14.911 at either small, or lar ge, trans- 
verse size r_L as compared to the saturation length 1 /Qs (t) arel^^ i) For 
a small dipole, r_|_ <C 1/Qs{r), the dominant contribution to the integral 
over k±_ comes from hard momenta k±_ ^ Qs{t), from the scattering off 
color sources located at \zj_ — b^ \ <^ 1/Qs(r) from the impact parameter 

of the incoming dipole. 

If 7 < 1 (in particular, for the BFKL value 7 ~ 0.64), this domi- 
nant behaviour has the power law form in Eq. 13.26L with "anomalous 
dimension" equal to 7. 

If 7 = 1 (the DGLAP value EHI, which supplants BFKL evolution at 
very small r^): 



ln^^+(l + 2^(2) + 21n2) 



(4.10) 



This is similar to the small-rj^ behaviour in the MV model, Eq. 12.521 . 
except that the infrared cutoff in the logarithm is now the saturation 
momentum, and not Aqcd- 

i) For a large dipole, r± » 1/Qs(t) (but r± < 1/Aqcd), the 
dominant contribution comes rather from momenta fc_L in the range 
l/rj_ <^ kj^ <^ Qs(t), that is, from scattering off saturated gluon 
sources, and reads : 

S.(r^)(xexp|-^(lnriQ?(r))'| , (4.11) 

in agreement with the results in Refs. ES^E^Hl ^s anticipated, the ap- 
proach towards the "black disk" limit Sr = predicted by the CGC is 
slower than that assumed in the saturation model, Eq. 14.81 . 



4.4. The Froissart Bound for dipole scattering 

Let us now address the fundamental question of the asymptotic be- 
haviour of the total cross-section at very high energies. Thus far, in the 
analysis of the dipole-hadron scattering, we have neglected the impact 
parameter dependence. We have assumed the hadron to be a homoge- 
neous disk with radius R. Then all the theoretical descriptions which 
include saturation, from the phenomenological "saturation model' in 
Eq. 14.81 to the QCD-based formalisms like the CGC or the Balitsky- 
Kovchegov equation, lead a dipolar cross-section which approaches a 

constant value aa = as s 00. Although consistent with the 

4 5 

Froissart bound—"*", this result cannot be right. From experiments, we 
know indeed that hadronic cross-sections (e.g., for pp collisions) keep 
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growing with s up to the highest energies that have been reached so 
far. At very high energy, the growth is rather slow, and can be fitted by 
either a small power of s (the "soft pomeron" ~ s"'^*), or some power 

81 

of In s, or a combination of them. A recent analysis - - of the data for 
several high-energy processes appears to favour a dominant behaviour of 
the log-squared type: (Jtot(s) — o"o In^ s, with an universal prefactor ctq. 
If true, this would imply that the Froissart bound is actually saturated 
in nature. 

Such a steady growth of total cross-sections is also expected on phys- 
7R 7Q 80 

ical grounds . As a quantum mechanical bound state, the hadron 

does not have a sharp edge but rather a tail. In a theory with a mass 
gap like QCD, this tail is an exponential whose width is fixed by the 
lowest mass in QCD, namely the pion mass. Impact parameters which, 
at some initial energy, are far away in the tail of the distribution will 
not contribute much to scattering, since the parton density is low there. 
But with increasing energy, the local gluon density will increase rather 
fast as a power of s, since gluon recombination is not effective when the 
density is low. Eventually, for sufficiently large s, the local gluon den- 
sity will become high enough for these impact parameters to contribute 
significantly to scattering. That is, with increasing energy, the effective 
interaction radius of the hadron is expected to grow as well, which then 
results in an increase of the total cross-section. In fact, since the local 
scattering amplitude A/t-(?'_l, bj_) cannot exceed the unitarity, or "black", 
limit A/t = 1, it is clear that, for sufficiently large energies, the increase 
of the cross-section with s will proceed via the expansion of the "black 
disk" (= the central area of the hadron where the unitarity limit has 
been reached already). 

This general discussion shows that a theoretical description of the 
dynamics of the black disk must combine two essential ingredients: i) a 
mechanism which ensures the unitarization of the scattering amplitude 
at fixed impact parameter, and ii) a description of the tail of the hadron 
wavefunction. In QCD, the second issue is certainly related to confine- 
ment, and is thus genuinely non-perturbative. But it has been unclear 
until recently whether the first issue, that of the unitarization, can be ad- 
dressed in perturbation theory or not. Indeed, since gauge interactions 
are a priori long-ranged, it could well be that soft, non-perturbative, 
interactions are responsible for the approach towards "blackness" at a 
fixed impact parameter. The "infrared diffusion" of the BFKL equation 
may be seen as in argument in that sense. 

However, our present understanding of the physics of high parton 
densities shows that the phenomenon of gluon saturation provides a 
perturbative framework for the study of unitarization. We have seen 
indeed that the perturbative, but non-linear, evolution equations yield a 
scattering amplitude which respects the unitarity bound (see Fig. llSl and 
Eq. lE^J). To study the expansion of the black disk, these equations must 
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r(q2) Rh b 



Fig. 20. Dipole-hadron scattering in longitudinal (left) and transverse right) projec- 
tions. 



be supplemented with some information about confinement and applied 
at impact parameters in the tail of hadron, namely, in the "grey area" 
at > R{t,Q'^). 

Here, R{t, Q^) is the radius of the black disk for an incoming dipole 
with transverse resolution = l/r\ and relative rapidity r = Ins/Q^. 
That is, Mrir ^,bi^ is of order one for < R{t,Q^), but it drops 
rapidly for b± > R(t,Q^). Equivalently, R(t,Q^) is such that the local 
saturation scale Qs{T,bj_) — which is largest towards the center, and 
decreases with b±, so like the gluon density — becomes equal to 
at b± ~ R{t,Q^) (see Fig. 1201 . Perturbation theory is valid as long as 
Qs{T,b±) > Aqcd, that is, for b± < Rh{t) in the plot in Fig. I^Hl This 
condition leaves us with a perturbative grey corona at R{t,Q^) < b± < 
Rh{''')^ which is sufficient, as we shall see, to perform a controllable 
calculation of the expansion rate of the black disk. 

Specifically, for impact parameters within this corona, the dipole 
undergoes mostly hard scattering, with transferred momenta fc^ ^ 
Q's('^!b±) (cf. the discussion after Eq. 14.91 '). This implies that it in- 
teracts predominantly with color sources which are relatively close in 
impact parameter space, within a saturation disk of radius l/QsiT,bj_) 
centered at b^. Physically, it is so because of the screening phenomena 
associated with saturation (cf. Sect. 3.5) : Color sources which lie further 
away create only dipolar fields at the impact parameter of the incoming 
dipole; such fields decrease rapidly with distance (more rapidly then the 
monopole fields due to the nearby color sources), and thus contribute 
less to scattering I^SEH. 
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Mathematically, this means that scattering in the grey corona is 
controlled by the linear BFKL equation 13.41 , but with an infrared cutoff 
of order Qsij, b±). This cutoff simulates the non-linear terms, which, in 
the full equation 113. 7II . would limit the range of the dominant interactions 
to 1/Qs(t, Clearly, with such a cutoff, the BFKL equation is not 
afflicted by "infrared diffusion" any more. As in the previous calculation 
of the saturation scale in Sect. 3.4, the black disk radius can be computed 
by first solving the BFKL equation, and then imposing a saturation 
condition similar to Eq. (13. 211 : 

A/'r(r_L, fex) ~ 1 for h ~ R{t,Q^). (4.12) 

The -dependence of the initial condition is determined by the 
non-perturbative physics of the confinement, and thus requires a model. 
However, what we need to know about confinement is quite limited, 
and can be inferred from general principles. Firstly, what is the typical 
scale for transverse inhomogeneity in the hadron: this is clearly I/Aqcq. 
Secondly, what is the 6 j_ "dependence of the scattering amplitude in the 
hadron tail. This is an exponential fall-off oc e^^™"'' where twice the 
pion mass enters since the long range interactions between the dipole 
and the (isosinglet) gluons require the exchange of at least one pair of 
pions (recall that pious have isospin one). 

The first requirement tells us that the inhomogeneity occurs over 
transverse scales much larger than the typical range of the interactions, 
^ l/(3s(i", &_!_). The scattering of the small dipole proceeds quasilocally 
in the impact parameter space^. As a consequence, the fo_L -dependence 
of the solution to the ( effec tive) BFKL equation factorizes, and is fixed 
by the initial condition 1^ : 

J^r{r±,b^) ^ J^r{r^)S{b^), (4.13) 

where Mt(t\^ is the solution to the homogeneous BFKL equation, as 
shown in Eq. 13.611 . The second condition provides the transverse pro- 
file function in the hadron tail: S{b±) oc g"^™'^^ (the proportionality 
coefficient is not important to the present accuracy). 

The previous arguments, together with Eq. 13.611 lead to the following 
estimate for the scattering amplitude valid at impact parameters in the 
grey corona'*" (with = l/r'\_) : 

A/V(Q^ &±)|^^^^ ^ exp|-2m^& + a;a5r- iln^l . (4.14) 

This equation, together with the blackness condition 14.1211 . imply: 

i.(.,Q^). ^ ilnf). (4.15) 

■I This would not be true in the absence of the infrared cutoff generated by the non-linear 
effects. 

'^At very high energies, we can also neglect the diffusion term in Eq. 13.61 . 
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We see that the black disk radius increases linearly with r. This be- 
haviour ensures the compensation between the exponential increase with 
T (due to the perturbative BFKL evolution) and the exponential decrease 
with (due to confinement). 

Since A/V is rapidly decreasing at b± ^ R(t,Q^), the total cross- 
section is dominated by the black disk: 

adipolc(T,Q') ^ 27vR\t,Q^). (4.16) 

Together with Eq. 14.151 . this yields the following dominant' behaviour 
of the cross-section at high energieJ^^ (see also Ref. 1^ for an early 
study): 

o-dipolc(s,0^) ~ ^ f ) In^s as s ^ oo. (4.17) 



2 y niTT 

This saturates the Froissart bound, with a universal coefficient for 
all hadrons and reflects the combined role of perturbative and non- 
perturbative physics in controlling the asymptotic behaviour at high 
energy. This behaviour is in qualitative agreement with the phenomeno- 
logical analysis in Ref. but the coefficient in Eq. 14.171 is too large to 
fit the data. The difference with respect to the data can be substantially 
reduced by using the RG-improved NLO estimate for the BFKL inter- 
cept, which decrease s th e leading-order value ujas by roughly a factor of 
^jjj.gg|9bl97| j^gfg_IH21^ tjjg result in Eq. 14.1711 has been generalized to 
7* — 7* scattering. 

From Eq. 14.141 . one can also estimate the radius Rh(t) where the 
saturation scale decreases to Aqcd (beyond which perturbation theory 
fails to apply) . This is the same as the black disk radius for a large dipole 
with ~ ^QCD- Therefore 

RHir)^l^r. (4.18) 

and the radial extent of the perturbative corona can be estimated as: 

Rh{t) - i?(r, q2) ~ -1- In |! , (4.19) 

This is independent of r, and much larger than l/mjr (because of the 
large logarithm ln{Q'^ /A^)), which demonstrates the self-consistency of 
the previous calculation: When t ^ t + dr with dsdr ~ 1 (the typical 
rapidity increment at high energy), the black disk expands from R{t, Q^) 
to R(r, Q^) + uj/2mTr <^ Rnij), cf. eq. 14.151 . and thus remains in the 
region controlled by perturbation theory. That is, the expansion of the 



'The subdominant behaviour is not corrmletely described by Eq. I4.15i . since it also 
receives contributions from the grey area''^. 
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black disk proceeds within the perturbative "grey" corona for rapidity 
intervals which are large enough to allow for the calculation of the rate 
of this expansion 

The perturbative evolution equation Eq. (13.711 should not be used for 
very large impact parameters bj_ ^ Rh{t), where Q^(r, 6) ^ ^qcd- 
Such equations lack confinement, so the long-range dipolar fields gen- 
erated by the gluons within the black disk can propagate to arbitrarily 
large distances, thereby creating power-law contributions to the hadron 
tail. At impact parameters within the grey corona, these long-range 
contributions are relatively small, and the evolution is driven by the 
short-range interactions, with the conclusions outlined above. But for 
sufficiently large impact parameters, well beyond Rh{t), the power- law 
contributions will eventually supersede the exponentially decreasing one, 
and the hadron will develop an unphysical power-law tail. (This is un- 
physical since, in the real world, it is removed by confinement.) If one 
pushes the perturbative expansion until the black disk enters this power- 
law ta il, t hen its expansion rate speeds up, and violates the Froissart 
boundl^. This is clearly an artifact of using perturbation theory outside 
its range of validity. The only way to circumvent this difficulty without 
introducing ad-hoc modifications in the evolution equations (to account 
for confinement) is to follow the perturbative evolution for only a limited 
interval of "time" At, such that perturbation theory remains valid. As 
explained above, a limited evolution in time is indeed sufficient for the 
calculation of the expansion rate of the black disk. The result of this 
calculation can be then extrapolated to arbitrarily high energies, even 
though the perturbative evolution eventually becomes meaningless. 

To conclude, while perturbation theory alone appears to be sufficient 
to describe unitarization at fixed impact parameter, one still needs some 
information about the finite range of the strong interactions in order 
to compute to tal cross-sections. This is reminiscent of an old argument 
by Heisenbergl^ which combines unitarity and short-range interactions 
(as modelled by a Yukawa potential) to derive cross-sections which satu- 
rate the Froissart bound. Fifty years later, our progress in understanding 
strong interactions allows us to confirm Heisenberg's intuition, and iden- 
tify short-range interactions with confinement, and unitarization with 
saturation. 

4.5. Saturation and Shadowing in Deep Inelastic Scattering 

Shadowing is the phenomenon where F-^{x,Q^)/AF-^ {x,Q^) < 1 at 
small X {x < 0.1). It is a large effect in the region where the coherence 
length of the probe (the qq-pair in DIS) l^oh ~ 2rnNX exceeds the intra- 
nuclear longitudinal distance between any two nucleons in the nucleus. 
The nuclear parton distribution is not merely the sum of nucleon parton 
distributions but also contains the intereference between the parton dis- 
tributions of the nucleons. When the coherence length is larger than the 
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nuclear diameter {Icoh ^ 2A^^^, or x << l/{Am]\fA^^^)), the qq-pair 
interacts coherently with the entire nucleus, and the collective effects are 
expected to be imporatant. 

There are several questions about shadowing that have not been 
unambiguously resolved in the framework of QCD. 

• Is shadowing a "leadig twist" effect, or it is suppressed by powers 
of Q^? An empirical answer to this question would he lp s ettle whether 
shadowing is an intrinsically leading twist phenomenon''^ or whether it 
is due to weak coupling, higher twist/high parton density effects^^^^5 

• What is the relation of shadowing to parton saturation? Does parton 
saturation provide a microscopic understanding of shadowing? 

• Does the shadowing ratio "saturate" at a minimum value for fixed 
and A with decreasing x? Does it saturate faster for quarks or gluons? 

• What is the relation of shadowing in nuclei to diffractive scatter ing 
of nucleons? The relation is well established at low parton densities'^ 
In an interesting recent exercise, it has been shown that diffractive nu- 
cleon data at HERA could be used to predict the shadowing of quark 
distributions observed by NMC . Significant deviations from the 
simple relation between shadowing and difi'raction, may again suggest 
the presence of strong non-linearities. 

• Is shadowing universal? For instance, are gluon structure functions 
extracted from p-A collisions at RHIC identical to those extracted from 
e-A in the same kinematic region? The naive assumption that this is 
true may be incorrect if higher twist effects are important. 

The answers to these questions will only be conclusively settled by 
the next generation of p-A collider experiments at RHIC and LHC '^ 
and by e-A collider experiments at RHIC (EIC/eRHIC) and/or DESY 
(HERA III)-see the following sub-section for a discussion of e-A collider 
plans. 

From the theoretical perspective, there are several shadowing models 

1 QQ QQ QQ 

which consider only leading twist shadowing' . Namely, shadow- 

ing effects are put in the non-perturbative initial conditions, which are 
then evolved in and x using the leading twist DGLAP equations ^1 
Modifications of these leading twist models to incl ude Mueller-Qiu type 
non- linear contributions have also been considered 

In the McLerran-Venugopalan saturation model, the gluon distribu- 
tion (defined as the integral over the unintegrated gluon distribution) is 
additive in A. However, the structure functions computed as discussed 
in section 4.1 will exhibit shadowing for < Q^. There have been a 
few attempts to compute shadowing from the non-linear renormaliza- 
tion group equation. One of these calculations, which is based on the 
non- linear equation derived in Ref. ^5 predicts that perturbative gluon 
shadowing will become large as one goes to smaller x^j's^^. Other more 
recent computations of shadowing have been performed within the con- 
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text of numerical solutions of the Balitsky-Kpvchegov equation for both 
inclusive and diffractive scattering"^. In both cases, geometric 

scaling of the nuclear distributions is claimed but it is not clear that 
both groups obtain the same ^-dependence for the saturation scale. A 
very interesting recent theoretical suggestion is that the scattering am- 
plitude for high energy scattering for nuclei, at fixe d im pact parameter, 
is the same as for protons at asymptotic energies! Integrated over 
impact parameter, this would suggest that shadowing saturates at very 
small X. The understanding of shadowing in the saturation picture is 
still preliminary-more detailed global fits of the non-linear equations (as 
for instance performed for the linear DGLAP fits) are needed. In addi- 
tion, different computations of the A-dependence of the saturation scale 
will likely converge as our theoretical understanding improves. Finally, 
as will be discussed in the following, the issue will likely not be resolved 
conclusively until DIS experiments off nuclei at small x are performed. 

4.6. Probing the CGC with an Electron Ion Collider 

A high energy electron-nucleus collider, with a center of mass energy 
^/s = 60-100 GeV, presents a remarkable opportunity to explore fun- 
damental and universal aspects of QCD. The nucleus, at these energies, 
acts as an amplifier of the novel physics of high parton densities-aspects 
of the theory that would otherwise only be explored in an electron- 
proton collider with energies at least an order of magnitude greater than 
that of HERA. An electon-nucleus collider will also make the study of 
QCD in a nuclear environment, to an extent far beyond that achieved 
previously, a quantitative science. In particular, it will help complement, 
clarify, and reinforce physics learnt at high energy nucleus-nucleus and 
proton-nucleus collisions at RHIC and LHC over the next decade. For 
both of these reasons, an eA collider facility represents an important 
future direction in high energy nuclear physics. 

We will briefiy discuss here experimental observables in deep inelastic 
scattering (DIS) which are signatures of the novel physics of the Color 
Glass Conde nsate ™. A more detailed discussion of the following can be 
found in Ref . EHU. 

The regime of small x^j's {xgj < 0.01) is easily accessed by an 
electron-heavy ion collider in the energy range ^/s « 60-100 GeV. These 
energies would be most natural if the Electron Ion Collider (EIC) were 
constructed at BNL-this particular realization is called eRHIC. The kine- 
matic coverage of EIC/eRHIC is shown in Fig. 1211 What is novel about 
these energies is that for the first time one can study the physics of 
XBj « 0.01 in a nucleus for » ^^cD^ where Kqcd ~ 200 MeV. 



™We will not cover here the interesting physics at intermediate and large xsj tha t can be 
studied with an eA collider. A nice discussion of these issues can be found in Ref.l^^^^^. 
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Previous (fixed target) experiments such as NMC and E665 and current 
ones such as HERMES and COMPASS could only access small x^j at 
small Q^'s. The center of mass (cm) energy of the Electron-Ion Col- 
lider (EIC) is a factor of 10 smaller than that of the current ep-collider 
at HERA (the proposed HERA HI plan-which includes e-A scattering 
would have a center-of-mass energy that's roughly 3 times greater than 
EIC). However, an eA collider has a tremendous advantage-the parton 
density in a nucleus, as experienced by a probe at a fixed energy, is much 
higher than what it would experience in a proton at the same energy. 
Since the parton density grows as A^^^, this effect is more pronounced 
for the largest nuclei-see Eg. 11.131 




Fig. 21. The x— range of the electron ion collider (EIC) compared to that of the 
HERA ep collider and fixed target experiments. The EIC's reach would encompass the 
fixed target regime as well as part of the HERA regime. 

Eq. 11.131 suggests that a^pioton ~ a^nuclcus/ ) • Since the nu- 
cleus is dilute and conservatively taking the effective A^^^ — 4, then, 
for A ~ 0.3, one finds iproton ~ a^nuclcus/lOO. Thus the same parton 
density in a nucleus at x^j ~ 10"'* and ~ a few GeV^ is attained 
in a nucleon at xgj ~ 10^^ and similar Impact parameter tagging 
is feasible by counting knock-out neutrons ^^Ljf go, the gain in parton 
density in eA relative to ep may be even more spectacular. 

In the following, we will discuss both inclusive and semi-inclusive sig- 
natures of the CGC. The latter in particular are very difficult to measure 
in fixed target DIS 
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Inclusive signatures of the CGC 

An obvious inclusive observable is the structure function F2{xgj , Q^) 
and its logarithmic derivatives with respect to xsj and Q^. The EIC 
should have sufficient statistical precision for one to extract the loga- 
rithmic derivative of F2 (and its logarithmic derivative!). Whether the 
systematic errors at small xgj will affect the results is not clear at the 
moment. The logarithmic derivative dF2/d\n{Q^), at fixed xgj, and 
large Q^, as a function of Q^, is the gluon distribution. QCD fits imple- 
menting the DGLAP evolution equations should describe its behavior at 
large . At smaller , one should see a significant deviation from lin- 
ear QCD fits-in principle, if the range is wide enough, one should see 
a turnover in the distribution. The at which the turnover takes plac e 
should be systematically larger for smaller x's and for larger nuclei 

At eRHIC one can extract the longitudinal structure function 
Fi^{xgj, Q^) = F2 — 2xBjFi at small xgj independently since the energy 
of the colliding beams can be varied significantly. In the parton model, 
F]^ = 0-thus F]^ is very sensitive to scaling violat ions. It provides an 
independent measure of the gluon distribution and in particular of 
higher twist saturation effects which may be prominent in both Fj^ and 
Ft but may cancel in the sum^^. 

The extended kinematic range of EIC may help determine whether 
shadowing is entirely a leading twist phenomenon, or if there are large 
higher twist perturbative corrections. As also discussed previously, there 
is a close relation between shadowing and difi^raction. At EIC the validity 
of this relation can be explored directly-different nuclear targets are 
available, and the diffractive structure function may also be measured 
independently. 

Semi-inclusive signatures of the CGC 

A striking semi-inclusive measurement is hard diffraction wherein 
the virtual photon emitted by the electron fragments into a final state 
X, with an invariant mass Mx 3> ^qcDj while the proton emerges 
unscathed in the interaction. A large rapidity gap-a region in rapidity 
essentially devoid of particles-is produced between the fragmentation 
region of the electron and that of the proton. In pQCD, the probability of 
a gap is exponentially suppressed as a function of the gap size. At HERA 
though, gaps of several units in rapidity are unsuppressed; one finds that 
roughly 10% of the cross-section corresponds to hard diffractive events 
with invariant masses Mx > 3 GeV. 

Hard diffraction probes the color singlet object (the "Pomeron") 
within the proton that interacts with the virtual photon. It addresses, 
in a novel fashion, the nature of confining interactions within hadrons. 
The mass of the final state is large and one can reasonably ask ques- 
tions about the quark and gluon content of the Pomeron. A diffractive 
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structure function F^^^ can be defined ^^^^^ (analogous to F2) as 



dxBjdQ^dx-pdt xQ^ [ 2[1 + R^^^\l3,Q^,x-p,t)]j 

F^j^\l3,Q\xT,,t), (4.20) 

where y = Q^/sxbj and ii^^'^^ = F^'^'^^ / F^^'^\ Also, = -g^ and 
^Bj = /(^P ■ q) have the usual DIS definitions and 

xr = . , = (p _ p'f ,f,= ^. (4.21) 

Here P is the initial nuclear momentum, P' the net momentum of the 
fragments Y in the proton fragmentation region and Mx the net mo- 
mentum of the fragments X in the electron fragmentation region. An 
illustration of the hard diffractive event is shown in Fig. 22. Unlike F2 
however, F^^"^^ is not truly universal-it cannot be applied, for instance, 
to predict diffractive cross-sections in p-A sca tterin g; it can be applied 
only in other lepton-nucleus scattering studies EiiSI 




In practice the structure function F^^^ ~ J F^^^ dt is measured, 
where ltmi„l < \t\ < \tmax\, where \tmin\ and \tmax\ and the limits of 
the empirically measurable momentum transfer to the nucleus. The ratio 
for two nuclei Al and A2 



RAl.A2{f3,Q\xT,) = ^^ , (4.22) 
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can be measured with high accuracy' — '. If Rai,A2 ~ !> the structure 
of the Pomeron is universal, and one has an j4-independent Pomeron 
flux. If Rai,A2 = /(^li^2), then albeit a universal Pomeron structure, 
the flux is A-dependent. Finally, if Pomeron structure is A-dependent, 
some models argue that Rai,A2 = ^2,Al/^2,A2- 

The ratio of Rj^) — Cdiffiactivc/'^total a-t HERA is ~ 10% for Mx > 3 
GeV. The systematic s of hard diffraction at HERA can be understood 
in saturation models ESI. For eA collisions at EIC energies, saturation 
models predict that the ratio R^ can be much higher-on the order of 
30% for the largest nuclei 108 

An important semi-inclusive observable in eA DIS at high energies 
is coherent (or diffractive) and inclusive vector meson production. For 
instance, the forw ard v ector meson diffractive leptoproduction cross- 
section off nuclei is ^3 



^\t=o{l*A^VA)^al{Q^) [Ga{x,Q^\ (4.23) 



for large and is therefore a sensitive probe of gluon satura- 
tion/shadowing. It is very important to measure inclusive and diffractive 
open charm and jets since they provide useful and complementary data 
to those of vector mesons. 

Due to large color fluctuations at small xgj, on can expect: a) a 
broader rapidity distribution in larger nuclei relative to lighter nuclei 
and protons, b) Enhanced anomalous multiplicity distributions where 
anomalous multiplicity in one rapidity interval in an event would be 
accompanie d by an anomalous multiplicity in rapidity intervals several 

units away ' and c) a correlation between the central mu ltipli city 

with the multiplicity of neutrons in a forward neutron detector' — -4 

In pA scattering at RHIC one also has the opportunity to study the 
high parton densitiy phenomena-see section 5.5 for a discussion. Some of 
the differences between pA and eA are as follows. In the pA Drell-Yan 
process, it is very hard to reliably extract distributions in the region 
below the tail- namely, one requires > 16 GeV^. In the x region 
of interest, one expects saturation effects to be important at lower 
of 1-10 GeV^ For = 16 GeV^, one might have to go to significantly 
smaller x's to see large saturation effects. Secondly, the survival prob- 
ability of large rapidity gaps is smaller in pA relative to eA. This is 
because in pA (unlike eA) the gap is destroyed due to secondary inter- 
actions between "spectator" partons in the proton and the "Pomeron" 
from the nucleus. Thus one expects that diffractive vector meson and jet 
production in pA should be qualitatively different than what one will 
see in eA. 
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5. Melting the CGC in Nucleus-Nucleus and 
Proton-Nucleus Collisions 

In general, the problem of high energy hadron-hadron collisions is very 
difficult. As we discussed in the introductory section, a full understand- 
ing of multi-particle production in high energy collisions is one of the 
outstanding problems in high energy QCD. The approach developed in 
previous sections however suggests that classical methods may be use- 
ful in studying multi-particle production in high energy collisions. In 
this section, we will discuss applications of this approach in describing 
the initial stages of very high energy heavy ion collisions, proton (or 
deuteron)-nucleus collisions and briefly, peripheral nucleus-nucleus col- 
lisions. 

Understanding the initial conditions for heavy ion collisions is of 
the utmost importance in studying the high energy heavy ion experi- 
ments performed at the Relativistic Heavy Ion Collider (RHIC) at BNL 
and (to be performed) at the Large Hadron Collider (LHC) at CERN. 
The initial conditions are crucial in determining whether the quark- 
gluon matter produced immediately after the collision will equilibrate 
to briefly form a thermalized Quark Gluon Plasma (QGP). Classical 
methods are applicable because the occupation numbers of gluons in the 
nuclear wavefunctions (and immediately after the collision) are of order 
1/as- One can therefore study the space-time evolution of partons which 
are "frozen" in the Color Glass Condensate through their release in the 
collision and their subsequent evolution. Unfortunately, it is not possible 
in this approach to follow their evolution all the way to equilibration. 
This is because the gluon occupation numbers become small well be- 
fore equilibration thereby signalling the breakdown of the classical field 
approach. 

Since the classical field equations are fully non-linear, they cannot be 
solved analytically The solutions can however be determined numer- 
ically and the initial energy and number distributions can be obtained. 
A particular feature of our formulation of the scattering problem is that 
the classical Yang-Mills field equations are boost invariant-the equations 
are independent of the space-time rapidity. The boost invariance leads 
to a significant simplification-the classical Yang-Mills equations are only 
functions now of the two transverse spatial directions and the proper- 
time. 

In order to maintain gauge invariance, the numerical problem is for- 
mulated on the lattice and solved as a function of proper time. In the 
following, we will discuss briefly the numerical procedure and the results 
of the simulations. The non-perturbative results can be compared to the 
results of analytical computations at large transverse momenta and we 
will show that they agree with these computations. 



"For recent analytical work in this direction, see Ref. 
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The results of these numerical computations can now be interpreted 
in light of the recent data from RHIC. We will discuss some of the 
phenomenological consequences of our results-in particular the implica- 
tions for the importance of final state interactions of partons (beyond 
the initial stage of high occupation numbers) and their subsequent ther- 
malization. Quantum effects, such as the non-trivial geometric scaling 
discussed earlier, may also play an important role in interpreting the 
very interesting and in some instances very puzzling RHIC data. 

As we will discuss, a contentious issue in deciphering the RHIC data 
is whether the phenomena observed are primarily due to initial state or 
final state effects. Experiments underway at RHIC on Deuteron-Gold 
collisions should help settle the issue since one expects final state effects 
in these collisions to be relatively unimportant ° Recently, there have 
been several papers studying various final states at RHIC energies in 
the CGC model. We will discuss these works and their predictions for 
the RHIC data. Finally, we will also discuss briefiy the predictions of 
the CGC model for particle production in peripheral nucleus-nucleus 
collisions. 

5.1. Classical Picture of Nuclear Collisions 

The classical picture of nuclear collisions was first formulated by 
Kovner, McLerran and Weigert 1^ in the framework of the McLerran- 
Venugopalan model for a single nucleus. At very high energies, ^ oo 
of one of the nuclei (and P~ —> cx for the other), the hard valence quark 
(and gluon) modes are highly Lorentz contracted, static sources of color 
charge for the wee parton, Weizsacker-Williams, modes in the nuclei. 
For simplicity, we will first consider only central collisions of cylindri- 
cal nuclei with uniform matter distributions. The case of realistic nuclei 
and non-central collisions will be discussed later. The valence sources are 
then described by the current 

r-^'ir^) = 5-'+pUr±)Six-) + S''-p^{r^)S{x+) , (5.1) 

where Pi(2) correspond to the color charge densities of the hard modes in 
nucleus 1 (nucleus 2) respectively. The 5-function in for the valence 
parton current of one nucleus (or 5-function in x'^ for the other) implies 
that we are literally assuming the nuclei to move at the speed of light. In 
reality of course this condition should and can be relaxed. We will discuss 
this issue later on in this section. In the collision, these valence partons 
are assumed to be Eikonal sources-they continue on their straight line 



"Scattering effects that are often thought of as "final state" effect s in proton- nucleus 
effects can, in a different gauge, be interpreted as initial state effects l2Zl. The final state 
effects emphasized here are those involving re-scattering of on-shell partons off each 
other. 
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trajectories along the light cones. They may acquire a phase due to the 
rotation of the color charge in the scattering . However, in the gauge 
we will use, no such phase will appear in our treatment of the classical 
scattering problem. 

For each color configuration in each of the nuclei, the classical field 
describing the small x modes in the Effective Field Theory (EFT) is 
obtained by solving the Yang-Mills equations in the presence of the two 
sources. We then have 

D^.F'"' = r . (5.2) 

To compute a physical quantity (O) , the gauge field configurations have 
to be averaged over the respective Gaussian path integrals of the two 
nuclei, 

[P?(rx)p?(r-±l+pi(r-x)p^(ri)]' 



(0)p = / dpidp2 0(pi,P2)exp ( - / rf^,^ i^iv-x;^u-x;^^P2V-x;P2V-^;j ) (5 3) 



ON^log(^|. (5.4) 



For instance, the small x gluon distribution is simply related to 
the Fourier transform Af{kj_) of the solution to Eq. ( 15.21 by < 
>P) where the subscript denotes the average above. Here 
have assumed identical nuclei with equal Gaussian weights g^^\, where 
aji is the average color charge squared of a nucleus, defined in Ea. 12.51 

2 4 2 

We will henceforth use the variable Ag = g It is simply related to 
the saturation scale Qs by the relation 

, ^QCD J 

In practice. As ~ Qs- 

In general, one can make the following ansatz for the gauge fields 
as a function of proper time r = \/2x+x^ in the different light cone 
regions: 

A' = a3{r,XT)0{x~)e{x+) + a\{T,XT)0{x~)e(-x^) 

+ al(r,XT)e{-x-)9ix+), (5.5) 
= ±x^a{T,x^)e{x^)e{x~^). (5.6) 

where Q3 and a are respectively the transverse and longitudinal compo- 
nents of the gauge field the forward light cone while ai 2(1' ±) (* = li 2 are 
the two transverse Lorentz indices) are pure gau ge fie lds defined through 
the gauge transformation parameters Aq(?7,r_L)' — ^ 



Ql.2(''i) = T -Pe ±''pr°j V Pe ±''pr°j . (5.7) 



Here r\ = ir^pi-oj =p log(x^/xJ^^j) is the rapidity of the nucleus moving 
along the positive (negative) light cone with the gluon field Q^x(2'\- 
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Ai 2(')j'"±) a-re determined by the color charge distributions Aj_Aq = pq 
(q=l,2) with Aj^ being the Laplacian in the perpendicular plane. 
We work in the Fock-Schwinger (or "radiation") gauge 

= x+A^ +x^A+ =0, 

which is the interpolation between two light cone gauges. Fixing this 
gauge however does not fix the gauge completely. The residual gauge 
freedom can be fixed by imposing the Coulomb gauge condition Vj^ • 
= in the two transverse dimensions. 

Substituting Eq. 15.51 in Eq. 15.211 . and re- writing the equations in 
terms of the transverse coordinates x_i_, the proper time r and the space- 
time rapidity rj, one observes that the Yang-Mills are independent of rj, 
namely, they are boost invariant. The fields q| and a are functions of 
xj_ and T only. The boost invariance of the solutions of the Yang-Mills 
equations is entirely due to our approximation-using 5-function sources. 
Smearing the sources in rapidity would destroy the boost-invariance of 
the solutions. Nevertheless, one expects that the boost-invariance ap- 
proximation is a reliable one especially at central rapidities in a nuclear 
collision. This point will be discussed further. 

The initial conditions for the solution of the Yang-Mills equations 
in the forward light cone are determined by matching the solutions in 
the space-like and time-like regions at r = 0. Requiring that the gauge 
fields be regular at r = 0, the coefficients of the singular pieces of the 
equations -D^i-F'" = and D^±F'^^ = (for x~ ,x'^ ~* 0) have to be 
set to zero. These give the boundary conditions at r = 0: 

al{0,x^) ^ a\{0,x^)+a'2{0,x^), (5.8) 

a{0,x±) = ^[a{{0,x^),al{0,x^)]. (5.9) 

These conditions, first formulated for infinitely large nuclei, are the same 
for finite nuclei as well. The boundary conditions remain the same even 
when the fields a\ 2 before the collision are smeared out in rapidity prop- 
erly account for singular contact terms in the equations of motion Qj^. 
Further, since the equations are very singular at r = 0, the only condi- 
tion on the derivatives of the fields that would lead to regular solutions 
are ^ralT^o, Sraxlr^o = 0. 

Perturbative solutions of the Yang-Mills equations to order in 
the color charge density (or equivalently to second order in As/k^) were 
found, and at late times, after averaging over the Gaussian sources, the 
number distribution of classical gluons was found to be 



29 116..115. 



dN 2 2gV jVc(jVg-l) 

= x4 71 L(k±,X), (5.10) 



dyd'^k^ (27r)4 fc4 

where L{kj^,\) is an infrared divergent function at the scale A. This 
result agrees with the quantum bremsstrahlung formula of Gunion and 



Bertsch' 
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The distributions are very sensitive to L{k^, A). What is novel about 
the classical approach is that, at sufficiently high energies, the non- 
linearities in the Yang-Mills fields self-consistently regulate this infrared 
divergence. To confirm this claim, one needs to solve the Yang-Mills 
equations to all orders in As/kj_. This is very difficult to do analytically. 
Fortunately, as we will discuss in section 5.2, the classical problem can 
be solved numerically. 

The discussion here has been extended to treat the collision of finite, 
ultrarelativistic nuclei with realistic nuclear matter distributions ' — 
For the case of finite nuclei, the issue of color neutrality, as discussed 
in sections 2.4 and 3.10, becomes very important. A global neutrality 
constraint at the level of the nucleus is insufficient to ensure that there 
are no large field strengths outside the nuclear radius-color neutrality 
must be imposed at the nucleon level. The practical implementation of 
color neutrality will be discussed further in the following. 

As mentioned earlier, a major simplification occurs in the classical 
approach when boost invariance is assumed. It is likely a good approx- 
imation, especially in the central region, but the effects of relaxing this 
condition are unclear. A first step in answering this question (before 
performing a fully 3-|-l-dimensional simulation) is to study the stabil- 
ity of the 2 -t-l-d imensional results with respect to a perturbation in the 
//-direction 

5.2. Numerical Gluodynamics of Nuclear Collisions 

Classical real time numerical solutions of gauge theories were first dis- 
cussed in the context of sphaleron transitions during the electroweak 
phase transition in the early universe Similar techniques can be ap- 
plied to discuss the problem at hand. It is most convenient to follow 
the Hamiltonian approach, namely, to construct the appropriate lattice 
Hamiltonian and solve Hamilton's equations on the lattice with the lat- 
tice analog of the initial conditions discussed in section 5.1. These will be 
discussed below. We will subsequently discuss the results of our classical 
numerical simulations. 

i) Numerical Solution of the Yang-Mills Equations 

The QCD action for gauge field in the r, 77, r coordinates reads 



Sqcd = / rdr]dT(fr 



4Tr(g'^"g''/'F„^F^.) 



(5.11) 



where F^y — dfi Ay — dyAfx — ig[Af_i, Ay] and the metric is diagonal with 
gTT ^ _gxx ^ _gyy ^ g^^^ gTiri = „ 1/^-2 = ^a^a^ and t° represent 

a gauge group matrices with the normalization of Tr(t^t'^) = 25ij. The 
Lagrangian density in = gauge is 

£ = Tr ('{{drA.f + ^{drA.f ^F?,i) , (5.12) 
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where i,j runs over transverse coordinate x and y. 

Now let us assume rj independence of the fields. As discussed pre- 
viously, the Yang-Mills equations have this property if the sources are 
strictly (5-function sources on the light cone. We have 

A,{r,i^,r)=MT,r), A„(r, 77, r) = <l>(r, r) , (5.13) 

resulting in F^^ = — 1)^$, where = di — ig[Ai, ■ ■ ■] is the covariant 
derivative in the adjoint representation. Defining the conjugate momenta 
Ei — rdrAi and pr/ = ^drAr^, one finds that the boost invariant Yang- 
Mills Hamiltonian i s th e QCD Hamiltonian in 2+1 dimensions coupled 
to an adjoint scalar 

H-J d'rTr{^E? + Ip2 + lFg + ^(Di<i>)^} . (5.14) 

In order to realize numerically the solutions to the equations of mo- 
tion in the previous section, while maintaining the gauge symmetry, we 
introduce the link variables at the site i 

Uj^i ^ exp[iaAj{i)] , {j ^ x,y), (5.15) 

where, a is a lattice spacing. Defining the plaquette Uu = 
UijUjyi j^lUj j^^U^ J, the Hamiltonian on the lattice is 

Hl = ^ Yl Tri5,2 + ^^(iVc-ReTrl7n) 

where the convention for the generators of the SU(3) color group is 
Trir^'r'') = 2^''^ For g = 2, one obtains the correct normalization of the 
Hamiltonian in the continuum limit. Lattice equations of motion follow 
directly from H]^ of Eq. 15.161 For any dynamical variable v, with no 
explicit time dependence, v — {H]^,v}, where v is the derivative with 
respect to t, and {} denote Poisson brackets. We take Ei, Ui, Pj, and $j 
as independent dynamical variables, whose only nonvanishing Poisson 
brackets are 

{Pi, ■I'i} = SijSab; {Ef, Urn} = -i&hnUlo"-; {Ef , E^} = 2Si,n£abcEl 

(no summing of repeated indices). The initial conditions for the trans- 
verse gauge field and the adjoint scalar field on the lattice can be ob- 
tained in complete analogy to the procedure followed in the continuum. 
The details of this procedure and the expression for the initial conditions 
can be found it Ref. ffl. 

We impose periodic boundary conditions on an A'^ x A'' transverse 
lattice, where N denotes the number of sites. The physical linear size of 
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the syst em is L = aN, where a is the lattice spacing. It was shown in 
Ref. I^that numerical computations on a transverse lattice agreed with 
lattice perturbation theory at large transverse momentum. 



ii) Numerical Method for Finite Nuclei 

In early studies l ^^ l ^^ l '^'j lMI nuclear collisions, for simplicity, were 
idealized as central collisions of infinite, cylindrical nuclei. The color 
charge squared was taken to be a constant for the uniform cylindri- 
cal n ucle i. Furthermore, color neutrality was imposed only in a global 
sense namely, the color charge distribution over the entire nucleus 
was constrained to be zero. While very useful in obtaining first estimates 
of the space-time evolution of the produced gluonic matter, these stud- 
ies did not make predictions for realistic nuclear collisions. In addition, 
studies of the distributions in peripheral collisions, in particular of the 
azimuthal anisotropy associated with elliptic fiow, require finite nuclei 
and realistic nuclear matter distribiitjonsmthin each nucleus. These re- 
quirements were discussed in Refs. ^ lo|119| 

The problem with color neutrality for a finite nucleus can be stated 
simply as follows. If we impose the simple and obvious constraint that the 
color charge distribution must be zero outside the nucleus, the solution 
of Poisson's equation can still give a non-zero gluon distributions outside 
the nucleus. In two dimensions, the fall-off of the gluon field is rather slow 
as shown in Fig. 15.21 This slow fall-off is a problem for a finite nucleus 
since the gluon field is associated with a non-zero field strength. Clearly 
the simple prescription for color neutrality is not sufficiently stringent. 



io' 

10' 



MV 

color neutral II 

nn^=1 fm"' 



2 4 6 8 10 12 
r(fm) 

Fig. 23. Gluon field as a function of radial distance. Original MV model is shown by 
circles, while squares correspond to the Color Neutral II prescription (see text). The 
triangles represent results from solution of Poisson's equation with a screening mass. 
The results are for AsqR = 37. 



A more realistic prescription would be to apply the color neutrality 
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constraint already at the nucleon level. Our numerical procedure to im- 
plement the constraint for finite nuclei is as follows. We first sample A nu- 
cleons on a discrete lattice requiring that they satisfy a Woods-Saxon nu- 
clear density profile in the transverse plane. Note that this procedure gen- 
erates the same distribution in the continuum as As(a;^) = A'^QTji{x±) 
where Ta{x±) = dz Hi{r) is a thickness function, x±^ is the trans- 
verse coordinate vector (the reference frame here being the center of the 
nucleus), k(t") is the Woods-Saxon nuclear density profile, and A^g is 
the color charge squared per unit area in the center of each nucleus. The 
only external dimensional variables in the model are A^o and the nuclear 
radius R. 

Next, Gaussian color charge distributions are generated on the lat- 
tice. The probability distribution of color charge in a nucleon is expressed 
as 



P[p] = exp 




(5.17) 



where A^ ^ is the color charge distribution squared, per unit area, of 
a nucleon at a lattice site j and N is the number of lattice sites that 



comprise a nucleon. A„ ^ is obtained from A^g by assuming that the 
color charges of the nucleons add incoherently. There are two versions 
of the subsequent step. In the first (which we term Color Neutral I), we 
subtract from every pj the spatial average pj/N in order to guarantee 
color neutrality (p) = for each nucleon. In the second (termed Color 
Neutral II), the dipole moment d of each nucleon is eliminated in a 
similar manner. 



^ 1.4 
ra, 1.2 

2 1 

^ 0.8 
0.6 
0.4 
0.2 





• MV 




■ color neutral 1 
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Fig. 24. Color Charge Correlator in momentum space. Original MV model is shown 
by circles, while squares correspond to the result from color neutrality condition I and 
triangles correspond to color neutrality condition II (see text). The color charge correlator 
is plotted versus ™ units of Asq. 
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In Fig. 15.21 we plot the Fourier transform ol the charge correlator, 
which in the continuum is defined as 



for the MV model and for the two variants which impose color neutrality 
on the nucleon level. In the MV model, p(p±) is a constant everywhere 
except at pj^ = where it is constrained to be zero from the global charge 
constraint. In the Color Neutral I (II) variant, we see that p{p±) ~ Pj_ 
(~ pj_) for small momenta pj_ < A^q and is constant at larger momenta. 
The oscillatory behavior seen for Color Neutral I and II is due to the fact 
that the correlator in coordinate space is not strictly a delta-function. 
In the coordinate space the charge correlator for the two models. Color 
Neutral I and II, falls off rapidly, as ~ l/a;j^ and ~ l/2;j_ respectively, 
at larger distances. 

It is an interesting coincidence that the behavior of p{p±) in our 
model is similar to the behavior expected from the renormalization group 
(RG) evolution of color charges in the McLerran-Venugopalan model as 

PTB 1 2S 

discussed here in Section 3.10. In Refs. , it is shown that the 

screening of color charges due to the RG evolution gives a behavior 
p(p_l) ^ p\ for pj^ < Aso (and p(pj^)=constant for p^ > A^o). 

iii) Numerical Results for Distributions in Energy and Number 
in Central Collisions 

We will now discuss results for central collisions of very large cylin- 
drical nuclei. We will see later that they are not very different from 
those for central collisions of realistic nuclei. There are only two free 
parameters for the problem of nuclear collisions as formulated int the 
EFT. One is the saturation scale As while the other is the nuclear ra- 
dius R f (For an infinite cylindrical nucleus, one has 7ri?^ = L^, where L 
is the lattice size.) Any dimensional quantity P, well defined within the 
EFT, can be written in terms of the physically relevant parameters As 
and R as As fp{AsR), where d is the dimension of P. The non-trivial 
physical information is therefore contained in the dimensionless function 
fp{AsR). On the lattice, P will generally depend on the lattice spac- 
ing a; this dependence can be removed by taking the continuum limit 
o — > 0. The broad range of physically relevant values of As for RHIC and 
LHC energies are ~ 1-2 GeV and 2-4 GeV respectively-corresponding to 
As-R « 30-120 approximately'^. Also, for central Au-Au collisions, we 



P Strictly speaking, the saturation scale i s not a free parameter since it can be determined 
from the gluon density as in Ref. I' I ' I However, for the momentum scales of interest, 
there is much uncertainity in the gluon density. Moreover, for nuclei, gluon shadow- 
ing contributions are not under control. Our results will therefore be non-perturbative 
formulae valid for a wide range of . 



'^li we extrapolate from the Golec-Biernat HERA parametrization, we get As = 1.4 GeV 
for RHIC and As = 2.2 GeV for LHC. 




(5.18) 
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obtain L = 11.6 fm as the physical linear dimension ol our square lat- 
tice. 
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Fig. 25. (a) er/Af as a function of tAs for A^R = 167.4. (b) er/Af as a function of 
Asa for As_R = 167.4 (squares) and 50(circles), where a is the lattice spacing. Lines are 
fits of the form a — bx. 



For the transverse energy of gluons we get, on purely dimensional 
grounds, 

The function is determined non-perturbatively as follows. In Fig- 
ure E21 (£'■)) we plot the Hamiltonian density, for a particular fixed value 
of AsR = 167.4 (on a 512 x 512 lattice) in dimensionless units as a func- 
tion of the proper time in dimensionless units. We note that er converges 
very rapidly to a constant value. The form of er is well parametrized by 
the functional form er = a + /3exp(— 7t). Here dEx / dri / n B?' = a has 
the proper interpretation of being the energy of produced gluons per 
unit area per unit rapidity, while = l/7/As is the "formation time" 
of the produced glue. 

In Figure 15.21 (b) , the convergence of a to the continuum limit is 
shown as a function of the lattice spacing in dimensionless units for two 
values of As-R. In Ref . this convergence to the continuum limit was 
studied extensively in an SU(2) gauge theory for very large lattices (up 
to 1024 X 1024 sites) and shown to be linear. The trend is the same for 
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the SU(3) results-thus, despite being further from the continuum Umit 
for SU(3) (due to the significant increase in computer time) a hnear 
extrapolation is justified. We can therefore extract the continuum value 
for a. 

Numerical results for the physically rele vant RHIC and LH C initial 
conditions were discussed in several papers I '^^ l ^^-^ l ^^'^ l ^^^ l ^-*^'^ ] Very re- 
cently, Lappi pointed out that the over all norm alization in Eg. 15.191 was 
inc orrect and the f r computed in Refs. I^^ l ^^l ig too large by a factor of 

two' ' As we will discuss below, the normalization for the number 

is correct. Wherever feasible, the corrected results will be discussed be- 
low. We find fsi^O) = 0.269 and /b( 167.4) = 0.248. The RHIC and LHC 
values likely lie in this wide range of AsR. The SU(2) value is approxi- 
mately half the SU(3) value. Note that the variation of as a function 
of AsR is extremely weak. The formation time tu = l/7/As is essen- 
tially the same for for both the SU(3) and SU(2) cases-for AsR = 167.4, 
7 = 0.362 ± 0.023 and . As discussed in Ref. it is ~ 0.3 fm for RHIC 
and ~ 0.13 fm for LHC (taking As = 2 GeV and 4 GeV respectively-for 
the values of As extracted using the Golec-Biernat fit, the corresponding 
times are of course larger). 

We now combine our expression in Eq. with our non- 

perturbative expression for the formation tim e to obtain a non- 
perturbative formula for the initial energy density 

s=^Ai (5.20) 

This formula gives a rough estimate of the initial energy density, at a 
formation time of r^) = l/'y/AsR where we have taken the average value 
of the slowly varying function 7 to be 7 ~ 0.3. For As = 1.4 GeV, one 
obtains e ~ 10 GeV/fm^. 

We now report our results for the initial multiplicity of gluons pro- 
duced at central rapidities. First consider a free field theory whose Hamil- 
tonian in momentum space has the form 

Hf = lj2 (\n(k)f+c.\k)mk)f) , (5.21) 



where 0(fc) is the fcth momentum component of the field, n{k) is its 
conjugate momentum, and uj{k) is the corresponding eigenfrequency. The 
average particle number of the fe-th mode is then 



N{k) = io{k){\mf) = ^Jmk)\'\Ak)\^), • (5.22) 

In our case, the average () is over the initial conditions. 

We use two different generalizations of the particle number to an 
interacting theory. We have verified that the two definitions agree in the 
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weak-coupling regime at late times Our first definition is b ased on 
the behavior of a free-field theory under cooling. We obtain 1^ 

where t is the cooling time (not to be confused with real or proper time), 
and V{t) is the potential energy of the relaxed free field after cooling. The 
relaxed potential V{t) is gauge- invariant-hence so is this definition of the 
particle number. This is an attractive feature of the cooling method. Un- 
fortunately, it presently only permits determination of the total particle 
rmmber and cannot be used to find the number distribution N{k±). 

Our second definition of the multiplicity will enable us to compute 
N{k^). We impose the Coulomb gauge condition in the transverse plane, 

• A±_ — 0, and substitute the momentum components of the resulting 
field configuration into Eq. 15.221 . We can determine N{k^_) from the 
rightmost expression of Eq. 15.221 : the middle expression of Eq. 115.221 
can then be used to obtain uj{k±^). 

In Fig. 26(a), we plot the normalized gluon transverse momentum 
distributions versus kj_/As with the value A^-R — 167.4, for both the 
SU(3) and SU(2) gauge theories " . Clearly, we see that the normalized 
result for SU(3) is suppressed relative to the SU(2) result in the low 
momentum region. In Fig. 26 (b), we plot the same quantity over a wider 
range in k±/As for two values of AsR- At large transverse momentum, 
we see that the distributions scale exactly as — 1, the number of 
color degrees of freedom. This is as expected since, at large transverse 
momentum, the modes are nearly those of non-interacting harmonic 
oscillators. At smaller momenta, the suppression is due to non-linearities, 
whose effects, we have confirmed, are greater for larger values of the 
effective coupling AsR. The SU(3) gluon momentum distribution can be 
fitted by the following function, 

^ ''^ \fnikT/As), (5.24) 



ttR'^ dricPk^ g 




k\ + m2/Tcff) - l] ' (fcT/A. < 1.5) 



(5.25) 



[kx/As > 1.5) 



with ai = 0.118, m = 0.034As, T^s = 0.47As, and a2 = 0.0087. At low 
momenta, the functional form is approximately that of a Bose-Einstein 



"^Our results for the number distribution ag ree w ith Ref. As As/2 in the ensuing 

discussion of the px distributions. See Ref. 11221 for details. 
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Fig. 26. Transverse momentum distribution of gluons, normalized to the color degrees 
of freedom, n{kj_) = /„/(Af^ - 1) (see Eq. I^H}) as a function of AsR for SU{3) 
(squares) and SU(2) (diamonds). Solid lines correspond to the fit in EQ. i5.25i . 



distribution in two dimensions even though the underlying dynamics 
is that of classical fields. The functional form at hig h monientum is 
motivated by the lowest order perturbative calculations ^^^^^^^1 

Integrating our results over all momenta, we obtain for the gluon 
number per unit rapidity, the non-perturbative result. 



1 dN 



--o = ^fN{AsR)Al. 
9^ 



(5.26) 



We find that fiv (167.4) — 0.3. Our results for fiu are in reasonable 

agreement with those of Lappi' ' '. We have checked for an SU(2) 

gauge theory that the results for a wide range of As-R vary on the order 
of 10%. The results from the Cooling and Coulomb methods also show 
very good agreement (less than 10%) especially for larger values of AsR- 
If we take the ratio of Eq. (15.1911 to Eq. 15.261 , we find that the initial 
transverse energy per gluon is dE^ /dri/dN/r]\^—Q = As = 0.88 A^. If 
we take As — 1.4 GeV, we find that the energy per gluon is ~ 1.23 CeV- 
about a factor of 2 larger than the value for charged hadrons measured 
at RHIC. 

The topological charge generated in the initial stages of a heavy 
ion collision can also be computed in the classical CGC framework. An 
interesting result is that if strict boo st inv ariance is a good assumption 
sphaleron transitions are suppressed This is because the Chern- 
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Simons number in this case is invariant under all rapidity independent 
gauge transformations. The primary mechanism for the generation of 
topological charge at the early stage is then by fluct uations of the color 
electric and magnetic fields. It was shown in Ref.^^that the topological 
charge generated in this manner, at the early stage, are small. These 
results may be relevant for the formation of P and CP-odd metastable 
states in the late stages of heavy ion collisions 

iv) Numerical Results for Centrality Dependence of Energy, 
Multiplicity and Elliptic Flow 

For realistic nuclei, the non-perturbative relations discussed in sec- 
tion 5.2.3 are less simple. One can write Eq. 1^^^ more generally as 

dri po 

where po = /5(0, 0) = 4.321/m~^ and Npart = / (fix^p{b,xj_). For 
explicit expressions for /jv(6) for different values of Aso, see Ref. UJ^. 
One can similarly compute dEg/drj for finite nuclei. 
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Fig. 27. The two components of the transverse pressure {T^x and Tyy) and the energy 
density e plotted as a function of r in dimensionless units. The results are for a impact 
parameter b = R and dimensionless coupling Aso^ = 37. Also shown is the sum of the 
two transverse pressures. 



The azimuthal anisotropy in the transverse momentum distribution 
is a sensitive probe of the hot and den se matter produced in ultra- 
relativistic heavy ion collisions IMIHI. A measure of the azimuthal 
anisotropy is the second Fourier coefficient of the azim uthal distribu- 
tion, the elliptic flow parameter «2. Its definition is^^ 



r ^rf0cos(2^)/rfV3 ^^ 



V2 = {cos(2<^)) = — ,3T' ■ (5-28) 
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The classical Yang-Mills approach may be applied to compute el- 
liptic flow in a nuclear collision. For peripheral nuclear collisions, the 
interaction region is a two-dimensional almond shaped region, with the 
X axis lying along the impact parameter axis and the y direction per- 
pendicular to it and to the beam direction. Even though large electric 
and magnetic flelds (and the corresponding transverse components of the 
pressure in the x and y directions) are generated over very short time 
scales r ~ 1/As, the signiflcant differences in the pressures, responsible 
for elliptic flow, are only built up over much longer time scales r ~ 7i. 
This can be seen in Fig. 27 where we plot the two transverse components 
of the pressure {Txx and Tyy) and the energy density as a function of 
proper time (in units of A^o) for a peripheral nuclear collision. Moreover, 
the elliptic flow is generated by soft modes px ~ As/8. Our result has 
important consequences for the theoretical interpretation of the RHIC 
data-these will be discussed later in the text. 

The elliptic flow, defined by Eq. 15.28II . can be computed, as in the 
case of the gluon multiplicity, in two different ways; directly in Coulomb 
Gauge (CG) and by solving a system of relaxation (cooling) equations for 
the fields. It is easy to show that V2 N, N being the total gluon number, 
can be reconstructed from the cooling time history of Txx — Ty y , just as 
TV can be reconstructed from that of the energy functionally^ 



y^N=^f^ /'°°-^(T^^(t)-T^2^(t)). 

V TT ^0 Vt 

This expression for V2 N is manifestly gauge invariant. 



(5.29) 



^8f 

CM 

2t 



CG v.s.Cooling 




corrected values 
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Fig. 28. The centrality dependence of V2 at the earliest times is computed using cooling 
(open symbols) and CG (filled symbols). Results are for Asoi? spanning the RHIC-LHC 
range, specifically, AsoR = 37 (squ ares'). 74 (triangles), and 148 (stars). Full circles 
denote preliminary STAR data LLaal The band denotes the estimated value of V2 when 
extrapolated to very late times. "Corrected values" denotes the late time cooling and 
CG result for AsoR = 37 at one centrality value. 
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In Fig. 1281 we plot «2 reconstructed from the cooling time history of 
only the potential terms in Txx — Tyy, along with the CG values (also 
including potential terms only) as a function of nch/ntot for different 
values of f^soR discussed in the figure. The systematic errors repre- 
sented by the band (for A^q = 37 are primarily due to limited resources 
available to study the slow convergence of the cooling and CG compu- 
tations. We have studied the late time behavior of V2 for one impact 
parameter-the results are shown in the figure. 

The asymptotic values of v^, as predicted by the model, undershoot 
the data. For a fixed impact parameter, the model predicts that, as 
Kg^R — > oo, the classical contribution to the elliptic flow goes to zero. 
This is because increasing A^o^ is equivalent to increasing R for fixed 
Aso and therefore reducing the initial anisotropy. 

In Fig EH V2{PA_) is plotted for h/2R = 0.75 for AsqR = 148. Our 
calculations show that the elliptic flow rises rapidly and is peak ed for 
p±_ ~ Aso/8 before falling rapidly. The theoretical prediction I 'is that 
for p_L S> Aso, V2{p±_) ~ Agg/p^. The lattice numerical data appear to 
confirm this result-better statistics are required to determine the large 
momentum behavior accurately. The dominant contribution of very soft 
modes to V2 helps explain why the cooling and CG computations differ 
until very late times. The soft gluon modes have large magnitudes and 
therefore continue to interact strongly until very late proper times. Con- 
comitantly, the occupation number of these modes is not small and the 
classical approach may be adequate to describe these modes even at the 
late times considered. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Pt/A, 



Fig. 29. f2(px) as a function of transverse momentum in dimensionless units for 
hsoB^ = 148. 
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Fig. 30. Comparison of the centrality dependence of the gluon distribution from SU{3) 
lattice results to data from experiments H^ ^ l -*^?^ ! The strong coupling constant is fixed 
to the value g^ = 4. The lattice results for As{0) = 1.18 GeV and As(0) = 1.95 GeV are 
multiplied by a factor 2.4 and 1.1, respectively. 



5.3. Melting the Color Glass Condensate at RHIC 

In Fig l5.3l we plot the computed centrality d epend ence of gluons together 
with the experimental data from PHOBOsll^and PHENIX ESfl. We 
assume here that the charged particle multiplicity is two thirds of the 
gluon number. The classical computation is performed for fixed Os; the 
centrality dependence, as seen from Eq. 15.27II . comes from the depen- 
dence of /tv on the impact parameter. In Ref. it was shown that 
/jv = /Ar(Asi?) increases slowly with As i?- hence one expects it to vary 
with impact parameter. We see that the results agree reasonably well 
with the data. 

The centrality dependence of the transverse energy is studied in 
Fig. 15.31 As in the case of the multiplicity, even though the absolute 
normalization is strongly dependent on one's choice of As, the central- 
ity dependence is very similar for the two A^'s and shows reasonable 
agreement with the data. 

Let us now compare our results with those derived previ ously by 
Kharzeev and Nardi in Ref. and discussed further in Ref. In 
these works, one obtains in terms of Qs, the average saturation scale, 
the result 



dN^ 
drj 



= Cat 



d X I 



;{b,xj 



as 



cj^xG(x,Qs{h)) 



part 



(5.30) 



In the leading logarithmic approximation, if c^r is constant, one obtains 
a logarithmic dependence on the centrality entirely from xG{x,Q^{b)). 
One could thus attribute the logarithmic behavior at the classical level 
to fixed as and leading logarithmic behavior of the gluon distribution 
function or equivalently, at higher order, to the one loop running of as ■ In 
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Fig. 31. Comparison of the centrality dependence of the gluon transverse energy dis- 
tribution from SU(3) lattice results to data from experiments 1 ^?? ! The strong coupling 
constant is fixed to the value g^ = 4. The lattice results for Aij(O) = 1.95 GeV are scaled 
by ^ while those for As(0) = 1.18 are scaled by ^. 



the physically interesting regime, it is difficult to distinguish between the 
two. A reasonable agreement with the data is also seen in this formulation 
of the problem. 

The relation between the two formulations is as follows. What jffi 



simulate numerically is the color charge squared per unit area, Aj 

The saturation scale, on the other hand, is a scale determining the 

27 

behaviour of the gluon number distribution ' . The relation between 
the two is given in Eq. EH The relat ion between cj\[ and /tv is then 
cn = 47r^/jv/(A'^ — 1)/ ln(Q^/AQp^). Therefore, if cj^ is to be a con- 
stant, /jv increases logarithmically with Qs- A weak rise in is seen in 
our simulations. If the infrared scale, as argued in Ref. is a number 
of order 0{1/Qs), we would have As ~ Qs- This is indeed the case in 
practice lii^. 

The initial transverse energy per particle is Ej^/N ~ 0.88As, which 
for t he Gole c-Biernat value of As — 1.4 GeV (also the value favored by 
Refs.E^l^, gives Et/N = 1.23 GeV. This is a factor of 2 larger than 

the value for measured charged hadrons^^ -: For this value of As, the ini- 
tial number of gluons is approximately half the multiplicity of hadrons at 
central rapidity. The excess ratio of Ej^/N in the initial state can be re- 
duced by a) inelastic fragmentation of partons, which increases A*' and b) 
thermalization and hydrodynamic expansion expansion which increases 
N and decreases Ex respectively. Thermalization increases the par ticle 
number because it is driven primarily by inelastic partonic processes'^. 
It is very conceivable that the multiplicity increases by a factor of two 
due to either a) or b) or both. Thermalization increases the par ticle 
number because it is driven primarily by inelastic partonic processes'^. 
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The p_L distributions from the melting of the CGC were discussed 
previously. The distributions obtained numerically from purely classical 
considerations can be used as initial conditions in parton transport equa- 
tions which describe the evolution of the system when the description 
in terms of classical fields is inadequate. The parton distributions ex- 
hibit geometric scaling-they are functions of p±/As alone. (At large p± 
(p_l_ >> As), one expects this classical geometric scaling break down and 
distributions to be described by perturbative QCD.) It has been argued 
that the RHIC data demonstrate scaling behavior in a broad kinematic 
region in p± A very similar scaling is seen in the string percolo- 

tion framework P*^^ ! In Ref. '2^, a novel mechanism for nuclear collisions 
was suggested, where, due to quantum evolution of saturation effects in 
the wavefunctions of the incoming nuclei, the pj_ window, where geomet- 
ric scaling holds, may extend to p^ > / Aqcd- This mechanism has 
precisely the same origin as the geometric scaling discussed at length 
in Sect. 3.4 and which was suggested in Sect. 4.3 to be an explana- 
tion the large scaling window observed in Deeply Inelastic Scattering at 
HERA. As a consequence of geometric scaling, the p± distributions in 
AA-collisions at RHIC (and more strikingly at LHC) scale with central- 
ity oc Npart- The boundary of this region is estimated to be ~ 4 GeV at 
RHIC and ~ 9 GeV at LHC. The CGC-|-geometric scaling picture there- 
fore explains the qualit ative behavior of the quenching seen in central 
collisions at RHIC^ I For peripheral collisions, the saturation scale is 
small: the window for geometric scaling shrinks and the binary scaling 
of pQCD is restored just as observed in the da ta. 

This CGC-|-geometric scaling picture of Ref.l^is in striking contrast 
to the QGP-|-energy loss d escription of the quenching of spectra seen in 
central AA-collisions EHl scenarios, at face value, make very 

different predictions for quenching in d-Au collisions. The former predicts 
quenching in roughly the same p± range with the centrality dependence 
(Npart)^ , the latter predicts a "Cronin-enhancement" just as seen in 
lower energy p-A collisions. However, multiple-scattering of the Cronin 
variety is also, in principle, present in the CGC picture. The magnitude 
of this effect must be quantified before we rush to any conclusions about 
the RHIC data. We will discuss p-A collisions further in Section 5.5. 

We now turn to the theoretical interpretation of the RHIC f 2 data 
in the CGC approach. It is clear from Fig. I28l that our result for V2 con- 
tributes only about 50% of the measured V2 for various centralities. Our 
p± distributions also clearly disagree with experiment -OS 1401 Naively, 
one could argue that the classical Yang-Mills approach is only applica- 
ble at early times so additional contributions to V2 will arise from later 
stages of the collision. While there is merit in this statement, it is also 
problematic as we will discuss below. The reason the situation is com- 
plex is as follows. We observed that it takes a long time r ~ i? to obtain 
a significant elliptic flow in the CGC. At these late times, one would ex- 
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pect that the classical approach would be inapplicable due to the rapid 

expansion of the system. On the other hand, we have seen that «2 in 

the classical approach is dominated by soft modes which are strongly 

interacting and don't linearize even at time scales r ~ 7?. Clearly, the 

soft modes cannot be treated as on-shell partons even at times r ~ J?! 

This conceptual problem will be discussed further in Section 5.4. The 
118 

observation in Ref. that the CGC alone is unable to explain the low 
behavior of the RHIC V2 data is, thus far, the strongest argument 

for the importance of final state effects at RHIC. 

Finally, we note that V2 is extracted from a variety of techniques-in 

addition to a rea ction plane method, two and four particle cumulant 

analyses are 

used inn. As we noted previously, the CGC p± distribu- 
tion at high p± decays as while the data is flat up to p± ~ 10 
GeV. Recently, non-flow two particle correlation s in the CGC model 
were suggested as an explanation of the V2 data li^. It is unclear at 
present whether non-flow correlations can explain other features of the 
measured azimuthal anisotropy. In the approach reported previously, a 
procedure very similar to the experimental approach c an be followed and 
two and four particle correlations can be determined It is not yet 
clear whether the results of this study would affect our above conclusion 
about the relevance of final state effects in heavy ion collisions. 

5.4. Equilibration and the Quark Gluon Plasma 

An outstanding issue in the physics of heavy ion collisions is whether 
the hot and dense matter formed equilibrates to form a Quark Gluon 
Plasma (QGP). Hydrodynamic models, which assume local equilibrium, 
have been used to study the evolution of th e QGP and are successful 
in describing some of the data at RHIC^^. It is however not easy to 
show from first principles that thermalization does occur in the extreme 
conditions of a high energy heavy ion collision. Even though the initial 
density is very high, the system is expanding and becoming dilute very 
rapidly; moreover, if the momenta of the gluons is large, asymptotic 
freedom dictates that the cross-section of the gluons is small. 

In the CGC approach, we have argued that the energy densities 
right after the collision are very high-£ ~ 10 GeV/fm^ at RHIC and 
e ~ 60 GeV/fm^ at LHC energies. These energy densities are formed 
very rapidly after the coUision-about Tf^j-^i^n^j^ — 0.4 fm for RHIC and 
'''formation = 0.25 fm at LHC usiug the values of Qs from the Golec- 
Biernat parametrization. At the high energies of interest, a heavy ion 
collision in the classical picture is boost invariant The dynamics of 
particle production at early times is therefore purely transverse. Consider 



°This assumption is modified by quantum effects-however, at central rapidities, boost 
invariance is a good approximation. 
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a narrow slice in space-time rapidity, around the central rapidity rj = 
0. Partons with any significant longitudinal momentum pz will not be 
found in this slice (over time scales of interest )-they have only transverse 
momentum pj_ ~ Qg- Our classical dynamics only describe the initial 
transverse dynamics of partons. In Fig. 27, we see that at a time of 
~ 20Aso''" (~ 2.8 fm), e = Txx + Tyy. This suggests that the partons 
that were strongly interacting in the transverse plane are free-streaming 
by this time. 

It is unlikely however that the classical approach is valid at these late 
times. Shortly after partons are formed ~ 1/Qs, they begin to scatter 
off each other and off the transverse plane. The physics of this process 
is not contained in the classical picture. A first guess is to see if small 
angle elastic 2 2 scattering will lead to thermalization. Analytical 
estimates^^ suggest and numerical solutions of the Boltzmann equation 
confirm that the time scale for equilibration due to small angle 
scattering is ~ exp{l / y/a2) / Q s parametrically. In practice, this time 
scale is not much smaller than the lifetime of the system. Also, as pointed 
out in Ref.ElHl, the pressure computed from these simulations is smaller 
than that of an equilibrated QGP. 

It was suggested by Baier et al. 1^ that inelastic 2^3 scattering 
may provide the dominant mechanism driving the system towards equi- 
librium. Naive power counting suggests that this process is suppressed 
relative to elastic scattering by a power of Qs-however, the 2 — > 3 process 
is more efficient in re-distributing momenta and changing particle num- 
ber. This scenario for thermalization, termed "bottom-up" on account 
of our particular initial conditions, can be outlined briefly as follows. As 
we discussed previously, the classical fields can be linearized into partons 
for T > {Qs)~^ ■ The parton distribution can be expressed as in Eg. 15.241 
with the functional form of the distribution given in Eq. 15.251 This de- 
scription, in weak coupling, is valid (parametrically) up to a proper time 
T ~ {as)~^^'^ IQs- This is the time at which, on account of the longitudi- 
nal expansion, the occupation number of partons / < 1. In a small slice 
in rapidity, pz ~ 1/t which is small for r > 1/Qs- However, pz is built 
up gradually through random multiple scattering: = m^N^ijii , where 
is the Debye mass and Ncoii is the collision rate. Simple estimates 
of these give self-consistently pz ~ Qs/{Qst)^^^ . Thus even though pz 
decreases with r, it does so at a slower rate than given by free streaming. 
The occupation number of the hard gluons (with p±^ ~ Qs) is then 

dNy, Qf 1 (Q.r)l/3 

(IzSaQIpz ^ CXs(Qst) Ql Qs ' ^ ' ' 

3/2 

which gives / < 1 for t > /Qs- This estimate was made assuming 

one had only elastic 2 — » 2 scattering. There is also inelastic 2 — > 3 
scattering going on but the number of these soft gluons Ns « Nf^ at 
these times. When / < 1, this inequality continues to hold but the soft 
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gluons begin to dominate the contribution to the Debye screening mass. 
One obtains A'^s ~ A'^^ at QsT = (as)"^/^. 

For QsT > (as)~^^^, most of the gluons are soft. These coUide fre- 
quently with each other and achieve an equilibrated distribution with a 
temperature T. However, there are still some hard gluons left which con- 
tinue to transfer energy to the soft gluons thereby acting as a heat bath. 
Indeed, as a consequence, the temperature of the soft gluons increases 
linearly initially even though the system as a whole is expanding. The 
temperature finally stops increasing when the hard gluons have lost all 
their energy. This happens at a time QsT = {as)~^'^'^ at a temperature 

2/5 —1/3 

of T = as Qs- The temperature subsequently decreases as r ' as 
one would expect for a fluid undergoing one dimensional expansion. 

The b ottom-up scenario is an attractive one and the authors in 
Ref. ED have show that the centrality dependence of the RHIC ex- 
perimental data can be fit successfully. Nevertheless, a fully successful 
application of the idea to phenomenology is still remote. For the various 
stages of this scenario to be realized, as must be very small-likely, much 
smaller than may be realized at RIfIC and perhaps even LHC energies. 
In the bottom-up scenario itself, the final temperature and the equili- 
bration time are not determined up to a constant-however this number 
can in principle be determined within the theoretical framework itself. A 
number that needs to be determined externally is the liberation coeffi- 
cient C]\[. As we described previously in section 5.2.3 and 5.2.4, we obtain 
C]\[ = 0.3 — 0.5 depending on what Qs is. From empirical cons iderations, 
a larger value cjy ~ 1 is favored in the bottom-up picture^^. Even for a 
favorable choice of parameters which fit the data, one finds r^quU ~ R/2 
where R = 6.8 fm is the radius of a Gold nucleus. While still smaller 
than the size of the system, a necessary condition for thermalization, it 
is not smaller by a large enough margin to make it appear inevitable. 

The long equilibration times discussed here are problematic for un- 
derstanding the RHIC data on V2- Hydrodynamical models that fit the 
data require early thermalization times of t^^^^^i ~ 0.6 fm. The V2 gener- 
ated by the CGC alone is not sufficient to explain the data. The correct 
way to treat the theoretical problem may be as follows. Hard modes 
with > Ks linearize on very short time scales r ~ 1/Qs- Their sub- 
sequent evolution is treated incorrectly in the classical approach, which 
has them free streaming in the transverse plane. In actuality, as dis- 
cussed here, they are scattering off each other via elastic gg gg and 
inelastic gg <- > ggg collisions which drive them towards an isotropic dis- 
tribution This dynamics would indeed provide an additional pre- 
equilibnum contribution to V2 and is calculable. An effect to consider 
here would be the possible screening of infrared divergences in the hard 
scattering by the time dependent classical field. More complicated is the 
effect of these hard modes on the classical dynamics of the soft modes 
and on their possible modification of the contribution of the latter to V2 ■ 
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One has here a little explored dynamical analog to the interplay of hard 
particl e and s oft classical modes in the kinetic theory of Hard Thermal 
Loops 1 ^1 1 147 1 The overlap between the classical field an d B oltzmann 
pictures in this context has been discussed recently in Ref. ^33 A practi- 
cal issue of interest is to use the results of the classical field simulations 
as the initial input to parto n cas cade models which simulate the later 
stages of heavy ion collisions EH 

5.5. p-A & Pempheml A A collisions at RHIC and LHC 

Proton-Nucleus collisions at RHIC and LHC will provide an excellent 
probe of the properties of the CGC. Deuteron-Gold collisions with the 
center of mass energy of -^s = 200 GeV/nucleon at RHIC started in 
January 2003 and Proton-Nucleus collisions are an im port ant part of 
the program at the upcoming LHC colUder at CERN 1^3. A big dif- 
ference between p/D-A collisions and A- A collisions is that final state 
interactions are more important in the latter than in the former. Thus 
p/D-A collisions provide an important benchmark to disentangle novel 
phenomena such as the Quark Gluon Plasma, which arise as a conse- 
quence of strong final state interactions, from initial state phenomena 
which may result from the Color Gluon Condensate. 

A case in point is the remarkable observation of the suppression of the 
single particle spectra, as a function of p_L , in An- Au collisions relative 
to that in pp-collisions (per binary collision) at ^/s = 130 GeV/nucleon 
and 200 GeV/nucleon . An interpretation is that high pj^ partons 
traversing a, QG P suffer significant energy loss resulting in fewer high 
p±^ partons' — ^. If correct, the observed suppression is evidence that a 
QGP has been created in high energy heavy ion collisions. In contrast, 
in p/D-Au coUsions, it is argued that one will see the Cronin effect- 
the ratio of d-Au relative to pp-collisions will show an enhancement at 
moderate p i. pe aking at p_L ~ 3 — 4 GeV before going down to unity 
at larger pj_ ' '. The argument is that energy loss in the "cold matter" 
of d-Au collisions is quite small and the pj_ broadening due to multiple 
scattering causes the Cronin effect. On the other hand, if no Cronin effect 
is seen and a suppression is seen instead it may be an indication that 
initial state effects are responsible for the phenomenon in both d-Au and 
Au-Au collisions. 

In this section, we will briefiy summarize recent work on proton- 
nucleus collisions in the CGC framework. Proton-Nucleus collision s jri 
this framework were first considered by Kovchegov and Mueller 
The proton in this case was modelled by a gauge invariant gluonic 
current. They showed that one obtains a Glauber- type formula with 
a saturation scale Qs (of the nucleus) for the inclusive gluon cross- 
section. The problem was looked at in more detail in Ref. The 
problem of proton-nucleus scattering was considered in a manner sim- 
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ilar to the classical fields treatment of nuclear scattering l-^^ l ^^ l except 
now they introduced two saturation scales Qsi and Qs2 for the nu- 
cleon and nucleus respectively with Q^i « Qs2- They were able to 
obtain analytical solutions for classical gluon production in the regions 
fc_L > Qs2 > Qsl as well as for Qs2 > fcj^ > Qsl- In the former, one ob- 
tains dN/d kj^d b oc QtiQi2/^± while in the latter case they obtained 
dN/d\d^h oc Ql-y\n{k\/Ql^)/k\. In this latter regime, one is solving 
the classical equations to all orders in 0^2/ No analytical solution 
is available for fcj^ < Qsl < Qs2, even when fc_L > f^QCD- The clas- 
sical problem discussed here can also be formulated numerically 
The A;_L dependence of the two kinematical r egio ns is seen clearly in the 
numerical result. It was also argued in Ref. ™"that since the rapidity 
distributions, for fixed k^, in one or the other regime are so different, 
these rapidity distributions could be used to isolate the Renormalization 
Group (RG)-evolution of the partons in both the proton and the large 
nucleus. 

The proton fragmentation region provides an excellent probe of the 
CGC. Final states measured in this region are produced by the scattering 
of high X partons in the proton off very small x partons in the nucleus. 
The scattering can be described as the convolution of the probability to 
find a quark in the proton q{xi , Q^) times the probability for the quark to 
scatter off the classical field of the nucle us ch aracterized by a saturation 
scale Qs2(a^2)^^''^- As discussed in Ref.^^, the kj_ distribution for this 
scattering is modified from the usual tree level pQCD distribution-this 
modified distribution, when convolved with the approp riate fra gmenta- 
tion functions will be refiected in hadronic final states H1II5DI Electro- 
magnetic final states such as photons and di-leptons are a particularly 
sensitive probe of saturation dynamics in the proton fragmentation re- 
gion EH. Interestingly, the modified distri butions lead to broad- 
ening of the final state and the Cronin eff ect This effect is also 
seen in other (dipole) models of saturation 1^^. 

Though the CGC picture gives rise to the Cronin effect (which can 
be represented by multiple scattering tree level diagrams), quantum ef- 
fects in the wavefunction can modify this picture significantly. In sec- 
tions 4.3 and 5.3 respectively, we discussed the CGC-|-geometric scaling 
mechanism which was first applied to understand the HERA DIS data 
and subsequently the suppression of p_L spectra of charged hadrons in 
Au-Au collisions at RHIC. In the geometric scaling regime, the anoma- 
lous dimensions for the evolution are very close to BFKL anomalous 
dimensions ^^-•^Ithis change in the anomalous dimensions (from the 
DGLAP one) is what causes the suppression of the pj^ spectra in the 
initial state CGC scena rio. A similar suppression must then persist in 
p/D-A collisions. In Ref.^^, it is predicted that the dependence of semi- 
hard processes on the number of participating nucleons of the nucleus 
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in D-A collisions will be ~ [^part)^^'^ ■ The number of participants are 
determined a la the Glauber calculation of Kharzeev and Nardil^. Thus 
more quantitatively, the suppression in D-Au collisions (relative to pp ) 
at ^/s — 200 GeV/nucleon due to the geometric scaling initial state effect 
translates into a 25-30% of moderately high p_i_ particles in the top 15% 
centrality range. This co nclusio n will be modified by the Cronin final 
state scatterings of Refs. I^^^these will change the behavior and one 
may even recover an enhancement. A quantitative calculation including 
both geometric scaling an d fin al state scattering is urgently required. 

Very recently, in Ref. qua ntitative predictions (along the lines 
of the qualitative ones in Ref. '-^ were made for hadron multiplicities 
in Deuteron-Gold scattering at RHIC-in particular for the rapidity and 
centrality dependence. The results of the current Deuteron-Gold run 
should therefore help determine whether a Color Glass Condensate is 
formed already at the moderately high energies at RHIC. 

An important part of the RHIC program is the study of peripheral 
nuclear collisions. Very intense electromagnetic fields are created at these 
energies and a variet y of photon-photon and photon-Pomeron final states 
can be studied Of particular interest to us is inclusive and diffrac- 
tive QQ-in the color field of a nucleus xhe problem is analogous to 
the photo-production of heavy quark pairs in Deeply Inelastic Scatter- 
ing. The transverse momentum and invariant mass distribution of quark 
pairs can be shown to depend sensitively on the saturation scale Qs. 
Thus ultra-peripheral nuclear collisions provide an independent method 
to extract properties of Color Glass Condensate. 
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